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MONTHLY NOTICES 


OF THE 


ROYAL ASTRONOMICAL SOCIETY 


Vol. 118 No. 5 





MEETING OF 1958 MARCH 14 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 
John Anderson, 22 Massey Park, Stormont, Belfast, N. Ireland (proposed 
by E. M. Lindsay); 
Charles David Austin, 49 Leysdown Road, London, S.E.g (proposed by 
C. W. Allen); 
Frederick Henry George Best, 24 Blake Road, Croydon, Surrey (proposed 
by P. Moore); 
Michael Gadsden, I.G.Y. Station, Invercargill, New Zealand (proposed by 
I. L. Thomsen); 
Eugene C. Robertson, Geological Survey, Silver Spring, Maryland, U.S.A. 
(proposed by A. H. Cook); 
Malcolm Walsh, ro Penylan Road, Newport, Mon. (proposed by C. W. Allen); 
and 
Jack Youdale, 21 Stanhope Road, Billingham, Co. Durham (proposed by 
D. Sinden). 


The re-election by the Council of the following Fellow was duly confirmed :— 
Victor Dumert, 88 Church Lane, London, N.2 (proposed by H. Bondi). 


The election by the Council of the following Junior Members was duly 


confirmed :— 








Robert John Dickens, Royal Greenwich Observatory, Herstmonceux Castle, 
Hailsham, Sussex (proposed by B. E. J. Pagel); 

Roger Anthony Gordon Smith, 184 King’s Hall Road, Beckenham, Kent 
(proposed by B. E. J. Pagel); and 

Ramon David Wolstencroft, 138 Butterstile Lane, Prestwich, Lancs. (pro- 
posed by C. W. Allen). 


Seventy-eight presents were announced as having been received since the 

last meeting, including :— 

T. Maloney: Other worlds in space (presented by John Calder, Publishers, 
Ltd.). 
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MEETING OF 1958 APRIL 11 
Dr W. H. Steavenson, President, in the Chair 


The President announced that the Council had elected the following 

Associates of the Society :— 

Francis Birch, Dunbar Laboratory, Harvard University, Cambridge 38, 
Mass., U.S.A.; 

Boris Vassilevitch Kukarkin, Chief of the Variable Star Department of the 
Sternberg Astronomical Institute, Moscow, U.S.S.R.; and 

Yngve Ohman, Director, Stockholms Observatorium, Saltsjébaden, Sweden. 


The election by the Council of the following Fellows was duly confirmed :— 

Donald Thomas Germain-Jones, British Petroleum Company, Beaufort 
House, Gravel Lane, London, E.1 (proposed by L. H. Tarrant); 

Alexander Leo Helm, Beech House, Holmewood, Langton Green, Tunbridge 
Wells, Kent (proposed by P. Moore); 

Ernest George Hill, Harbour St. Bride, Durlston Road, Swanage, Dorset 
(proposed by H. G. Hughes); 

R. Ingram-Brown, Brown, Son & Ferguson, Ltd., 52/58 Darnley Street, 
Pollockshields, Glasgow (proposed by C. W. 'T. Layton); 

Roger Jackson Lee, Box 2391, Brown Station, Providence 12, Rhode 
Island, U.S.A. (proposed by L. M. Milne-Thomson); 

Martin F. McCarthy, Specola Vaticana, Palazzo Pontificio, Castel Gandolfo, 
Italy (proposed by W. P. Bidelman); and 

John Grahame Sutherland, 172 Overhill Road, London, S.E.22 (proposed 
by H. J. Lewis). 


Fifty-eight presents were announced as having been received since the last 


meeting, including : 


E. O. Tancock (ed.), Philip’s Chart of the Stars (presented by G. Philip & 
Son); 

F. B. Wood (ed.), The present and the future of the telescope of moderate size 
(presented by the University of Pennsylvania Press); 

N. Roman (ed.), /.4.U. Symposium No. 5. The large-scale structure of the 
galactic system (presented by the International Astronomical Union); and 

M. Migeotte, L. Neven, J. Swensson and E. Vigroux, An atlas of nitrous 
oxide, methane and ozone infra-red absorption bands (presented by the 
Institut d’Astrophysique de |’ Université de Liége). 


MEETING OF 1958 MAY g 
Dr W. H. Steavenson, President, in the Chair 


The President announced the death of Knut Emil Lundmark, an Associate 


of the Society, and paid a tribute to his memory, the Fellows standing. 


The election by the Council of the following Fellows was duly confirmed :— 


Donald Murray Edwards, c/o Mapping Branch, Lands and Surveys, Perth, 
W. Australia (proposed by C. R, Edwards); 








cor 


si 











No. 5, 1958 





Meeting of 1958 May 9 409 





Hubert Basil Lewis, 8 Baker Street, Moorabbin, Australia (proposed by 
E. B. Walton); 

Donald Malcolm, 18 Caledonia Street, Paisley, Scotland (proposed by A. E. 
Roy); 

Seyed Mohamed Razavy, 51a Blenheim Crescent, London, W.11 (proposed 
by H. Wildey); 

Thomas E. Sadler, Gravenhunger, Woore, Staffs. (proposed by C. D. 
Grimwade); 

Cyril Alfred Edgar 'l‘omes, 65 East Dulwich Road, London, S.E.22 (proposed 
by P. Moore); 

Ernest W. Wilkes, Allendale Road, Hexham (proposed by C. Gilbert); 

William Williamson, 1 Gold Street, Longton, Staffs. (proposed by C. D. 
Grimwade); and 

Ronald William Girdler, Corpus Christi College, Cambridge (proposed by 
R. Stoneley). 


The re-election by the Council of the following Fellow was duly confirmed :— 


Herbert Gilgryst, 5 De Vere Gardens, London, W.8 (proposed by C. M. 
Botley). 


The election by the Council of the following Junior Member was duly 
confirmed :— 


Lionel William Moore, 190 Pickhurst Lane, Hayes, Kent (proposed by 
C. W. Allen). 


ADDITIONAL MEETING OF 1958 JULY 11 


in the Physics Lecture Theatre of the University College of North 
Staffordshire, Keele 


Dr W. H. Steavenson, President, in the Chair 


The President announced that this Additional Meeting was the ninth to 
take place outside the rooms of the Society in London. He expressed the 
thanks of the Society to Sir George Barnes, Principal of the University College 
of North Staffordshire, and to the College authorities, at whose invitation the 
Meeting was being held at Keele. ‘The thanks of the Society were also extended 
to Professor H. D. Springall and Professor F. A. Vick for the use of the Depart- 
ments of Chemistry and Physics respectively, to the Registrar, Mr. J. F. N. 
Hodgkinson, for his co-operation in making the necessary arrangements for the 
Society’s stay, and to Dr J. N. Hodgson who, in acting as the Society’s Local 
Secretary, had undertaken much of the detailed work involved in the local 
organization of the various meetings and associated activities. 


One hundred and thirteen presents were announced as having been received 
since the last meeting, including :— 

G. Alter, J. Ruprecht and V. Vanysek, Catalogue of star clusters and associa- 
tions (presented by G. Alter); 
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A. Dauvillier, Le volcanisme lunaire et terrestre (presented by Editions 
Albin Michel); 


G. de Vaucouleurs, La photographie astronomique (presented by Editions 
Albin Michel); and 

United Nations Educational, Scientific and Cultural Organization, The 
Earth as a Planet (presented by U.N.E.S.C.O.) 


Other meetings and activities in connection with the Society’s visit to Keele 
included a colloquium on Interplanetary Matter, held in the morning of July to. 
In the afternoon of the same day, at the kind invitation of the Director, Professor 
A. C. B. Lovell, and with the co-operation of his staff, there was a conducted 
tour of the Jodrell Bank Experimental Station. In the evening a Public Lecture 
entitled Why is it dark at night ? was delivered in the College Conference Hall 
by Professor H. Bondi. On July 11 a number of excursions were arranged to 
places of scenic and historical interest in the neighbourhood. 
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THE CONTRIBUTION OF THE IMPERSONAL ASTROLABE 
TO FUNDAMENTAL ASTRONOMY 


George Darwin Lecture delivered by Professor André Danjon on 1958 May 9 


I wish to recall that Sir George Darwin, to whom this lecture is dedicated, 
devoted a part of his work to rotating celestial bodies. The title of his first 
paper was On the Influence of Geological Changes on the Earth’s Axis of Rotation. 
Since then the discovery of the movements of the poles and, later, of the variations 
in the length of the day have aroused a renewed interest in the problems of 
practical astronomy related to these phenomena. My aim in this lecture is to 
describe the programme of observations now being carried out at the Paris 
Observatory and to mention the further developments that can be anticipated 
in the near future. 

A few years before the outbreak of the last war, I had published the description 
of a new type of transit instrument for the determination of right ascensions. 
The results obtained at Strasbourg with a prototype of small dimensions (aperture 
6 cm) having been very encouraging, a new series of tests was commenced at the 
Paris Observatory in 1946. Meanwhile, an instrument of similar design but of 
12 cm aperture was constructed. The early observations made with this larger 
instrument revealed unforeseen difficulties, which could have been overcome 
by completely redesigning the mounting. But, for the time being, the study | 
of this instrument has been discontinued and all the facilities available in Paris 
Observatory, including scientific and technical staff, laboratories and workshop, 
etc. have been concentrated on the construction of an instrument of an entirely 
different type: the impersonal astrolabe, which will be the main topic of this 
lecture. 

I had moreover another reason for giving priority to the astrolabe, of which 
an experimental model had already been constructed at the Observatory in 1950 
and progressively improved. This reason is that, in 1952, the International 
Council of Scientific Unions decided to organize the International Geophysical 
Year and to incorporate in its programme a determination of world longitudes 
and latitudes. The observatories participating in the operation were recom- 
mended to observe stars selected in the Fundamental Catalogue FK3 (or in the 
improved catalogue FK3R). It was agreed that preference should be given 
to modern instruments giving simultaneously and with the highest possible 
accuracy the two coordinates, time and latitude. ‘The zenith telescope, indeed, 
gives both coordinates but requires a special catalogue limited to zenith stars, 
in which therefore non-fundamental stars have to be introduced. ‘The astrolabe 
is the only instrument which permits both coordinates to be determined by the 
observation of fundamental stars alone. 

A large size astrolabe was then ordered from the Société Optique et Précision 
de Levallois (OPL). The design of the instrument was worked out in close 
co-operation with the technicians of this firm. The new instrument differs 
from the experimental model mainly in its dimensions, which are larger, and in the 
improved arrangement of the mechanical parts. 
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The instrument was delivered in 1956 and put into service in July of that 
year. Since then it has been in use without interruption. More than twelve 
thousand transits have been observed, which is adequate for evaluating the 
qualities of the instrument, if not for establishing definitive results. The 
observations have been carried out mainly by L. Arbey, B. Guinot and Mlle 
Débarbat. 

Under the most favourable atmospheric conditions, the OPL astrolabe 
permits the transit of a star, at a well determined zenith distance, differing very 
slightly from thirty degrees, to be recorded with an accuracy characterized by 
a standard deviation of 0”:og, assuming that the errors of the catalogue have been 
corrected. ‘The deviation increases with atmospheric scintillation, and_ its 
most frequent value is 0”-17._ For only six per cent of the nights is it larger than 
0-26. But what is most important is that the observations have proved to be 
free from systematic errors, a result of outstanding importance for the use of 
the impersonal astrolabe as a fundamental instrument. 
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Indeed we very soon realized that this instrument was useful not only for the 
determination of time and latitude but that it could be used as well for the 
improvement of the FK3 catalogue and for the determination of the fundamental 
astronomical constants. The accuracy of the results is due partly to the intrinsic 
qualities of the instrument and partly to the application of the so-called ‘‘ method 
of equal altitude of stars’’. 

As for the instrument itself I wish first to recall its principle. Suppose that 
we fix an equilateral glass prism in front of a horizontal telescope in such a 
way that its edges are horizontal and one of its faces vertical (Fig. 1). A beam 
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of light coming from a star entering the face AB at right angles to the edge of the 
prism, is then reflected from AC and emerges through the lower half of BC. 
Another portion of the beam falls on AC after having been reflected from the 
surface of a bath of mercury. It is reflected by AB and emerges through the 
upper half of BC. Thus, looking through the telescope, two images of the star, 
S and S’, are seen which appear to move in opposite directions, one moving 
upwards and the other downwards. ‘They appear to coincide when the star is at 
a zenith distance of exactly 30° and the instant of coincidence is noted. As the 
prism angle is not exactly of 60°, there are three unknowns in the problem, 
namely: the correction of the clock, the latitude, and the excess of the zenith 
distance over 30°. ‘Therefore at least three stars, in different azimuths, must 
be observed. In fact, many more than three stars are observed, and in general 
each of our groups contains twenty-eight stars. The plotting of position lines, 
a method familiar to navigators, gives an approximate graphical solution but 
as our observations are very accurate a solution by least squares should be applied. 

The astrolabe in its original form was designed by MM. Claude et Driencourt*. 
It dates back to 1g00 and it was then a very simple instrument which consisted 
merely of a telescope, a prism, a mercury horizon and adjustment devices. The 
observer noted the time of coincidence of the two images of the star by the eye 
and ear method or by the hand tapping method. ‘Thus, each star gave only a 
single determination of time, which was affected by personal equation. But 
it was intended to be an instrument for use in the field and not an instrument of 
high precision. ‘There was indeed another source of error resulting from the 
images being displaced when the position of the eye-piece was changed even 
by a slight amount. The reason is that the images are formed by two pencils 
of rays the axes of which are not strictly parallel but are inclined to each other 
at an angle of two or three degrees. If the focusing is altered, the distance 
between the two pseudo-images changes, which produces the same effect as a 
variation inthe prism angle. [| will later show you that this angle has an astonishing 
stability, but, owing to the defective arrangement of the instrument, no full 
advantage of this remarkable property was taken. 

Various astronomers have been successful either in correcting this optical 
defect or in rendering the astrolabe an impersonal instrument. ‘The excellent 
results obtained by Paul Muller, the well known double star observer, with a 
micrometer containing a doubly refracting prism, led me to the concept of 
a more complete solution which I will now describe and which is the essential 
feature of the impersonal astrolabet. 

A doubly-refracting prism is placed at the focus of the telescope. For 
certain reasons, which would take too long to describe, I have selected a double 
Wollaston prism made of quartz. The arrangement of this prism is shown in 
Fig. 2. The separation angle of this prism is precisely equal to the angle between 
the axes a and 5 of the two pencils transmitted by the objective. The splitting 
of the pencil a gives two pencils a, and a,, as is shown in the figure, while the 
splitting of pencil b gives pencils 6, and 6. It is at once seen that a, and b, are 

* A. Claude and L. Driencourt, Description et usage de Il’ Astrolabe a Prisme, Paris, Gauthier- 
Villars, 1910. 

J. Ball and Knox Shaw, A Handbook of the Prismatic Astrolabe, Cairo, 1919. 

See also: E. Chandon and A. Gougenheim, Les Jnstruments pour l’ observation des Hautcurs Egales 


en Astronomie, Revue Hydrographique, 12, 1, 1935- 
+ A. Danjon, L’Astrolabe impersonnel de l’Observatoire de Paris, Bull. Astro., 18, 251, 1955. 
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parallel and that therefore the distance of the two corresponding images o and 
o’ no longer depends on the adjustment of the eyepiece. The two other images 
are eliminated by suitably placed screens. Moreover, the coincidence of the 
two images o and o’ can be obtained by a displacement of the quartz prism parallel 
to the axis of the telescope. The coincidence thus obtained can be maintained 














Fic. 2. 


by traversing the prism at a rate proportional to the sine of the azimuth and to 
the cosine of the latitude. The micrometer screw which drives the prism, 
carries a disk with contacts connected to a chronograph. ‘The astrolabe, already 
corrected for its optical defects, then becomes impersonal. ‘The observer 
needs only to correct the movement produced by a motor, so as to maintain the 
coincidence of the two images. 

In fact this coincidence is obtained not by superposing the two images, but 
by placing them side by side along the same horizontal line between two illuminated 
wires as shown in Fig. 3. 














Fic. 3. 


It would take too long to describe, in all details, the arrangements adopted 
for the component parts of the astrolabe*. Figure 4 gives a very simplified layout 
of the instrument. The objective of the OPL model has an aperture of ten 
centimetres and a focal distance of one meter. ‘Two equal lune-shaped holes 
are made in a diaphragm, giving diffraction images elongated along the vertical, 
with axes in the ratio of three to five. The faces of the prism are rectangles of 
10cm by8cm. The dimensions of the micrometer screw are rather unusual: 
25 mm in diameter and 5 mm in pitch. If the star is observed exactly to the 
East or to the West, the screw makes a complete revolution in 2°61 seconds. The 
speed of the motor is constant and is adjusted to suit the latitude of the station. 


* A. Danjon, “ L’Astrolabe impersonnel modéle OPL,’’, Bull. Astro., to be published in 1958. 
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A speed reducer automatically varies the speed of rotation of the screw pro- 
portionally to the sine of the azimuth, when the astrolabe is turned around its 
vertical axis to observe another star. To make the instrument as compact as 
possible, the beam of light is broken twice by mirrors made of fused quartz 
and located between the objective and its focus. 





Fic. 4. 


The accuracy of astrolabe observations is due to a large extent to careful 
manufacture by the Société OPL but equally to the working conditions and to the 
intrinsic properties of the equal altitude method. I shall consider these two 
points in turn. As the instrument rotates around a vertical axis, its flexures are 
not liable to vary with the azimuth. Its design provides protection against 
sudden changes of temperature. The results of time and latitude observations 
do not even depend on the absolute value of the prism angle; they could be 
affected only by a variation of this angle during the observation. But results 
obtained since July 1956 have clearly demonstrated that there is no systematic 
variation of this angle in the course of the night, not even in the course of a year, 
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and that the zenith distance remains invariable within a few hundredths of a 
second. Our observation programme includes twelve groups, each of twenty- 
eight stars chosen from the FK3 or its supplements. These groups are identified 
by letters a, b,..., k, 1. 

A single observer can observe, during a single night, at least two consecutive 
groups, often three, or even more*. ‘The computed zenith distance, corrected for 
refraction, would have the same value for all these groups if the prism angle were 
constant and if the catalogue were not affected by any systematic error. In fact, 
two consecutive groups observed in the course of the same night, give different 
values as is shown by Fig. 5. In the course of a year, differences a—b, b-c, 

., La are obtained. The distribution of these differences does not reveal 
any systematic variation in the nature of a function of the right ascension, which 
indicates that the prism angle is independent of the annual variation of temperature 
and that errors in the catalogue are solely responsible for the differences obtained. 


“ON AA Wh 
ce VT 


Fic. 5 


+010 





Even more, we have been able to verify that the zenith distance of transits does 
not vary systematically during the course of the night, although, as a general 
rule, the temperature of the instrument decreases from the evening to the morning. 
This is illustrated by the smallness of the closing error of the differences that 
B. Guinot has computed for a complete year. For instance, for group k this 
closing error was found as small as —o”-017+0”-104. At any rate it is not 
significant. It must be emphasized that the prism has been very carefully 
annealed, that its mounting does not exert any prejudicial stresses, and that these 
two conditions seem to be essential. 

Let us now consider the method of equal altitude of stars and its advantages. 
Devised by Gauss about a hundred and fifty years ago for field operations, its 
application to problems of fundamental astronomy has so far been restricted 
to special cases such as, for instance, the determination of latitude by the Talcott 
method, which is still in use in the International Latitude Service, or the observa- 
tion of transits in the prime vertical by the method of W. Struve. However, it 
is very valuable, in its general form, for linking together stars whose 
declinations differ by several tens of degrees and whose right ascensions differ 
by several hours. 

The declinations of the stars that can be observed with the astrolabe are 
between ¢—30° and ¢+30°, ¢ being the latitude of the station. In 
Paris, these two limits are theoretically + 18°50’ and + 78°50’, or practically 
+20° and +78°. Their hour angle at the time of the transit is comprised 
between the limits —H and +H, H being the solution of the equation 
sin H=}sec¢. For Paris 2H=6" 36™ and for Herstmonceux, 6" 59™. ‘The 
observation of one group requires about one hour and a half so that it is possible 
to link together stars whose right ascensions differ by as much as 8 hours. 


* B. Guinot, “ L’Astrolabe Impersonnel. Réduction des Observations. Etude des Résultats.”’ 
Bull. Astro., 18, 283, 1956. 
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When two transits of the same star are observed, one to the east and the 
the other to the west, the mean of the observed times gives its right ascension, 
and its declination is easily computed from the difference of the same transit 
times. ‘The determination of the right ascension is not very accurate when the 
declination of the star is near the limits and, on the other hand, the determination 
of the declination is not possible when the hour angle of the transit is close to H. 
To fill these gaps it would be advantageous, and for more than one reason, to 
set up astrolabe stations in several latitudes. In Paris for instance we obtain 
good values of the right ascensions for declinations comprised between + 20° 
and + 78°, and good values for the declinations in the intervals between + 20° 
and +50° and between +68° and +78°. In these intervals the application 
of the method of equal altitude with a high precision instrument discloses at 
once the indvidual errors of the catalogue, as these errors are important compared 
to the accuracy of the observations. I have already mentioned that the average 
value of the standard deviation for one transit was 0”-17 and that, under very 
favourable atmospheric conditions, this deviation is as small as 0”-og. As each 
group is observed about 20 times each year, the plotting of the position lines is 
correct within a few hundredths of a second. Ina few years the precision of their 
determination will be as accurate as 0”-o1. Individual errors being expressed 
by a dispersion of the position lines in the diagram, some of them intersect the 
circle, whose centre and radius are obtained from a least squares adjustment, 
and some others remain outside the circle. As for the systematic errors of the 
catalogue they are sometimes revealed by an asymmetry of the diagram. 


D 


4850 














459 50 571 aiid 4°9 S°O 5% *” 
Fic. 6 
AT= 48:987 + 0°-006 AT= 45987 + 08-002 
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The diagram in Fig. 6(a) was drawn from the results obtained on 1958 January 
22 ongroupd. The same group was observed in 1957 and the displacements to 
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be given to the position lines to make them tangent to the circle were already 
known. Fig. 6(b) shows the same diagram corrected from the observations 
obtained in the preceding year. Their quadratic mean error decreases from 
o”-21 to 0”-og and the asymmetry disappears. This example illustrates clearly 
the effectiveness of the method of equal altitudes in correcting the position of 
stars, but it must be remarked that this operation is merely an internal adjustment 
affecting each group independently and that it is by no means sufficient to correct 
the catalogue. 

It is now essential to make a fundamental remark. Whatever be the instru- 
ment used, a single observation of a group of stars cannot disclose systematic 
errors of a catalogue which are equivalent to a change in the coordinate system 
and which do not affect the relative positions of the stars. Such are those produced 
by an error in the value of the precession or of the nutation. As the relative 
positions of the stars remain unchanged but the position of the pole and the 
origin of the coordinate system are erroneous, the diagram of the position lines 
is translated without distortion and therefore neither its symmetry nor the radius 
of the circle of the position lines is affected. But of course the coordinates of 
the centre of the circle are altered, and consequently so also are the time and 
the latitude. The method of equal altitude of stars, which reveals so clearly the 
individual errors and the localised errors affecting the various regions of the 
sky, cannot therefore provide any evidence of a systematic error in the reference 
system by a single observation of a given group. Errors of this sort can be found 
out only by a chain programme of observation, associated with group to group 
connections, covering at least one entire year. These errors will then appear in 
the group errors in the form of a variation depending on sin «. 

I have said a moment ago that there was in this respect no privileged instru- 
ment. Suppose that we observe with a meridian circle the transits of a group 
of stars of limited range in right ascension. If the pole of the catalogue is at a 
distance Av from the true pole of rotation, in consequence of an error in the 
precession, and if am is the observed mean right ascension of the group, the 
azimuth of the meridian circle, computed from the observations, will have an 
error equal to Av sin am sec ¢, but nevertheless the representation of the observa- 
tions by means of this faulty azimuth will still be excellent. Therefore no 
improvement of the catalogue can be expected. Here again the ‘‘ chain adjust- 
ment ’’ alone will enable the errors of the catalogueto be determined and 
corrected. 

Theoretically the method of equal altitude of stars and the method of meridian 
transits are therefore equivalent in this respect. But in practice their possibilities 
are very different. In the first place, the stability of the astrolabe can be 
guaranteed with much greater certainty than that of the meridian circle. In 
particular, the astrolabe enables the value of the prism angle to be accurately 
controlled, and that is the only instrumental constant whose variation would be 
prejudicial to the accuracy of the observations. We now know by experience 
that the changes in the prism angle are negligible compared with the accidental 
or systematic errors of the fundamental catalogue. Moreover, the observations 
with the astrolabe are referred to a single terrestrial reference: the vertical. 
The discussion of the latitude observations have proved, as I will show you, 
that the daily variations of the zenithal refractions are imperceptible, at least in 
Paris. 
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The conditions are quite different in the transit instrument. The solid 
body which serves as the direction of reference is not now a simple glass prism 
but a telescope with its trunnions. The line of sight is determined by the optical 
centre of an objective, consisting of two lenses which cannot be linked rigidly 
together nor to their mounting, and by the spider’s web carried by the movable 
frame of a micrometer. ‘The objective and the micrometer are linked together 
by a metallic tube exposed to unceasingly variable conditions of temperature, 
as this tube radiates either from one side or the other as the star to be observed 
is to the north or to the south. The axis of rotation is defined by two trunnions 
carried by a heavy metallic part exposed to the same variations of temperature. 
Consequently the position of the instrument changes with respect to the vertical 
from one star to another, and as the flexures are not rigorously symmetrical with 
respect to the meridian plane, lateral deflections are unavoidable. To measure 
the collimation of the instrument and the inclination of its axis of rotation, the 
telescope is pointed towards the nadir in a direction quite different from the 
positions it occupies for the observation of stars. Who could say what happens 
to this collimation and inclination when the telescope is pointed towards a star ? 
To obtain an accuracy of 1/10th of a second of arc, it is necessary that the mech- 
anical or thermal deformations of the instrument are less than 1 micron. 
Experience shows that they are in general much larger than this limit. More 
than 70 years ago, Chandler already asserted that the method of equal altitudes 
of stars was better suited to give accurate results than meridian astronomy*. 
We have reached this very conclusion by comparing the physical properties of 
the astrolabe and the transit instrument. But there is another reason, of a 
geometrical nature, which confirms this conclusion. With the method of equal 
altitudes, transits of stars can be observed all around a complete circle 
(almucantar) on the celestial sphere ; and that is why the radius of the circle and 
its pole, which is the zenith of the station, are so well determined. Meridian 
observations provide different conditions, as the observations can only be made 
over 45° of zenith distance, that is along an arc which covers only one quarter 
of the circumference on the meridian; thus, the position of the pole of this arc, 
which defines the direction of the rotation axis of the transit instrument, is 
practically less accurately determined than that of the zenith of the astrolabe. 

Even if it were perfect the meridian circle by itself could not conveniently 
provide for the ‘‘ chain adjustment ’’, because such adjustment cannot be 
effected without referring, as an intermediate step, to an azimuth mark or a 
collimator. I know nothing more disappointing than the observation of a mark. 
I have often observed those of the Strasbourg Observatory which were set up 
very carefully by Winnecke. Even when at first sight they appeared to be stable, 
they suffered every now and then from sudden variations amounting frequently 
to one second of arc. If the mark is located not far from the instrument then 
its azimuth is not well defined. If it is far from the instrument its direction 1s 

* §. C. Chandler, ‘‘ The Almucantar’’, Ann. of Harvard Coll. Obs., 17, 1887, 179-181. 

‘* T think it can be successfully maintained that the almucantar system of observation constitutes 
a more simple and direct solution of the problem of determining the absolute time from the pheno- 
mena of diurnal motion, and consequently one more likely to be free from the systematic error, than 
that of the meridian transit.’’.... ‘“‘ As regards the adaptation of the instrument to purely astro- 
nomical purposes, one line of usefulness seems to be clearly marked out, in the application to the 


class of position measurements for which meridian instruments have heretofore been exclusively 
employed,”’ 
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frequently deviated by lateral refractions. But in any case a terrestrial mark is 
far from being as reliable as the prism of the astrolabe for defining a direction 
of reference. Its stability cannot even be assure: through the course of a night, 
from evening to morning, and far less during day time. 

In my opinion, better results can be expected from a collimator than from a 
mark, because the collimator can be installed in the interior of the meridian 
room and can then be efficiently protected against variations of temperature, as 
well as its pillar and even the ground all around the pillar. That is why I have 
used two collimators, one to the east and the other to the west, with the transit 
instrument that I mentioned at the beginning of this lecture and of which I will 
now recall briefly the principle (Fig. 7). In front of the objective of a telescope 


whose axis lies in the meridian plane, is mounted a reversing prism, which covers 
its surface with the exception of two segments which remain uncovered. In the 
field of the telescope two images of the star appear, which move in opposite 
directions. ‘The two images coincide when the star transits the plane of the 
reflecting face, which differs but by a very small amount from the direction of 
the meridian plane. ‘The position of the prism is then determined by auto- 
collimation, by means of the lateral collimator pointed towards its reflecting 
face. The readings made with this collimator provide the necessary correction 
for reducing all the transits to a same plane. ‘The second collimator is used 
when the instrument is reversed. I will not mention the various operations 
which enable the horizontality of the collimator to be controlled by means of a 
mercury horizon and its azimuth to be controlled by the observation of 
fundamental stars. What I wish to stress is that the line of sight is much better 
defined by the prism than by the optical axis of the conventional meridian 
telescope. The major difficulties occurring in the meridian telescope are 
automatically eliminated, but indeed other problems arise which I am now 
endeavouring to overcome. If there is a noticeable difference in temperature 
between the interior and the exterior of the room, the reflected image of the 
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cross wires as seen in the collimator is unsteady. This is but a minor 
inconvenience which can be easily corrected. The most troublesome effect is 
caused by the position of the prism in front of the objective. If its temperature 
differs from that of the surrounding air even by a few tenths of a degree, it is 
enveloped by a layer of air about 1 cm thick within which there exists a marked 
temperature gradient. The part of the pencil of rays which crosses this layer 
becomes divergent if the prism is warmer than the air and convergent if it is 
cooler, and in both cases the image given by the free aperture is altered. If the 
air around the prism is stirred by means of a fan, this layer of air disappears and 
the image resumes its normal aspect. But if the fan is stopped, then this layer 
is at once formed again. 

A similar effect is also produced in instruments in which metallic parts are 
placed across the path of the light. I have observed for instance, by a very 
sensitive interferometer method, that the bright extensions which surround the 
stellar images in a reflecting telescope are caused by the presence of the supports 
of the secondary mirror which, cooling off during the night by radiation, are 
surrounded by a layer of warmer air producing a refraction effect, and not a 
diffraction effect as is generally assumed. André Couder gets rid of this effect 
by covering the supports of the secondary mirror with sheets of polished 
aluminium. 

The case of the reversing prism of my transit instrument cannot be treated 
so easily and I have had to adopt a more radical solution : the prism of the 
projected instrument will be placed in an atmosphere of helium, this gas having 
a refractivity ten times smaller than that of air and a conductivity six times 
higher. The idea is not a new one : Bernard Lyot had elaborated plans for a 
reflecting telescope of an aperture of 150 cm without a dome, to be installed at 
the Observatory of the Pic du Midi, the tube of which, sealed by a glass plate 
at its upper part, would have been filled with helium. 

I have given these detailed explanations in order to recall, in a specific case, 
how different can be the actual properties of an instrument and its theoretical 
properties. When an inventor forgets this important truth, the actual use of 
the instrument soon reminds him of it. But let us come back to our subject 
and take up again the ‘‘ chain adjustment ’’. Years will elapse before we know 
if this new instrument is as well adapted as the astrolabe for this particular 
purpose. 

I will not wait any longer to answer objections raised by various authors 
about this chain programme, the legitimacy of which has been disputed. ‘Their 
criticisms are based, in particular, on the closing errors, sometimes quite 
large, found for certain latitude stations equipped with zenith telescopes*. 
The mean value of the closing errors for the years 1956 and 1957, for 
the four stations of the International Latitude Service, was —0”:124 between 
the evening group and the middle group and —0”-101 between the middle group 
and the morning group. But the individual values show a very large scattering. 
The extreme annual means were —0”:345 for Ukiah and +0"-162 for Kitab. 
If these closing errors had the same value for all the stations they could 
be eliminated by correcting the aberration constant but they are very different 
from one station to the other. ‘This dispersion can be assigned only to 


*G. Cecchini, Relazione sul attivita del Servizio Internazionale delle Latitudine, 1.U.G.G, 
General Assembly, Toronto, 1957. 
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instrumental or to meteorological factors. I am myself convinced that the first 
occur much more frequently and are in general more harmful. 

It may also occur that the local terrain or the presence of large buildings in 
the immediate vicinity of the instrument will give rise during the night to atmos- 
pheric turbulence which produces appreciable accidental refractions. Such 
cases are no doubt exceptional as the site of an observatory is as a general rule 
selected with full knowledge of these facts. But a fall in the temperature in the 
course of the night affects the level and the micrometer screw of the zenith 
telescope, and these effects contribute to the closing error of the chain programme. 

The astrolabe does not have such disadvantages. Dr Guinot, who has 
determined the connections of our groups, has found for the first year of observa- 
tion a very small closing error. He first smoothed the observations, as I have 
already stated for the group of January 22, in order to free them from individual 
or local errors. Then he computed the differences between consecutive groups 
observed by the same observer in the course of the same night. He started with 
group k and came back to the same group after a complete year. Table I gives 
the differences in time and latitude with the related standard deviations. 


TaBLe I 
1956-57 
Groups Time Latitude 
unit : o%-oo1 Unit: 0”-o1 n 
j-k +15'0 $4°1 —2°9 ae 5 
k-l + 3:2 2°2 +0°5 4°0 12 
l-a — 80 hy | +3°0 3°6 10 
a-b —10°0 2°5 —3°1 2°9 II 
b-c +23°0 7 +5°7 18 10 
c-—d — 12°! 1 —38 2°0 12 
d-e + 3°0 ry +1°9 2°0 13 
e-f + O5 2 —26 1°2 13 
f-g — O'5 I°5 —7°3 16 14 
g-h —12°8 I's +7°4 2s 13 
h-i + 1°97 18 —4'9 1°4 15 
i-j — 16 I's +2°6 2°5 16 
Closing error + 1°4 say —3°5 +8°7 


These results have been confirmed by those of the second year of observations 
which will end by July next. 

Certain differences, mainly in right ascensions, reach surprisingly high 
values. In the course of the same night, groups b and c give clock-corrections 
differing by an amount of 23 milliseconds. But what is still more remarkable is 
the smallness of the closing errors which are not significant. We can therefore 
assume that the chain method of adjustment is permissible, and we can conclude 
that not only the astrolabe but also the atmosphere of Paris is very stable. If, 
contrary to my opinion, the closing errors observed elsewhere were not chiefly 
due to instrumental causes, we would then have to admit that, by bad fortune, 
these other stations had all been located in exceptionally unfavourable sites and 
that Paris enjoyed uniquely privileged conditions. Flattering though this last 
assumption may be for the pride of the French, I strongly refuse to believe in 
the poor quality of the other latitude stations and I persist in accusing the zenith 
telescopes, the control of which seems to me very difficult. By control I mean a 
permanent and absolute control such as that which is provided for the astrolabe, 
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at least twice per night, by the determination of the prism angle and not a mere 
testing of the instrument at long intervals. 

Having conclusively established that the application of group corrections is 
permissible, Dr Guinot computes the corrections required to refer the 
determinations of time, latitude and zenith distance to a provisional system 
more consistent than that of the FK3. 


TasB_e II 
a b c d e f g h i j k | 
Time +o'1 —10°0 +1299 +0°7 +36 +3°9 +3°3 —9°6 —80 —9°7 +5:°1 +82 


(08-001) 

Latitude +1°5 — 1°33 + 4°7 +12 +36 +12 —5°8 +2°0 —2:11 O10 —26 —1°8 
(o”-o1) 

Zen. Dist. —7°3 + 1°4 + 2°44 +08 —7°4 +2°3 +18 47°4 —3°2 +2°4. +2°9 —3°5 
(o”-o1) 


Fig. 8 represents the group corrections for time and latitude reduced to 
nearly the same scale, as an arc of 1 millisecond at the zenith in Paris is nearly 
equal to o”-o1. It appears at first sight that the right ascensions of the FK3 are 
less homogeneous than the declinations. These differences may be partly 
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explained by the uncertainty of the proper motions, the computation of which 
depends to some extent upon meridian observations made before the introduction 
of the travelling wire micrometer. But that such deficiencies still exist in the 
FK3 is a proof of the inability of conventional methods to link together stars 
with very different right ascensions. 

The analysis of group corrections discloses an irregularity of small amplitude 
in right ascension depending upon sin «. ‘The corresponding error in the declina- 
tions of the form cos « is still smaller. The pole of the provisional system, which 
practically coincides with that of the FK3, therefore differs but little from the 


29 
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true pole. But the systematic errors of the FK3 are much smaller than the local 
errors, which concern areas of an extent comparable to that of the constellations. 

We shall know in a few years whether the corrections given in Table II are 
stable or whether they vary in course of time. In the latter case their variations 
will enable corrections to the constants both of precession and nutation to be 
derived. But we shall have to wait for 19 years until we can accurately sort out 
these two constants. 

It is well known that the constant of aberration k can easily be deduced from 
the closing error in latitude, but the value obtained from the first year of 
observations, equal to 20”-48, has only a low weight. B. Guinot has shown that 
the time determinations could also be used for this purpose. Group d can be 
observed in Paris for about six months. If the conventional value of k is too 
small by 0”-04, as has been suggested, then the clock correction given by this 
group will show a sensible variation compared with the determinations made on 
normal groups at the same dates. It will at first decrease, by a total amount 
of 3 milliseconds, and will then increase by the same quantity. It seems possible 
to obtain evidence of this effect. 

It is not surprising that the determination of the constants requires long 
series of observations but apart from this long-term programme the data given 
in Table II have immediate applications: (1) to the determination of the 
irregularities of Universal Time; (2) to the study of the variation of latitude; 
and (3) to the establishment of a catalogue of stars referred to the provisional 
system of the astrolabe. 

The first of these applications has been studied by Dr Arbey, the two others. 
by Dr Guinot*. My two associates have refrained from making any selection of 
the observational data. They have rejected no observations and have retained 
incomplete groups as well as those observed by personnel under training. But 
they have weighted the results according to their standard deviations by means 
of a suitably decreasing scale. They have then computed the means, including 
from 10 to 19 groups but more generally 13. All these means have nearly the 
same weight and they are shown by the points on the following figures derived 
from astrolabe observations grouped as just explained. 

The diagram g(a) was derived by Dr Arbey and represents the irregularities. 
of the Universal Time TUO in Paris from 1956-5 to 1958-0. 

These observations have been referred to the Uniform Time given by Dr 
Essen’s atomic standard at the National Physical Laboratory which is at present 
the best tested of these standards. The time given by the crystal oscillators of 
Paris Observatory has been compared with Essen’s standard through the 
intermediary of the MSF time signals. Dr and Mme Stoyko have shown that 
the definitive terrestrial time from 1955-5 to 1957-7 requires a correction of 
the form A#*® to compensate for the effect of a progressive increase in the length 
of the dayt. The Universal Time TUO represented on Fig. g(a) has been 
corrected for this term. It is then directly comparable with the curve in Fig. 9(b), 
which represents the difference TUO-TUz2 computed by extrapolation by 
Dr Stoyko, in conformity with the instructions of Commission 31 of the Inter- 
national Astronomical Union for the rapid computation of the provisional time 


* A memoir by B. Guinot is in the press, to be published in Bulletin Astronomique. 
+ A. et N. Stoyko, C.R. Acad. Sci. Paris, 246, 235, 1958. 
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TU2*. Fig. 9(c) represents the difference of the ordinates of the curves g(a) 
and 9(b), that is the provisional time TUz. 

It is not possible at present to give an interpretation of this curve, as the 
final coordinates of the pole are not yet known, neither are the resuits obtained 
by other observatories which would enable the seasonal variations of the Earth’s 
rotation to be computed. It can, however, be remarked that the rate of 
provisional ‘TU2 is systematic and in particular an annual variation can be 
recognized. I have endeavoured to smooth out the TU2 irregularities and 
finally I have applied the following correction expressed in milliseconds : 


A=6sin 2n(t+0°15)+2°5 sin 4m(t +0°05) + 5(t—1957°5)” 
t being the time expressed in years. 
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The residuals, shown on Fig. 9(d), show no systematic trend and their mean 
quadratic value is of 3-0 milliseconds. I wish to stress that the diagram includes 
38 values for a period of 18 months, that is an average of 2 values per month. 

What is the significance of the quantity A? The first two terms will possibly 
be justified when the final value of the difference TUO-TUz2 is known, though 
it is possible that a portion of these terms should be attributed to errors in the 
group corrections and to the annual variation of the vertical. A slow variation 
would not affect in any way the closing error and therefore could not be detected. 


* N. Stoyko, Bull. Hor. du B.I.H., Série 4, 10, 203, 1956. 
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On this particular point the question remains open. As for the third term of 
the expression of A, that is the term in #?, it is to be added to the term identified 
by Dr and Mme Stoyko, the existence of which is confirmed by the present 
work. For the period covered by our observations, the gradual progressive rate 
of decrease of the terrestrial time is 0-52 ¢? milliseconds, ¢ being here expressed 
in months (instead of 0-484 ??). 

It is not without interest to have been able to establish these empirical 
corrections, which make it possible to pass from the TUO time of Paris to 
the uniform time of the Essen’s standard. When Ephemeris Time has been 
determined for the same period, the connection of the Essen standard to the 
conventional second of time will become easy. If we use the results of all the 
astrolabes now in operation throughout the world, including those to be 
delivered in the near future, it can be ascertained that this connection can be 
performed within 5 years with an accuracy of 10~°. 

But of all the problems dealt with in this lecture that of the time is indeed 
the most difficult. To obtain good time determination it is not enough to have 
available a good atomic standard and a good optical instrument. The chrono- 
graph, the time signal receivers, the propagation time of radio waves can each 
give rise to sources of errors or uncertainties. Great progress remains to be 
made in the transmission of time. In addition, with the degree of accuracy that 
we endeavour to obtain, the personal equations of the various observers are 
not negligible. 

I have endeavoured to show the personal equations by plotting separately 
the diagram of Fig. 9(d), for the three principal observers, Dr Arbey, Dr Guinot 
and Mile Débarbat. For certain periods, the same details appear on these 
curves as, for instance, the low values of April-May 1957, which very likely are 
real. But for other periods the curves do not agree so well. For instance, in 
January 1957, we find —7 milliseconds for the difference A-G while this 
difference is +4 milliseconds at the end of June 1957. In general, a systematic 
difference seems to exist between the observers, as for the whole of the 
observations of 1956-1957 the following values are obtained: 


A-G = — 2:9 milliseconds 
G-D = — 0-7 milliseconds. 


But these figures should not be interpreted as invariable personal equations. 
Indeed these differences are not stable; they vary fortuitously. Nowadays the 
name of “‘ noise ’’ is’given to functions characterized by such random variations. 
We could therefore speak of ‘‘ personal noise ”’ instead of ‘‘ personal equation ’’. 

It is very difficult to explain the persistency of such a marked difference as 
that of January 1957 between two well-trained observers. On one hand it is 
rather unlikely that the deflection of the vertical takes different values for odd 
and even days. On the other hand the observations are carried out in such a 
way that a systematic error in the appreciation of the coincidence of the two 
images affects only the radius of the circle defined by the position lines, but not 
the position of its centre. It must therefore be assumed that the observer does 
not time the coincidence in the same way on the east as on the west and that the 
error therein involved varies from week to week or from month to month. For 
the time being this effect remains without a satisfactory explanation. I have 
asked M. André Lallemand to study a photo-electric receiver designed 
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especially for the astrolabe in order to increase the accuracy of the observations 
and to obtain with confidence that of one millisecond. 

The determination of the latitude offers less complexity than that of the 
determination of time, as it is derived straight away from the observations 
themselves after the group junctions are completed. Fig. 10 was derived by 
Dr Guinot and represents the latitude of the Paris astrolabe since the middle of 
the year 1956. Our groups have been selected in such a way that they will be 
observable without any change until the year 2000. Homogeneous data will 
then be available for the determination of the secular wandering of the pole 
which is not the case at the moment. 
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Fic. 10. 


We have still to indicate what is the order of magnitude of the errors 
involved in both determinations. ‘The average standard deviation of a group- 
observation is 4:3 milliseconds in time and of 0-062 in latitude, including 
incomplete groups (the best determinations have given much smaller values, 
respectively 2-3 milliseconds and 0”-025). ‘These figures refer to the consistency 
of the values within one group. If we compare the observations of groups 
between themselves by referring them to the smoothed curve, then the standard 
deviations take larger values, particularly for time, viz., 7-6 milliseconds and 
0”:073. 

I shall now present the results of B. Guinot relating to the establishment of 
a catalogue of fundamental stars. In 1955 he presented to the General Assembly 
of the International Astronomical Union a note on the error Aw, of the FK3 
derived from the results of the prototype astrolabe which was the only one in 
service at that time. Fig. 11 shows the similarity of the results obtained since 
then with the OPL astrolabe, although the observed stars are in general different 
in the two cases. 

But B. Guinot has gone further. He has undertaken to determine with the 
astrolabe the coordinates of as many stars as possible, referred to the provisional 
system derived from the group corrections. He first computed the translation T 
which should be applied to each of the position lines to take into account : 

(1) The internal smoothing of each group. It has been checked that these 
corrections were not affected by magnitude equations or by colour equations. 

(2) The corrections of Table II. 
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He then derived from T the correction to the tabular coordinates of the star. 
If it has been observed at its two transits, east and west, and if 7), and Ty are 
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Fic. 11.— Aag (Astrolabe—FK3). 
Continuous line: experimental astrolabe (1954). 
Broken line: OPL astrolabe (1957). 


the corresponding values of 7, A is the azimuth of a transit reckoned from the 
south, positively towards the west, and 7 is the parallactic angle, then the 
corrections to apply to the coordinates of the catalogue are 


TF, 
Aa (Astrolabe — FK3)= ——2—-W 
a ee 3) 2|sin A|cos 
Aé (Astrolabe — FK3) = — Sgt sy. 
cos 7 


In certain special cases, advantage can be taken of a single transit, for instance 
when the star is close to its maximum elongation. In that case it can be assumed 
that cosy=o and the right ascension is then given by the following relation 


ie T + A6 cos 7 
Aw (Astrolabe —F K3) i ee 
sin A cos¢ 
The same formula gives Aé for stars close to the southern limit in declination as 
it then can be assumed that Aa sin A =o. 

B. Guinot compiled in that way a first catalogue of 115 stars of which 47 are 
known by their two coordinates, 56 by their right ascension only and 12 by their 
declination only. The Institut Géographique National has lent to Paris 
Observatory a second, OPL astrolabe, and special groups of stars, different from 
that of the routine programme, are now being observed in order to extend this 
catalogue to 300 fundamental stars, that is, more than half of those of the FK3 
and of its supplement contained in the Paris zone. Already, interesting remarks 
can be made. 

As for the declination, the astrolabe catalogue and the FK3 are in good 
general agreement. ‘The average of the differences is 


AS (Astrolabe — FK3) = — 0”-002 


and is not significant. This is only a check of the adjustment. Fig. 12 shows 
the behaviour of A6, for stars comprised in the zone + 20° to + 30°. Systematic 
differences are small but the statistics relate only to a limited number of stars. 
It is obvious that a constant error could affect the declinations of both catalogues. 
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In right ascension, the differences are larger and the distribution is more 
systematic. The chart in Fig. 13 gives a very expressive representation. Each 
star is there represented by a circle, whose diameter is proportional to the 
difference Ax (astrolabe — FK3), black when the difference is positive and open 
when it is negative. When the difference is smaller in absolute value than 
5 milliseconds the circle is replaced by a cross. We notice at a first glance a 
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region of high declination, comprised between 23 hours and 11 hours, where the 
differences are very considerable and positive. It is indeed because of this 
anomaly that the difference Ax, increases so rapidly beyond +50°. ‘This 
illustrates the disadvantage of dividing Ax into two terms. The FK3 system is 
homogeneous between 11 hours and 23 hours. In this interval it would be 
sufficient to correct the individual positions to obtain a good catalogue with 
very small systematic errors for the coordinates as well as for the proper motions. 
On the contrary, to apply the correction Aa, in this interval would be injurious. 
If the correction were applied it would be necessary to compensate for it by 
applying another correction Aw, of the opposite sign, this correction being 
a function both of « and of 6. 
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One may wonder why the astronomers have not hitherto realised the 
existence of these regional errors. ‘The reason is very likely because it had 
become the practice to analyse the errors in an arbitrary way by means of a 
formula instead of considering very simply only Ax and Aé. In my opinion, 
it would be advantageous to abandon this method which can in no way be 
justified and which is merely confusing. For the compilation of the FK4, it 
would be wiser to consider rather the individual errors, the local errors, and 
the systematic errors, the latter mostly arising from errors in the adopted 
positions of the pole and of the equinox. 

Fig. 14 represents the mean values of Ax cos 6 between +60° and +75°. 
There are two conspicuous discontinuities at about 11 hours and 23 hours. The 
representation of such residuals by trigonometrical series would of course be 
nonsensical. ‘The conclusion is that the conversion of the FK3 to the FK4 
cannot be derived from harmonic analysis. 















ms 
+30 r 
+20 Ax cosd 
{+ 60° to+75°) 
+10 
fe) nw 
Oh 12h 24h 
Fic. 14. 


1 will now make a last remark. In the case of meridian observations, the 
accuracy of which decreases from the equator as cos 6, the large discrepancies 
observed for high declinations can be allocated to random errors. But the 
accuracy of the measurements of right ascensions with the astrolabe depends 
only on sin A. It is the same at Paris for the zone from +20° to + 22° as for the 
zone +68°-5 to +73°- Hence if the Ax are much larger in the second zone 
than in the first, it is the catalogue that is responsible. The uncertainty of the 
group corrections could explain discrepancies of the order of 5 milliseconds but 
not of 20 milliseconds or more. ‘The anomaly of the FK3 is found again, but to 
a lesser extent, in the N30 catalogue of H. R. Morgan. The discontinuities are 
smaller and less abrupt. However, positive corrections should still be applied 
to the right ascensions between 23 hours and 11 hours beyond +50°. But the 
first astrolabe catalogue has as yet too few stars to allow for further detailed 
discussion. 

If I have taken the liberty of presenting so many provisional results, perhaps 
too many, it is not so much for their intrinsic value, for I do not ignore that they 
have to be confirmed*, as to emphasize the advantages of the method of equal 
altitude of stars. When this method is applied with the impersonal astrolabe, it 
provides for a permanent test, the zenith distance, which enables the quality of 
the results to be objectively evaluated. The several unknowns, namely time, 
latitude, corrections to the catalogue or to the fundamental constants, can be 
clearly separated from each other in a manner which cannot otherwise be 


* All the reductions will soon be revised using the high precision ephemerides which Dr Fricke, 
Director of the Rechen-Institut, has kindly undertaken to compute for the stars of the Paris zone. 
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attained in meridian astronomy. Instead of amending the catalogue hour by 
hour, it enables large areas to be dealt with, which are easily linked together. 
The results take an expressive form and their discussion leaves hardly any 
obscurities. It also gives to the computer an impression of confidence which I 
have observed with all my associates. I am convinced that this impression is 
not fallacious and that favourable prospects are really opened in the difficult 
field of fundamental astronomy, which is by no means a dead subject. 





A STUDY OF THE DETERMINATION OF 
SOLAR ATMOSPHERIC ABUNDANCES 
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Summary 


A study is made of the determination of solar atmospheric abundances 
from the observed intensities of faint Fraunhofer lines. ‘The method of 
weighting functions, which utilizes the results of the observed limb- 
darkening, is considered in detail. An alternative method is proposed, the 
method of the Planckian gradient, which makes no recourse to limb-darkening 
observations and which is entirely dependent upon the physical characteristics 
of the solar atmospheric model considered. 

Throughout the paper reference is made to the case of nitrogen (ionization 
potential = 14°54 volts) primarily to compare results obtained by the above 
two methods. Nine absorption lines are considered and the nitrogen 
abundance is determined from each of the nine lines. Appreciable dis- 
crepancies are found between the results obtained by the two methods: 
‘Theoretical considerations lead to the conclusion that in the case of nitrogen, 
and in general for all elements of high ionization potential, the method of 
weighting functions is quite unreliable and the Planckian gradient method 
should be used. 

Three atmospheric models are considered, and, for nitrogen, the mean 
results obtained by the latter method are (on the basis log A},;= 12-00) 
Claas model : log dy =7°71, Vitense model : log Ay =7°81, Swihart model : 
log Ay =7°93 (cf. Hunaerts 9:02, Unséld 8-61). 





Introduction.—A method for the derivation of solar atmospheric abundances 
from observed faint Fraunhofer line intensities was proposed by Unsold (18) 
in 1932 and later modified and improved by Minnaert (10). This method, 
that of weighting functions, has been used extensively for abundance determina- 
tions (e.g. Claas (5)). Here the information provided by the observed limb- 
darkening of the Sun is used and the method may be subject to certain 
limitations. It may be, for example, that the limb-darkening results, which are 
reliable only for the outer layers of the atmosphere, do not give adequate 
information concerning phenomena occurring in the deeper layers. The theoretical 
basis of the method of weighting functions is again considered in detail in 
Section 4. 

An alternative method is proposed, here called the Planckian gradient method, 
which is determined entirely by the physical characteristics of the particular 
solar atmospheric model considered. ‘This procedure makes no reference to 
the results of limb-darkening and is therefore free, in some degree, from the 
above limitations (Section 5). 

The two methods are applied to the case of nitrogen, which, due to its high 
ionization potential, will provide appreciable absorption in the lower layers of 
the atmosphere, and the numerical results may then be compared. A theoretical 
comparison is also made in order to obtain a measure of the reliability of the two 
methods. 
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1. The solar atmosphere.—The energy which is radiated from the surface 
of the solar atmosphere, together with its angular and spatial distribution, gives, 
at least in principle, complete information concerning the radiation from each 
elementary layer in the photosphere. Most of the radiation from the upper, 
cooler, layers will reach the surface of the photosphere; the radiation from the 
lower layers will however be modified by absorption in the intervening strata. 

We consider a plane parallel atmosphere in local thermodynamic equilibrium, 
i.e. to each layer in the atmosphere we assign a local temperature 7 such 
that, for this temperature, Kirchoff’s law is valid, i.e. j,/k,=47B,(T) at 
any frequency v where k, andj, are the mass absorption and emission coefficients 
respectively and B,(7) is the Planckian function at temperature 7. The 
emergent intensity of radiation in the continuum at an angle @ to the outward 
normal is given by 


Ic(0, 0) = [“ B(T)e-r-s° *sec dr, (1.1) 
J0 
where, in material of density p, 


— [ Apa (1.2) 
0 


expresses the transition from the linear depth x to the optical depth 7, in the 
continuum at frequency v. Here B,(T) is the source function for the radiation. 

In the case of the solar atmosphere the continuous absorption coefficient k, 
is almost entirely determined by the photodissociation of the negative hydrogen 
ion and by photo-ionization from discrete levels of the neutral hydrogen atom 
(bound-free transitions). In the present investigation for faint absorption 
lines it is assumed that the condition of local thermodynamic equilibrium is 
satisfied, i.e. the frequent interchange of energy between the absorption electrons 
and the free electrons ensures that the population of the energy states will be 


ces in accordance with Boltzmnann’s law, and scattering is unimportant. ‘The 
8) formation of the Fraunhofer lines is then due to selective absorption followed 
od, by continuous absorption. If /, is the mass absorption coefficient for line 
1a- absorption in frequency v then the emergent intensity in the line, at an angle @ 
‘b- to the outward normal, is given by 
ain ls er oe a 

I ,(0, 0) = B ,(T )e-*v 8e¢ ° sec 6 dt,, (1.3) 
are 4 
ate where 
cal re 
in t= | (k,+1,)pde (1.4) 

“0 


is the optical depth in the line. 


rd, 2. Solar atmospheric models.—Using equation (1.1) it is possible in principle ' 
lar to obtain the temperature distribution with optical depth for any wave-length 
to from the observed limb-darkening of the Sun. Chalonge and Kourganoff (4) 
he and Barbier (2) have made most careful studies of the limb-darkening observations 
and many models are based on these results. 
gh For example, we may represent the Planckian function in the form 
. B,(r,)=a,+6,7,+¢,7,2 (2.1) 
cal ws ee ee y 
-_ and derive the values of the constants a,, 6, and c, in accordance with the method 


outlined in Section 4.1. For a given wave-length the Planckian function B,(7,) 
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is a function of 7' only and hence equation (2.i), a quadratic in 7,, gives the 
solution 7,=7,(T). 

Since we may obtain the relationship 7,=7,(7') for any wave-length, we are 
able to compare the optical depths in the several wave-lengths corresponding to 
a single value of 7, i.e. toa single layer in the atmosphere. Since 7, ,/7, 2=h,,4/k,.» 
we are also able to determine the ratio of absorption coefficients in several wave- 
lengths. As we are able to correlate the optical depths in the various colours 
we may select one wave-length as typical (A= 5000 A) and express the temperature 
variation in the photosphere in terms of this optical depth 7) (see Section 4.2). 

Since it will be of interest, in the abundance determinations which follow, 
to examine the sensitivity of the results to the model atmosphere chosen, three 
models will be used with slightly different characteristics (especially surface 
temperature). 

(a) Model adopted by Claas.—The model adopted by Claas (5), which we 
will refer to hereafter as the ‘‘ Claas model ’’, is based on the limb-darkening 
calculations of Barbier (2). With regard to the structure of the atmosphere 
only two assumptions are made, i.e. the layers are in radiative and local thermo- 
dynamic equilibrium. The calculations of the gas pressure P, and the electron 
pressure P, were carried out in accordance with the method developed by 
Barbier (2) except that a fourth approximation based upon Barbier’s third 
approximation for P, and P, was made. 

‘The model is defined by Claas from the surface down to optical depth 7, = 2-2 ; 
this was found to be sufficient for the application of the method of the weighting 
function (Section 4), but quite insufficient for the method of the Planckian 
gradient (Section 5). Below this depth the relationship 7, = 7,(7') derived from 
limb-darkening observations is rather uncertain. Lyttkens (g) has extended the 
model to greater optical depth using the relationship 


T= 37 Alr + q(7)], G4) 
where 7), is the effective temperature of the Sun (7.,=5713 K) and 7 is the 
Rosseland mean optical depth. ‘The function q(7) varies between 1/4/3 at 
7=0-'0 and o-710 at infinitely large optical depth, for the present model the 
value 0-710 is adopted for the whole range below ty=2:2. Lyttkens further 
calculated the P,=P,(7) relationship utilizing the tables of Stromgren (15). 
The electron pressure can now be found directly from the same tables (‘Table 8 
in Strémgren’s publication). 
By conversion to the standard optical depth 7) in wave-length A= 5000 A 
by the relationship 
Ty [ear (2.3) 
Jok 
utilizing tables of k,/k=k,/k(0,P,) kindly supplied by Dr B. Strémgren 
(unpublished), and by numerical interpolation in the results of Lyttkens, the 
Claas model was extended in regular numerical interval in 7, to the depth 7) = 4:2. 
(b) Vitense model (20).—The method adopted for the determination of the 
T= T9(7’) relationship is essentially the same as the method used in the foregoing 
model. Here however the limb-darkening results of Chalonge and Kourganoff 
were preferred to those of Barbier and corrections were made to conform to 
the more recent values of the radiation constants in the Planck function. Vitense 
(20) however took into account that the measured values of /%(0, @) are available 
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the only for 0-3 <cos@< 1-0, and therefore the derived function B,(7,) is reliable 
only in the corresponding region of 7,, i.e. approximately 7,>0-2. The 
are temperature stratification of the uppermost layers, o< 7, < 0-2, was found from 
to the limb-darkening variations of the residual intensities of the strong Balmer 
lines of hydrogen and the strong lines of iron. ‘This resulted in a much lower 


hg value of the surface temperature, viz. 7,= 3800+ 100 K, compared with that 
urs used by Claas (7y=4900°K). Using this value of the surface temperature, 
ire the intermediate temperatures, i.e. 0 < 7, < 0°2, were found by interpolation. 

2). In the model used in this investigation (Model (ii) in Vitense’s publication) 
w, a slightly steeper temperature gradient, found by D. Labs, was adopted, resulting 
ree in a slightly higher surface temperature (7,= 4006 “K). It was found however 
ace that the uncertainties in the temperature and pressure in the superficial layers 

do not influence the distribution in the region 7, > 0-2. 
we Again the model derived is defined only from the surface to optical depth 


ng t,=2°0. The results of Lyttkens were again used to extend the model to much 
lower layers, the fitting being made at log P,, = 5-10 for #=0-706 (a small error 


ere 

10- may have been introduced here). 

on (c) Swthart model.—Swihart (16) produced, initially, a standard model of 
by the solar atmosphere in which the observed integrated flux of radiation is 
rd reproduced by the assumed += 7(7’) relationship. 


The flux FY at the surface is obtained by integrating the surface intensity 
2: over all solid angles, i.e. 











nv n/2 : 
ma F= [ Fe, 0) cos 6 dw =2n | I’(0, 0) cos @sin 0 d0=nF. (2.4) 
” 0 
Introducing equation (1.1) we have for the monochromatic flux 
1¢ 
F,=2 | B(T)dr, | ery see sin Odd, (2.5) 
Z) . 0 0 
_ If we now introduce the variable 
ot secO=y, ie. sin@d0=y~* dy, 
he and the exponential integral function 
Ler "OD @-LY 
5), E, (x)= i; se ys (2.6) 
> 8 equation (2.5) becomes 
" F,=2]| B,(T)E,(r,)dz,. (2.7) 
“0 

More generally for the flux at any depth 7,, where there is an inward flowing 
3) contribution from the layers above 7,, we have 
en Fi(7,)=2| B,(T)E.(t,—7,)dt,-—2 | "BA(T)E(7,—t,) dt,, (2.8) 
he J Ty “0 
2. and the monochromatic flux may then be determined once the 7 = (7) relationship 
he is specified. ‘lhe integrated flux is then found from the expression 
ng © 
ff | M(r,)= | F,(r,)dy, (2.9) 
to : ‘ e wile 
pa and if the assumed 7=7(7’') relationship is incorrect then the derived flux will 
le also be incorrect. 
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If the r=7(7)) relationship is written in the form 
T4=3TA[a+br+cE,(7)+dE;(7)], (2.10) 


(cf. equation 2.2) the coefficients a, 6, c, and d may be found from the condition 
that the observed flux is reproduced by equation (2.8). 

Swihart first constructed a standard model, a plane parallel atmosphere in 
radiative and mechanical equilibrium, by the above procedure, taking into 
account the uncertain blanketing effect by increasing the flux at all depths by 
approximately 10 per cent. He later (17) modified this model by introducing 
an assumed, plausible, convective transport in the lower layers (~1 per cent 
at T= 1-5 increasing to 40 per cent at r=3°0). This had the effect of reproducing 
the desired flux rather more accurately than the standard model but increased 


TABLE I 


Solar atmospheric models 


Claas Model Vitense Model Swihart Model 

To 6 Log P, Log P, 6 Log P, Log P, T 6 Log P, Log P, 
O°OI 1°022 3°92 0°09 1'072 3°95 oO'10 O'o!1 1'098 3°93 —0°03 
0°02 I°OI5 4°07 0°24 1'°046 4°15 0°32 0°02 «1°088 6409 o'12 
0°03 «6T'OrI: «64°16 «0°31 1°035 4°23 O°'4! 0°03 «:1°078)«=—-4"19 0°22 
0°04 1°007 4°23 0°38 1029 4°30 0°48 0°04 1°070 4°25 0°28 
0°05 1°003 4°28 0°42 1°022 4°36 0°53 0705 1°062 4°31 0°34 
0°06 =0°999 4°33 0°47 1°O15 0°57 0°06 §=1'054 «44°35 «= 0°39 
0°07 I°OIO 64°43 0°60 0°07 1°048 4°38 0°43 
0°08 0990 4°39 0°54 1°004 0°64 0°08 1'°040 4°42 0°47 
0°09 0°999 0°68 0°09 1'034 4°44 0°50 
0710 o'980 4°45 0°61 0°995 4°53 O°'7I O10 «=1°028 84°47) 0°53 
0:20 0°942 4°61 0°81 0°949 4°68 0o-go 0:20 0980 4°63 0°75 
0°30 O'g10 «64°71 «20°98 o'915 4°78 1°05 0°30 0°946 4°73 0°89 
040 0°884 4°78 #I'1l 0°887 44°86 1°17 0°40 0916 4°80 1°02 
0°50 0860 4°82 1°24 0°864 4°91 I°3I 0°50 600°894) «4°85 1°12 
060 0840 4°86 1°37 0°845 4°94 1°43 060 0°874 4°89 1°22 
0°70 0829 64°97. 1°55 0°70 0°854 4°92 1°32 
080 o808 4:91 1°58 0813 4°99—«11:°65 080 8 0°838 = 4°95 1°42 
0"90 °0°799 +~5°OI 1°70 0790 0°824 4°97 I°51 
1700 60781) «=644°94-—S «1°77 0787 «5°02 «1°78 1°00 6—0'810 499 1°60 
1°20 0°757 4°96 1°94 0°765 5°05 1°93 1°20 0°784 5°02 1°77 
1°40 0°736 4°97 2°09 0°746 5:07 2:06 1°40 0'760 5:04 1°94 
1°60 0720 4°99 2°22 0°729 2°17 1°60 0°740 5°05 2°08 
1°80 0°706 5:00 4°32 O'715 5°10 2°27 1°80 0°724 5:06 2°20 
2°00 0694 5°01 2°41 *o-702 5°10 2°40 2°00 O'710 5°07 2°30 
2°20 *0°684 5°01 2°49 0°693 «S'II) = 2°47 2°20 0697 5°08 2°40 
2°40 0674 5:02 2°56 0°683 5°12) 2°55 2°40 0'687 5:08 2°48 
260 0665 5:03 2°63 0674 «5°12 2°61 260 0679 5°09 2°54 
2°80 0658 5:03 2°68 0°665 5°13 2°68 280 0672 5:09 2°59 
3°00 o°651 5°04 2°74 0°658 5°13 2°74 3°00 0°666 5:10 2°64 
3°20 0°644 5°04 2°80 o'651 5°14 2°86 3°20 0662 5:10 2°67 
3°40 0°637 5°04 2°85 0°644 5°14 2°85 3°40 0659 511 2°70 
3°60 0631 5°05 2°90 0°637 5°14 2°90 3°60 0656 511 2°72 
3°80 0°624 5°05 2°95 0°631 «5°15 2°95 3°80 0°653 5°11 2°74 
4°00 oO619 5°05 2°99 0°624 5°15 2°99 4°00 0650 5°12 2°77 
4°20 0613 5°05 3°03 5°00 0°642 5°13 2°83 

6:00 ©0°632 «5°14 2°92 
* Extension of model to greater depth by Lyttkens. 7°00 0°620) 5°15 3°00 
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the divergence of the predicted limb-darkening results when compared with the 
observed limb-darkening. Swihart concluded that the actual convective transport 


10) could not be much larger than assumed in this model. 

ion ‘The model was recomputed with regular numerical interval in 7 using the 
values of a, b, c and d obtained by Swihart. ‘Tables of the first exponential integral 

in function are available and functions of higher degree were calculated from the 

nto relationship 

by nE,,,.,(7)=e— — TE, (7). 

ing The complete models are given in Table I. 

ent 3. The nitrogen absorption lines.—The electron configuration of the ground 

ing state of the nitrogen atom is 1s?2s?2p%. The electrons in addition to the closed 

sed shells, i.e. the equivalent electrons 2p*, give rise to possible total spin quantum 
numbers S=}, 3 and therefore to multiplicities, (2S+1), of doublets and 
quartets. By the Pauli exclusion principle the three equivalent electrons may 
give rise to the terms ?P, 2D, and 4S only, and, as the term with the greatest 
multiplicity designates the lowest energy state, the ground state of the nitrogen 
atom is designated 1572s?2p*4S},.._ Here the subscript refers to the value of the 

P, inner quantum number /, and the superscript ‘‘o”’ refers to the fact that the term 

3 is odd, i.e. }/; is odd where the summation is extended over all the electrons 

> of the atom. Excited states are obtained when one of the 2p electrons moves to 

8 a higher orbit; in addition a 4P term is indicated where one electron moves from 


the 2s to the 2p orbit. ‘The electron configurations, term designations, and energy 
levels of all the states of nitrogen are available in the publication Atomic Energy 
Levels (13) and all the relevant information is taken from this source. 

Transitions between two atomic levels, producing a spectral line, are subject 

to the selection rules (for dipole radiation): total angular momentum AJ =o, +1 
with the restriction J =o->J=o forbidden; spin angular momentum AS=o, 
i.e. only terms of equal multiplicity will combine; orbital angular momentum 
AL=o, +1. Inaddition for the electron making the quantum transition A/= +1, 
this may be interpreted as A(>J;)= +1, i.e. even terms will combine with odd 
terms only, and odd terms with even terms only. Table II gives the possible 
multiplet transitions for nitrogen, together with their component lines, within 
the wave-length range of the present investigation. 

The identification of the lines is based upon measurements made by 
Professor Edlén and supplied to the author by Dr B. Stromgren. The strongest 
Fraunhofer line, A=8680-24A, is seen to have a considerable wave-length 
discrepancy when compared with the solar spectrum (—0-14A), thus placing 
the identification in some doubt. It will be seen later, however, that the omission 
of this line has no effect upon the abundance determinations. The excitation 
potentials refer to the lower level of the transition and have been corrected to ° 
correspond to an ionization potential of hydrogen of 13°595 volts (i.e. by 
multiplying the level by the factor 0-00012395). 

The oscillator strengths f of the transitions in question can be shown, by the 
well-known analysis of Einstein, to be related to the strength of the spectral 
line S, by the equation 

2 
a ee (3-1) 
3 he gia 


where m, c, e, and h have their usual significance, and g; is the statistical weight 
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3s'P-3p!S° 
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TABLE II 


Nitrogen 1 


Vol. 118 


Ionization potential 14°54 volts 


Desig. J Level cm=! Config. Desig. J Level cm! 

2p>4S® 13 ° 28°2p?("P)3p 3p°S” 3935823 

2p>*D® 23 19223 2s*2p°(°P)3p 3p*D® 3 94772°2 
I} 19231 13 94794°8 

2p>*P° 1s, 2$ 948321 
i} om 33 948831 

3s *P 5 —-83285°5 2s°2p*(*P)3p 3p°P® 3 95476°5 
13 83319°3 1 954949 
23 83366-0 2  95533°2 

3s *P 5 861314 2s°2p°(*P)3p 3p 4S? 13 = 96751°7 
1} 86223°2 2s°2p°(®P)3p 3p7D® 145 96788-2 

2p**?P 24 88109°5 24 =. 9 6864-2 
1s 881534 2s*2p*(*P)3p 3p*P® 3 97770°1 
j 88173°0 1} _g7805°8 

N 11(°P°) Limit... 117345 
Absorption lines 
Laboratory Sun 

AJ Wave-length Excit. Pot. Sun — Lab. Identification 
A(A) (Volts) AA 

t-1} 7423°63 10°32 —o'12 Si Masked 

14-1} 7442°28 10°33 —0'05 Nitrogen 

24-1} 7468:29 10°33 —0'02 Nitrogen 

13-2} 8184-80 10°33 —o-oz2 Atmos. N 

dr} 8187-95 10°32 —o-'10 Atmos. N Masked 

I 3 8200°31 10°32 —0°32 Atmos. N 

14-1} 8210764 10°33 0°32 Atmos. N 

24-2} 8216:28 10°33 002 = Nitrogen 

1i- } 8223°07 10°33 —o'19 Atmos. N 

23-15 | 824234 10°33 o-o2 ~=©Atmos. N 

ty} 8567°74 10°67 0°04 Fe 

j- 4 8594°01 10°67 ; Absent 

1}-1} 8629-24 10°68 —o0:08 —— Nitrogen 

13-3 8655-88 10°68 0°17 Blended. Wing Ca 

24-33 8680:24 10°33 —o'14 Nitrogen 

13-2} 8683°38 10°33 0:00 ~=—— Nitrogen 

t-y} 8686-13 10°32 0°24 Nitrogen Masked 

j- 3 8703°24 10°32 —0'09 Nitrogen 

14-14 8711-69 10°33 —o'02 #8 Nitrogen 

23-2} 8718-82 10°33 —0:06 = Nitrogen 

1}- } 8728-88 10°33 —0:'28 Fe 

23-1} 8747°35 10°33 og = Fe 
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of the lower level concerned in the transition. The statistical weight is equal to 
the degree of degeneracy, (2/ +1) of the level in question. Introducing the 
numerical constants in equation (3.1), expressing the wave-length A in Angstrém 
units and S, in atomic units (a,”e”, where dy is the radius of the first Bohr orbit) 


we have 
3 _ 3°04 x 1072S, 
cem7} f giA ; (3-2) 
32°3 The line strengths S, may be expressed in the following form: 
12,°2 
4°8 S, = S(M)S(L)o?, (3-3) 
32°1 
33° where S(M) is a factor which depends upon the particular multiplet concerned 

P P P P 

16°5 in the transition and S(L) is a factor which depends upon the particular line of 
edi the multiplet, o? is a function of the radial wave functions of the initial and final 
a“ “ states of the active electron and upon the greater of the two azimuthal quantum 
ade numbers. S(M) may be found from tables in two papers by Goldberg, the 
38-2 a 7 * PP y § 
14°2 first (6) giving the relative strength of the multiplet and the second (7) a factor 
orl for the conversion of relative strengths to absolute strengths. S(L) is obtained 


5 “8 from Russell (14) where the line strength of a particular transition is expressed 
as a fraction of the sum of line strengths comprising the multiplet. An 
approximate evaluation of o” has been made by Bates and Damgaard (3) for a 
large number of transitions in light elements and it is found that the value of o? 
is the same for all multiplets of the 3s—3p transitions which are of interest in the 
present investigation. The oscillator strengths may therefore be calculated 
directly from equation (3.2). 

The equivalent width W, of a spectral line is by definition the width of a 
on perfectly black line, i.e. of rectangular profile, which, in the same region of the 
spectrum as the line itself, would absorb the same amount of energy from the 
continuum as does the spectral line, i.e. 





+O AT, 
W={ etd (3-4) 
sked eats 
= The nine nitrogen lines identified in Table II are measurable in the Utrecht 
Photometric Atlas (11). The nine lines, however, are faint and difficult to 
measure with accuracy (Table III). 
TasB_e III 
Atomic parameters and equivalent widths 
g Ca 
Transition 
A (A) 38-3p S(L) S(M) a o S; gi Z W, (mA) 
ked 
7442°28 — *P3jo-*S% 5/0 0333 2a 0-2 «9°76 = 1300 4 0°133 3 
7468:29 = *P5/2-*S° 3/0 0°500 2a 2 «= 9'76-Ss 19°52 6 0132 3°5 
8216-28 = *#P 5 jo—*P 5/0 0°350 6a 2 9°76 40°99 6 0°253 6 
8629°24 = *P3jo-*P®52 «= 0556 3a 2976 3256 4 0287 4 
8680-24 *P5/2-*D° 7/2 07400 10a 2 9°76 78:08 6 07456 17 
8683°38 = *P3/2-*D® 5/2 0'210 10a 2 9°76 40°99 4 0°359 7 
8703°24 = *Py/2—*D°/2 0°083 10a 2 9°76 16°20 2 0°284 5 
8711°69 = *P3/2-*D 3/2 0107. I0a@ 2 9°76 20°89 4 o182 4 
8718-82 *P5).—-*D° 5/2 07090 10a 2 9°76 17°57 6 o'102 4 
30 
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4. The method of weighting functions.—We wish, initially, to consider the 
well-known method of weighting functions for the determination of solar 
atmospheric abundances from weak Fraunhofer lines. The relationship between 
the intensity of an absorption line and the structure of the solar atmosphere is 
simpler in the case of weak absorption lines than in the case of strong lines. 
For any given solar model we may assign a weight g,(7) to each elementary layer 
of the atmosphere which expresses the contribution of that layer to the emitted 
radiation, in any wave-length, either in the continuum or in the absorption line. 
For a faint line in pure absorption we have therefore 

AI, _ I&(o,@)—I,(0,8) _ 
T(0,0) ~*‘I°(0, 8) =| gAepee- (4-1) 

In his development of the theory of weighting functions Minnaert (10) 
obtained the form of the weighting function essentially in the following way. 

The intensity in the continuum, at the centre of the solar disk, is given by 
(cf. equation (1.1)) 








1,"(0,0)= | B,(r)e*r dz, (4.2) 
J0 
and the intensity in the line is 
1,(0,0)= | B,(r)e-tedt,. (4-3) 
0 
The optical depths 7, and t, are related by 
dt, _k,+l, ; 
alias eae where /, <k,, (4.4) 
i.e. 
t=r,+ ["» dr,. (4.5) 
v y Je k, v 


We may write this in the form 
t=7,+7, where 7,= \’ > dr,. (4.6) 
0 


v 


Equation (4.3) then becomes 
I ,(0, 0) = [- B,(r)e vt? (: + ) dz, 
=| B(r)e77 (aa) (14+ ? Jar, 
s [ B,(r)e-*»dr,— | * B}e~"(71- :) dr,, 
0 0 k, 


to the first order. The intensity dip in the continuum is nennene 
AI, =I*(0, 0) —1,(0,0)= i B,(r)e-"»7,d7, — in B,(r)e-*» a “+ dr, (4.7) 
and transforming this equation by use of the relationship 


| y(7)2(7) dr = {- dz | ydr, 


we have 


= [dr _B A(t)e-7» dr, — J. y B (rev dr, 


tr isl ~B,(r)e-" |ar, 4.8) 
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he an equation which was also obtained by Unséld (18). Expanding the first term 
lar of equation (4.8) we have 
en - 00 ~ 00 
is | Bucjeedz,=— | B,(r)d(e~) 
es si ™ 
; * dB °. oe. 
aes =— | By (rje7»]  — |= er] + | aie "vdr, 
eS Ty v Ty Ty v 
dB, . aB 
ne. — “ > +... Eris 
[ B+ dr, + dr3 + |e 
) Equation (4.8) then becomes 
” ‘= dB, . dB 
Al, = “ “+... |e- dz, . 
ay. Similarly the intensity in the continuum is given by 
by es) © 
Is(0,0)= | B,(r)erdr,= — | B,(x)d(e~») 
0 0 
}-2) dB, dB, 
= [8.4 +5 te] me (4.10) 
i.e. . 
t-3 +00 
- | (B,’+B,"+B,"+...Je-"vdr, 
v a 0 
Kec)” (8,48, +8, +... ese) 
+4) Unsold indicated how this integral may be evaluated assuming the Planck 
function to be an approximately linear function of optical depth. Minnaert 
((10), note also (12)) improved upon this and considered the form of the Planck 
4-5) function as a power series, i.e. 
B,{(7,) =a,+6,7,+¢,7,3+4,7,>+e,7,4, (4.12) 
1.6) giving 
; B,(7,) =b,+2¢,7,+34,7,? + 4e,7,5 
B,"(t,) =2c,+6d,7, + 12e,7,? 
B,"(z,) = 6d, + 24e,7, 
B,!"(7,) ie 24e,. 
Introducing the following constants 
A,=6b,+2c,+6d,+24e, } 
B,=2c,+6d,+24e, | 
C,=3d,+ 12é, } (4.13) 
D,=4e, 
4-7) and N,=a,+b,+2c,+6d,+24e, | 
equation (4.11) simplifies to the form 
[- [A,+B,7,+C,7,2+...Je-» dr, - 
ei = 1 = [os (r,)dz. (4.14) 
E(0, 0) N, — 
It is clear therefore that a knowledge of the coefficients a,, 6,, c,, d, etc. of the 
expansion representing the Planckian function is sufficient completely to 
8) evaluate the weighting function at any depth 7,. This information may be 
4 obtained from the observational investigation of the limb-darkening of the Sun. 
30% 
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4.1. Limb-darkening.—The intensity of radiation emitted by the solar 
atmosphere in the continuum at frequency », in a direction at an angle 6 with the 
outward surface normal is given by 

I5(0,0)= | B,(z)e-7» 8¢ 8 sec 0 dr,. (4.15) 
This has a solution of the form 
I°(0,0)=A,' +B,’ cos8+C 


i.e. direct integration of equation (4.15) yields 


’cos?0+..., (4.16) 


v 





I°(0, 8) = \- B(r)e ru OT, where p= cos 6 
p 


=~ |" By(r)der) 


) 
2 dB, , @B, : 
= [8,40 dr, +p dr,” tee] (4-17) 


which is of the same form as equation (4.16). Introducing the Planck function 
in the form given in equation (4.12), we have 


I¢(o, 0)=a,+ 6, cos@+2c,cos*6+.... (4.18) 


Observations of J¢(0,@) have been made by a large number of workers. The 
results obtained by Abbot and Raudenbusch have been most carefully analysed 
by Chalonge and Kourganoff (4), who, expressing /¢(0,4) (or more exactly 
¢, =I1%(0, 4)/I<(0,0)) in the form of equation (4.16), obtained the values of 
the coefficients A,’, B,’, and C,’ for a large range of wave-lengths. 

From equations (4.16) and (4.18) we see that A,’=a,, B,’=6,, and C,’ =2c, 
etc., and therefore the constants of equation (4.14) may be evaluated and the 
value of the weighting function obtained at any depth 7, in the atmosphere. 

An advantage of using the observational results of limb-darkening is that th's 
minimizes the number of uncertain assumptions which are required to be made, 
e.g. no uncertain factors depending upon turbulence or the blanketing effect 
need be considered. The obvious disadvantage is that, from the nature of the 
observations made, the values of the coefficients will be accurate for only the 
outer layers of the atmosphere, and for large optical depths the method may be 
unreliable. 

4.2. Transition to the standard optical depth.—In order to be able to compare 
calculations made in different regions of the spectrum it is advantageous to use 
one variable for the optical depth in all spectral regions. ‘The optical depth 7» 
in wave-length A= 5000 A is chosen as the standard variable; this also facilitates 
numerical calculations since many of the model atmospheres are defined with 
T= 7 (9) expressed in regular numerical interval. (In the case of the Swihart 
model however the temperature distribution is referred to the Rosseland mean 
optical depth t= 7(@) and, for this model only, all spectral regions are referred 
to this optical depth.) 

In some particular wave-length A, corresponding to frequency v, equation 


(4.14) gives 








es) - | #r)dr= [so ear (4.19) 
where " , . 
&,(7)=8,[7(70)]- (4.20) 
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Equation (4.20) expresses the transition of the weighting function for wave-length 
A=x, referred to optical depth 7, (i.e. in wave-length x), to the weighting function 
for A=x referred to optical depth 7, (i.e. in the standard wave-length A = 5000 A). 
Corresponding values of 7) and 7, are easily calculated, since 


te k 
Tj;)= | ZU (4.21) 
0 v 


and tables of the absorption coefficients k,=k,(6,P,) are available (Strémgren 
(15)). ‘The weighting function may thus be expressed with 7, as depth variable 
and this function is designated g,(t9). If, following Claas (5), we write 





& To) (7) =h (Te) (4.22) 
Ry 
then equation (4.19) may be written in the form 
Al, = l, k, (2 
T}(0, 0) e \, 8.70) k,” ke al J, feral dro. (4-23) 


4.3. The line absorption coefficient l,.—If x, is the selective absorption coefficient 
per absorbing atom, and N, the number of absorbing atoms per cm’, then N;/N my 
is the number of absorbing atoms per gram. That is, due to the overwhelming 
abundance of hydrogen in the solar atmosphere, it is assumed that the density 
of the atmospheric gas is N,m,gcm~*. The selective absorption coefficient per 
gram is therefore 

— 

= (4-24) 

It has been established (e.g. Aller (1)) that the total absorption by an atom 

radiating as a classical dipole is given by 
+ © mez 2 

[= M= Zo. (4.25) 

Quantum theory however ‘requires that the absorption of N; atoms is the same 
as the absorption of N,f classical oscillators, i.e. equation (4.25) becomes 








mre")? 
i | x, d= Sf, (4.26) 
and we may therefore write 
+ © +o WN. qre2r2 ra N; 
l,dvA= ix d\= — a ' 
i o- 3 o Numy " me* * my Ny (4-27) 


4.4. Introduction of Boltzmann’s law.—The absorption line intensity depends 
not only on the particular atomic transition probability but also upon the number 
of atoms in the initial absorbing state, N;.. Boltzmann’s equation expresses the 
number of atoms in absorbing state 7, N;, in terms of the number of atoms of the 
same element in all states n, i.e. the total number of atoms N: 

—EJKT —EKT 
7 « fees « (4.28) 
dg,e~%n u 
U 





where g,, is the statistical weight of the level n, u= > g,e~"n/*” is the partition 


function of the element in question, and £; is excitation energy of the level 7. 
The population of a level may also be influenced by the perturbations produced 
by neighbouring particles; this however will be of importance only for states 
of large principal quantum number and these levels produce spectral lines in 
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the far infra-red, i.e. far beyond the wave-length range considered here. There 
is therefore no perturbation correction required for the population of the initial 
level in the present investigation. Claas (5) has obtained values of the partition 
function of nitrogen as a function of temperature and these values are used in 
the present discussion. The introduction of the parameter 6 = 5040/T enables us 
to express the excitation potential directly in volts, for if E, is the excitation us 
energy (in electron volts) then 
log,)¢~#/*7 = — V8, 

(V in volts) and equation (4.28) becomes 

N; _ g;x 10-8 

oe ill (4.29) 

N u 


4.5. Determination of relative abundance.—The equivalent width of a spectral 
line is expressed by 


+o AT, 
a 


Introducing equation (4.23) this becomes 


W,= | is | ” h,(t9)l, d75 aA, (4.30) 
—aov0 


h,(79) is of course practically constant over the width of the line and therefore 
by equation (4.27) we have 


da. 





| 





22 f oe) N. 
W,=~ a h, (79) <= dro. , 

A mc2 my Jo (70) Nu To (4 31) 
If A is the relative abundance of nitrogen with respect to hydrogen, i.e. A = N/N,, 
(at temperatures prevailing in the solar atmosphere the degree of ionization ci 
nitrogen, I.P.=14:54 volts, is very small indeed and we may consider all the 
nitrogen atoms to be neutral), we have then 


N; N; 
7“? (4.32) 
Utilizing equations (4.29) and (4.32), equation (4.31) becomes 
me*r? f A poate so ¥ 


mc? my ° u 





w- dry. (4-33) 
The numerical value of the constant ze?/mc?m, is 5:28 x 1011; if however the 
equivalent width W, is expressed in mA and the wave-length A in Angstrém 
units we have finally 


- 00 —Ve 
W,=5:28 x 10° . Af g; A | Sfre)x zo" dt. (4.34) 
0 


This equation expresses the equivalent width of the absorption line in terms of 
the weighting function, various other parameters, i.e. A, f, g;, V which are character- 
istic of the particular atomic transition, and the relative abundance of nitrogen. 
The weighting function may be calculated from the observed limb-darkening 
(Section 4.1), the parameters of the atomic transitions and the equivalent widths 
are known (Section 3), and hence the relative abundance may be determined. 
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4.6. Numerical calculations.—In the wave-length range under consideration, 
i.e. 7442-8719 A, it was found sufficiently accurate for the determination of the 
weighting function g,(7,) to adopt the wave-length A= 8206 A as representative 
of the range. 

(a) g,(7). For A=8206A the investigation of Chalonge and Kourganoff (4) 
gives the results 


A,'=a,=0°'501, B,' =b,=0°666, C,’ =2¢, = —0°172. 


v 


The weighting function g,(7,)=(A,+B,7,)e~*»/N, may then be evaluated at 
any value of the optical depth 7, since 


A,=b,+2c,=0°494, B,=2¢,= —0°172, N,=4,+6,+2c,=0°995. 


Corresponding values of 7, (A=8206A) and t) (A=5000A) are determined 
from equation (4.21) and hence the weighting function is determined for a large 
range of values of 7) with irregular interval. Interpolation in these results enables 
£,(7) to be obtained for a range of values of 7) with regular interval (see Fig. 1). 


() (Q)= [* hax 10-” 


° u 
from the results of Section 4.6 (a) and the values of ky =k,(0, P,) obtained from 
tables provided by Dr B. Strémgren (unpublished). The excitation potentials V 
of the lines considered are listed in Table II, t9=79(@) is determined by the 
model, and the nitrogen partition function u is obtained from Claas (5). The 
function h,(7)) x 10~-"°/u may then be calculated at regular intervals of 7) (see 
Fig. 1) and the integration (Q) determined. The values of (Q) obtained for the 


dr). The parameter h,(7)) =g,(79)/Ap is obtained 














— h,(7») x 10— V8 Cnn, Uno Fe 
u “ 
o—o—o— 9 (To) —e—e— = 9,(7) 
4.04 
15 
ole _-CLAAS MODEL 
110" x10" a 
30 _-VITENSE MODE! TR HER 
1.0 
2.05 
0. 
(0 
L | | ~J, 
0.5 1.0 1.5 2.0 0.5 1.0 1S 20 


zr 
° (d= 5000A) (ROSSELAND MEAN OPTICAL DEPTH) 


Fic. 1.—The contribution to the total nitrogen absorption as a function of optical depth in the 
atmosphere by the method of weighting functions. (8=10°33 volts ; v refers to A=8206A.) 
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excitation potentials of interest, for each of the models, are shown below: 


TaBLe IV 
Model (Q) 10°33 volts (Q) 10°68 volts 
Claas 1°768 « 107!” 0°923 X 10°19 
Vitense 1°398 x 107!” 0°728 x 107!” 
Swihart 0921 X 10719 0°471 x 10719 


(c) (Y)=5:28 x 10°. A*fg;. All the parameters are given in Table III. We 
may then calculate W,/A =(Y)(Q) and hence determine the relative abundance 
A from each of the nine nitrogen lines, since 


log A = log Wag, — log (Y)(Q). 
The final results are given in Table V. 


TABLE V 
Nitrogen abundance by weighting function method 
log (Q) log A 
A log Winn, log (Y) Claas Vitense Swihart Claas Vitense Swihart 
7442'28 0°4771 14°1920 10'2475 10°1456 11°9645 4°038 4°140 4°321 
7468:29 © 0°5441 = :14°3678 10°2475 10°1456 11°9645  5°929 4°031 4212 
8216:28 07782 14°7333 10°2475 10°1456 11°9645 = 5°797 5°899 4:080 
8629°24 0-6021 14°6546 11'9651 11°8622 11°6729 5°982 4°085 4:275 
8680:24 =1'2304 ~—:150371 10°2475 10°1456 11°9645 5°946 4:048 4:229 
8683-38  o'8451 14°7571 10°2475 10°1456 1179645 = §°841 5°942 4°124 
8703°24 0°6990 14°3564 10°2475 1071456 11°9645 4°095 4°197 4°378 
8711°69 ~=©06021 14°4649 10°2475 10°1456 11°9645 5°890 5°992 4°173 
8718:82 06021 14°3903 10°2475 10°1456 11°9645 5°964 4:066 4:247 
Mean log A=5-942 4:044 4'227 
(on the basis log Ay = 12-00) 
Claas model Vitense model Swihart model 
log Ay=7°94 log Ay =8-04 log Ay =8:23 


The results show a substantially lower nitrogen abundance than those 
obtained by Unsdld (19) (log A, =8-61) and by Hunaerts (8) (log 4, = 9:02). 
The method of weighting functions is however subject to several uncertain 
factors which will be examined in Section 6. In the following section a second 
method for the determination of the nitrogen abundance is developed and, in 
Section 6, a comparison is made between this second method and the method of 
weighting functions. 

5. The method utilizing the Planckian gradient.—The purpose of the present 
section is to develop a method of determining the equivalent width of an 
absorption line as a function of the relative abundance of the element, the 
characteristic parameters of the particular atomic transition, and other factors 
which, in contrast with the method of weighting functions, are determined 
purely by the physical characteristics of the model solar atmosphere considered. 
In this method, suggested by Dr B. Strémgren, no recourse is made to the 
observational results of limb-darkening. 
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The intensity of radiation in the continuum, at the centre of the disk, is 
given by 











Ik(o,0)= { B,(r)e-"» dr, (5-1) 
i.e. . 
I5(0, 0)= ~ | B,(r)d(e") 
0 
~ dB, 
“ |B. 1 +f meer, 
-© dB 
B,(t):, otf a dr, (5.2) 
Again the intensity in the line is 
1,(0,0)= {- B,(r)e~tedt,, (5.3) 
where , 
t= \° (s+ ‘ zt), =T +7, (5.4) 
0 k, . 
and 
n= | aa (5.5) 
/ . k, v° 5-5 J 


We may develop equation (5.3) in the following manner : 


1,(0,0)= ~ |" B,(r)d(e) 











a ” dB, 
= »—t, 
[Bere + a 
and introducing equation (5.4) we have 
°dB 
= ——” e-AT y +7) 
T,(0,0)=B, (tr), -0 + |. Be trndr, 


°dB 
=B,(r),.0+ | FteTedr,(1~7), 
“0 v 





to the first order, i.e. 
T,(0, 0) - Biz), =() a 











oa) _ ‘ — = —tv7, dt, (5.6) 


“0 dr 


From equations (5.2) and (5.6) we see that the dip in the continuum is 


v 





AI,=1,(0,0)~1,(0,0)= [°F e~en dr, (5.7) 


which, by equation (5.5), becomes 


“4B, 
‘ L dr, |d 
Al, = i. = anf; r, |ar, 7, (5.8) 


k, 


The equivalent width of the weil line is therefore given by 


[ cage ened 


« 17(0,0) [- B,(r)e-*» dr, 
0 


a 





(5-9) 
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All factors, other than /,, are of course sensibly constant over the width of the 
line and therefore, utilizing the development of Section 4.3, equation (5.9) 


gives 
p22 Ty he 
\, ° dB, . ty 7é “ all Ni az, | dr, 
W,= dr, mc? my\Jo Nuk, : 


Tae 
0 
Since N,/N,=(N;,/N)A, where the relative abundance A=N/N,,, we may 


introduce Boltzmann’s equation (cf. Section 4.4), i.e. 











(5.10) 











N; _ gx 107"? 
N _ 
and obtain the equivalent width in the form 
° dB, er | ‘hi 107! 5 ar, Jar, 
=e wf o a, \, uk, (5.11) 


ioe] 
mc my | B,(r)e-*» dr, 
0 
For ease in numerical integration we again transpose the integral to a function 
of optical depth 7) in the standard wave-length A= 5000A by the relationship 
dz ,=(k,/Rko) dt» (cf. Section 4.2). As in the method of weighting functions, 
if we express the equivalent width in mA and the wave-length in Angstrém 
units the value of the numerical constant is 5-28 x 10%, and the final form of the 


equation becomes 
re 00 dB, T (70) _ (i) 
“ai dry | dt. 
~~ 


\° B,(r)e-"» os 
0 0 











W,=5'28 x pe | (5.12) 


For a given atomic transition (A, f,g;, V known) the equivalent width is therefore 
expressed as a function of the relative abundance A and the characteristic 
T)=7 (9) of the particular solar atmosphere under consideration. 

5.1. Numerical calculations.—In the wave-length range under consideration, 
i.e. 7442-8719 A, it was again found to be sufficiently accurate for the determination 
of k,=k,(0,P,.) and B,=B,(@) to adopt the wave-length A=8206A as repre- 
sentative of the range. ‘Tables of the Planckian function and absorption 
coefficients were provided by Dr B. Strémgren (unpublished). 


ated. Mided — , ; , 
(a) I(t))= j. ~ dz. The excitation potentials V are given in 





Table II and values of the partition function u=u(6) are taken from Claas (5). 
The value of the integrand was then computed for a large range of values of 7) 


with regular interval and the integration performed at each value of 7) throughout 
the atmosphere. 


(4) (2Z)= 


known as heites of the optical depth 7), the regular interval in t) could be 
maintained thereby greatly facilitating the integration. For a given wave-length, 
in this case 8206 A, the Planck function is a function of temperature only, i.e. 
B,=B,(@), and, as the atmospheric model provides 7) = 7)(@), 


dB, |dt = (dB,|d9)(d0/d7,) 


.— 





e~**I(t))dt». Since all the factors of the integrand are 
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may be calculated in the following manner: 
2hv8 











B, a (e*—1)7, 
where ie Pd — pO, 
1.€ Be = —Be(e—1)9 
and = =f, 
i.€. dB, dB, da — BB, 


d0 da “dO (1—e7*) 

(the constant f has the value, for A=8206A, B=hv/KT@=3-4795) and d@/d7, 
may be determined directly from the model. Values of the integrand obtained 
are illustrated in Fig. 2. It was found necessary to extrapolate graphically the 
integrand to zero at large optical depth; any reasonable change in the extrapolated 
values, however, produces little change in the result of the numerical integration. 

(c) If(o,0)= [- B,(r)e77» re dt). Only in the case of the Swihart model 

0 0 

is the atmosphere defined to sufficiently great depth to allow exact evaluation 
of the.integral, graphical extrapolation of the integrand is necessary for the Claas 
and Vitense models, but again the numerical uncertainty is small (Fig. 2). The 
values obtained for these integrations are shown below. 
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Fic. 2.—The contribution to the total nitrogen absorption as a function of optical depth in the 
atmosphere by the Planckian Gradient method. (@=10°33 volts; v refers to X=8206 A.) 
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Taste VI 
Model (Z) erg cm-* sec! IS(o, 0) erg cm~? sec! 
10°33 volts 10°68 volts 
Claas 124°69 x 10716 67°60 x 10716 4°172 x 1075 
Vitense 96-60 x 10716 52°41 x 10°16 4°049 x 107° 
Swihart 65°27 x 10716 34°44 X 10718 3°605 x 10°° 


(d) (Y)=5-28 x 10°A2fg;. All the parameters are given in Table III. We 
may then calculate W,/A=(Y)(Z)/I¢(0,0) and hence determine the relative 
abundance A for each of the nine nitrogen lines since 


log A = log Wi oy,. 108 (Y)(Z)/L5(0, ©). 


The final results are given in Table VII. These results show an appreciably 


Tasie VII 
Nitrogen abundance by Planckian gradient method 
log (Z)/K(0, 0) log A 
A log W, Abs. log (Y) Claas Vitense Swihart Claas Vitense Swihart 
7442°28  0°4771 = 14"1920 10°4755 10°3777 10°2579 5810 5°907 4027 
7468:29 0°5441 = :143678 10°4755 10°3777 10°2579 5°701 5°799 §°918 
8216°28 07782 =: 147333 10°4755 10°3777 10°2579 5°569 5°667 5-787 
8629°24 0°6021 14°6546 10°2095 I0°I1I20 11-9801 5°738 5§°836 5°967 
8680°24 1°2304 15°0371 10°4755 10°3777 10°2579  5§°718 5°816 5§°935 
8683°38 08451 14°7571 10°4755 10°3777 10°2579  5°613 5°710 5830 
8703'24 o'6990 1473564 10°4755 10°3777 10°2579 5°867 5°965 4:085 
8711°69 = 06021 14°4649 10°4755 10°3777 10°2579 5°662 5°760 5879 
871882  o-6021 14°3903 10°4755 10°3777 10°2579  5°736 5°834 5°954 


Mean log A=5'713 5°810 5§°931 
(on the basis log Ay = 12-00) 


Claas model Vitense model Swihart model 
log Ay=7°71 log Ay=7°81 log Ay=7°93 


lower nitrogen abundance than that obtained by the previous method of weighting 
functions (cf. Table V). The reasons for this numerical discrepancy will be 
considered in the following section where a comparison between the two methods 
is made. 

The accuracy of the results would obviously be increased if more detailed 
profiles of the faint nitrogen lines were available, as measurements using the 
Utrecht Photometric Atlas are difficult and rather uncertain. In spite of this 
difficulty it is pleasing to note (Table VII) that results for each of the nine nitrogen 
lines give only a small spread about the mean value. 

6. Comparison of methods used.—In the Planckian gradient method (Section 5) 
it has been established that for a weak line in pure absorption the dip in the 
continuum is given by the equation 


°dB, —[ rl, 
Al, = i. 7° lf ar. | ar. (6.1) 





Minnaert, in the development of the theory of weighting functions (Section 4), 
obtained the equation 
ro) 2 
Al, = [ S, , €2, tee] er itdr, (6.2) 
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It has been found however that a large numerical discrepancy exists between the 
two methods, the former method using purely the characteristics of the model 
atmosphere, and the latter method utilizing the observational results of limb- 
darkening. We seek then to establish the identity of the two equations and then 
to consider the source of the numerical discrepancy 


io ee 8 ? aT c* é, ak 
Oa ee ee 


0 v v 


| [ea fen] + | nl 
a 




















dr, \|J, k, , 9 a7, k, 
ihe Jene 
“laa ee oe 
-f- i a Fei te lergear, (6.3) 


Again, in the establishment of equation (6.2) Minnaert obtained the equation 
(previously established by Unsdld) 


Al, = |" [ ” Bien» dr,— Be | de (6.4) 
0 Ty k, 





(see equation (4.8)); this, of course, should also reduce to the form (6.1), 1.e. 


0 v 1, a7, 
7 fT [- Feerras, | ear, 
0 Ty Ty v 

= : ‘tee 2 ws "ae. TF rrt, 
= IL 7 dr, e ax, |[-4+ I tae, || + dr, e | ' az, | dt, 
“dB, 
+Ay =. dr, 
is [ een [year, |r. (6.5) 


o dr, ok 


The equations are equivalent. We must therefore look to the information used 
in the numerical integrations for the source of the numerical difference in the 


results obtained by the two methods. 
Minnaert, in the method of weighting functions, considered the precise form 


of the Planckian function, i.e. 


ao 
Be + | 
T 














v 


B,(7,)=a, +b,7, +¢,7,2+d,7,3+e,7,4. 


However, in introducing the observational results of Chalonge and Kourganoff 
only a,, b, and c,, i.e. the first three coefficients, could be determined. In effect 
the form of the Planckian function used is then 


B,(7,)=a,+6,7,+¢,7,?, (6.6) 
i.e. B,/ =b,+2¢,7, (6.7) 
and B,” =2¢,, (6.8) 
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The equation (6.2), representing the dip in the continuum, therefore reduces to 





© /dB a?B l 
= v v —Ty ¥ 6. 
os j, (= . 7i)e * asia (6.9) 
= [- (b,+2c,7,+2¢,)e""» pdr, (6.10) 
0 v 


It may easily be shown that for the special case of the Planckian function of the 
form (6.6) the equation (6.1), used in the method of the Planckian gradient, 
also reduces to the form (6.10). Minnaert has then assumed the first differential 
of the Planckian function to be a linear function of optical depth, and the 
coefficients 6, and c, determining the differential are obtained from observational 
limb-darkening. The values of these coefficients are reliable only for the outer 
layers of the atmosphere and therefore linear extrapolation of the first differential 
to large optical depths (where the contribution of the undetermined coefficients d, 
and e, would be appreciable) is justifiable only if the contribution of the lower 
layers to the absorption is small. The weight of the deeper layers is in fact small 
for all elements of low ionization potential. For nitrogen however the weight 
of the deeper layers is large (maximum contribution at 7) ~ 1-6 and approximately 
40 per cent of the total absorption originates at 7) > 2-0) due to the high ionization 
potential of the element. (See Fig. 2.) 

In the case of the method of the Planckian gradient the precise form of the 
Planck function is used,-i.e. 
2hv* I 
B= =F oar" 
the atmospheric model provides +,=7,(@) and hence dB,/dr, may be computed 
directly as a function of optical depth. The value of the first differential is therefore 
determined entirely by the temperature dependence on the optical depth, which 
is characteristic of each atmospheric model considered. 
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Fic. 3.—Comparison between dB,|dr, as predicted by the solar atmospheric models and as 
determined from observed limb-darkening results. 
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A comparison is made in Fig. 3 between dB,/dr, calculated in accordance 
with equation (6.7) and dB,/dr, as predicted by the atmospheric models (in this 
latter case dB,/dr,=(dB,/dt))(Ro/k,) and the values are then reduced to the 
same scale, dB,/dr,+I$(o,0), used in the method of weighting functions). 
It is seen that there is quite close agreement in the upper layers, except immediately 
below the surface where the weight is small, and the difference is not important 
in the range o<7,<1°5. For 7,>1-5, however, the difference is large and the 
divergence becomes important. It will be noted also that the method of weighting 
functions assumes the second differential constant, equation (6.8), and for 
\=8206A the limb-darkening results give d?B,/dr,2= —0-172 and therefore 
equation (6.9) is equal to zero for dB,/dr,=0-172, i.e. at approximately 7,=2°8 
(t)=2°0). By the method of weighting functions therefore, no contribution to 
the total absorption is possible below t)2-0, while the method using the precise 
form of the Planck function shows that approximately 40 per cent of the total 
absorption, in the case of nitrogen, originates below this level. 

It may be concluded that in cases where the absorption originates almost 
entirely in the uppermost layers, i.e. for elements of low ionization potential, 
the method of weighting functions may be used and may, in extreme cases, be 
superior to the method of the Planckian gradient (since the atmospheric models 
are rather uncertain very close to the surface). In the case of nitrogen, however, 
and in general for elements of high ionization potential, the weight of the deeper 
layers is large and the method of weighting functions is quite unreliable. Here 
the method of the Planckian gradient should be used. For these reasons, in the 
present investigation for nitrogen, the results obtained by the method of weighting 
functions are discarded. 

Since now additional weight is placed upon the results of the Planckian 
gradient method it is obviously desirable to consider equations whereby these 
results may be checked. This is done in the following section. 

7. Equations for use in numerical checking 

(A) Reconsidering equations (5.3), (5.4) and (5.5) we see that we may write 
the intensity in the line in the form 


T (0,0) = I B (r)e*» o(1 + l, i) 
0 k, 
= | *By(r)erefx — (1 e-")] (14 ‘) dr,. (7.1) 
0 » 
Since the intensity in the continuum is given by 
I5(0,0)= | B,(r)e-*rdr,, 
0 
the dip in the continuum may be written 
Al, =I¢(o0, 0)—I,(0, 0)= [- B,(r)e-*» [ _ z +(1-—e77) (: + ¢)| dr, (7.2) 
0 v v 
From equation (5.5) we have tol 
cS ae | % dr,, 
and equation (7.2) beco mes 0 
AI,= [Bune - 2 +| ie ( 
0 


= ‘went 4 : 
\. (r) ; J, 


i) ar Jar, 


ah 
dr, | dr,, (7-3) 
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to the first order. We now require to show this is in fact an alternative form of 
the equation used in the method of the Planckian gradient (equation (5.8)), 
i.e. equation (7.3) may be written 


Al,= - Sg. e+ [fae Jae) 
--|B 1-z + Gear r Jen” + [OS ~p+f Bede, |eredr, 
0 v 
ji «005... 
+ B, k, rarely ei 


k, 
l, i= l eh ae | 
— - ge —trd 
“e : elt =( J, dr, i k, a i k, |e i 
B# 
7 


sae “dB l 
e777 aa ane eae —trd 
(S). dr, 2.48 | +i << 7, 
hes d [l —— 
+ [Bz (Je de, 


. {- oe f . de fer ra (7.4) 
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An attempt was therefore made to use equation (7.3) as a numerical check on 
the accuracy of the computations made in the method of the Planckian gradient. 
Comparing equations (5.9) and (5.12) we find 


Are) J T (Te) 10-V8 
I. ho dr» «| az dt» 

and /,/k, is therefore easily obtained since /,/k,=(1,/ky)(Ry/k,), all other factors 
have been previously computed and hence equation (7.3) may be evaluated. 
A feature of the computations is the negative contribution to the absorption in 
the upper layers which subtracts from the ‘‘ integrated ’’ absorption in the lower 
layers. Here the contribution from the very deep layers is most important and 
it was found that the method could only be applied in the case of the Swihart 
model, which is defined to great depth. The equation (7.3) is moderately 
insensitive to small variations in the parameters and should provide a good 
numerical check for models defined to sufficiently great depth. 

(B) Since the models of Vitense and Claas are not defined to sufficiently 
great depth to allow the application of equation (7.3) recourse was made to Unsdld’s 
equation (equation (4.8)) 


Al, = [elf Bera Be | de. (7.5) 
0 k, Ty ji 


It has been established that this equation is equivalent (equation (6.5)) to the 
form (5.8) used in the method of the Planckian gradient. It was found however 
that the results obtained by equation (7.5) are extremely sensitive to small 
variations in the values of § P 
e 


Ty 
entedr, and | B,e-dr,. 
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The results obtained by the above methods are compared with the Planckian 
gradient results in the following table. 


ie 
wh 
gra 
sen 


an 
ex 


le: 








XUM 


118 No. 5, 1958 " the determination of solar atmospheric abundances 455 


: of ‘Tasce VIII 


3)), Model Planckian Gradient Check (B) Check (A) 
method method method 
Claas (Z) = 124°7 x 10716 126°9 x 10718 
Vitense (Z)= 96:6 x 10716 86-8 x 10716 aa 
Swihart (Z)= 65-3 x 10716 64°0 X 10716 61:8 x 1071* 


Equivalent width W,=const. (Z)A. 


It must be emphasized however that only in the case of the Swihart model, 
where the model is defined to large depth, are the numerical results of the check 
methods reliable, and even in this case are not to be preferred to the Planckian 
gradient calculation (i.e. in the check (B) method the results are extremely 


sensitive to small variations in | B,e~*» dr,, and the check (A) method shows 


an appreciable contribution even at very large depth). Again the apparently 
excellent agreement obtained for the Claas model can only be regarded as 


me) 
extremely fortuitous in that a reduction of o-oro in the value of | B,e~*» dr,, 


0 
which is easily obtained by a reasonable variation in the extrapolation of Fig. 2, 
.4) reduces the value of (Z) to approximately 108 x 10716. For the Vitense model 
the agreement is only moderate but here again a similar increase of o-oro in 
“0 
on B,e~*» dz, increases the value of (Z) to approximately 102 x 10~'®. We may 
nt. 


cnnibuie therefore that, for the Swihart model, the agreement between the 
Planckian gradient result and the results obtained by the check methods is very 
good. For the Claas and Vitense models the Planckian gradient results are well 
within the degree of error possible in the results of the check methods. 

8. Discussion of results—The method of weighting functions developed by 
d Unséld and Minnaert leads-to the results, for the nitrogen abundance, shown in 
‘Table V. The method of the Planckian gradient developed in the present paper 


or leads to the substantially lower values of the nitrogen abundance (Table VII): 
ad Claas model Vitense model Swihart model 

a log Ay =7°71 log A, =7°81 log Ay =7°93 

~ (on the basis log A,, = 12-00) 

The comparison of the two methods, as presented in Section 6, leads to an 
ly argument which, in the opinion of the author, is sufficiently in favour of the method 
’s of the Planckian gradient, in the case of nitrogen, to necessitate the exclusion of 

the results obtained by the method of weighting functions. 

Although the above results indicate that the derived nitrogen abundance is . 
5) moderately insensitive to the atmospheric model chosen, we would require, in 

estimating the most probable value of log A,, some measure of the reliability 
1€ of each of the models. This will be considered in a subsequent paper. 
or (i) Previous determinations.—Hunaerts (8), in 1950, determined the abundance 
I] of nitrogen from measured intensities of rotational lines in the NH molecular 
band and obtained the value log A,,=9-02. The molecular abundance has first 
to be determined by a curve of growth method and this requires a knowledge of 
the oscillator strengths of the molecular transitions. A great deal of doubt 
n however still exists with regard to these f-values, and Hunaerts points out that 
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the values adopted may be incorrect by appreciable factors. The atomic 
abundance may now be determined by consideration of the dissociation 
equilibrium, but here additional uncertainties arise due to the inadequate data 
concerning the dissociation energies. Hunaerts concludes that the estimated 
atomic abundance is very uncertain. 

Unséld (19), in 1948, in an extensive analysis of the solar atmosphere, 
obtained the value log A, =8-61 by a curve of growth method. Slightly low 
values of ionization potential were adopted for the five lines considered although 
the error introduced is probably slight. The method does however suffer 
from the fact that the detailed structure of the atmosphere is not considered and 
an idealized atmosphere is assumed. 

(ii) Comments.—A major source of uncertainty in the abundance calculations 
is the measurement of the equivalent widths of the weak nitrogen lines as 
presented in the Utrecht Photometric Atlas. However, since the nitrogen 
abundance, as determined for each of the nine lines, shows only a small spread 
about the mean value (Table VII), it is probable that this restriction is not too 
severe. When more definite measurements become available, simple 
recalculation in Table VII will provide any necessary correction in the final 
estimate of the nitrogen abundance. 

The hydrogen/metal ratio, o, in each of the models has been taken to be 
log.  =3°8. This is probably rather low and a value of 3-90 or 3-95 may be a 
better estimate. Further, the helium abundance in the atmosphere has been 
ignored but, as shown by Claas (5), abundance determinations will be rather 
insensitive to this factor. 

It may be considered to be a disadvantage of the method of the Planckian 
gradient that no recourse is made to the accurate determinations of the limb- 
darkening of the Sun, and that consequently great stress is placed upon the 
accuracy of the theoretical model. It must be pointed out however that the 
Claas and Vitense models are themselves based upon the observed limb-darkening 
results in the regions where these results are most reliable. The Swihart model 
however is more at variance with the limb-darkening results. 


The author would like to thank Dr B. Strémgren for suggesting this problem 
and for the time he most generously spent in many helpful and stimulating 
conversations during the author’s visits to the Yerkes Observatory. The author 
is also indebted to Dr Swihart for permission to utilize his results in advance 
of his own publications and to Professor C. W. Allen for the interest he has 
shown in the progress of this research. 


Department of Theoretical Physics, 
The University, 
Manchester 13 : 
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EXCITATION TEMPERATURES FOR THE CN EMISSION FROM 
THE LOW CHROMOSPHERE 


D. V. Thomas 
(Received 1958 May 2) 


Summary 


Microphotometer tracings have been made of the region of the CN 3883 A 
band system, which is a prominent feature of the flash spectrogram obtained by 
Redman at the 1952 eclipse. The CN bands are seriously affected by 
blending, but enough unaffected lines have been found to enable the rotational 
intensity distribution in the o—o band to be measured at three heights below 
about 400 km. Self absorption is important, particularly at the band heads, 
and its effects have been estimated. ‘The measurements indicate a negative 
excitation temperature gradient, the temperature at about 400 km being 
~ 4500°, and the value of the gradient uncertain. However, the CN emission 
may not be inconsistent with the positive kinetic temperature gradient of 
de Jager’s most recent chromospheric model, if departures from local thermo- 
dynamical equilibrium and the imperfection of the observational data are 
taken into account. 





1. Introduction.—Until quite recently it was thought that there was strong 
evidence both for a high ( ~ 30 000°K) and a low ( < 10000°K) kinetic temperature 
for the low chromosphere. The conflicting evidence has been summarized and 
discussed by van de Hulst (1953). While a low temperature is now generally 
accepted, its exact value and gradient are still uncertain. One of the most 
prominent features of the flash spectrogram obtained by Redman at the total 
eclipse of 1952 February 25 is the 3883 A band system of cyanogen. In principle, 
these bands supply an excellent means for obtaining accurate temperatures, as 
the relative intensities of the lines composing the bands and of the bands themselves 
depend strongly on the temperature. However, in practice, measurements of the 
intensity of the CN emission in the flash spectrum yield directly only the relative 
populations of different energy levels of the molecule, and hence an effective 
excitation temperature. The corresponding kinetic temperature is then a matter 
of interpretation in terms of the physical conditions assumed to exist in the 
chromosphere. 

Determinations of molecular excitation temperatures in the low chromosphere 
have already been published by Blackwell (1955), Pecker and Athay (1955) and 
Parker (1955). Using Redman’s 1940 eclipse spectrogram (Redman 1942) 
Blackwell measured lines of the weak 4216 A CN and 4300 A CH bands, neglecting 
the bandheads. Pecker and Athay’s discussion is confined to the 3883 A CN band 
head on the High Altitude Observatory’s 1952 slitless spectrograms, and is open to 
serious criticism. Parker, working with the Mount Wilson 75-ft spectrograph 
outside eclipse, had the advantage of a high resolution, but needed to employ a 
difficult correction for the contribution of photospheric radiation (both scattered 
and direct) to the intensities of the Swan C, bands. These authors’ results are 
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discussed in the light of the present work in Section 7. The excitation tempera- 
tures determined in this paper are based upon the largest dispersion slit spectro- 
gram yet obtained at an eclipse. In the analysis, both band heads and individual 
lines have been considered and allowance has been made for self-absorption. 

2. Observational material.—Only the essential details of the instrument used 
and the procedure followed at the 1952 eclipse need be noted here, asa full account 
has already been given by Redman (1953, 1955). ‘The slit spectrograph employed 
a 15000 lines per inch 6-in concave grating of 21-ft radius of curvature in a Wads- 
worth mounting. The resulting dispersion in the wave-length range 3400A 
to 4100 A was 2-4A/mm. The instrumental profile was determined at the eclipse 
site by means of a mercury isotope lamp. The excellent flash spectrogram used 
in the present work was obtained during an exposure of 1-5 sec immediately 
preceding third contact, the coelostat being rated to follow the Moon. An 
attempted relative intensity calibration on the eclipse plate itself was unsuccessful. 
Four photospheric spectrograms were therefore obtained on another plate from the 
same batch with a Hilger rhodium-on-quartz step wedge in front of the slit. The 
two plates were processed identically and may be assumed to possess the same 
characteristics. ‘The transmission factors of the step wedge were determined at 
the National Physical Laboratory. 

A portion of the flash spectrogram in the region of the principal CN band heads 
is reproduced in Plate7. The variation in intensity perpendicular to the direction 
of the dispersion is due partly to irregularities on the limb of the Moon and dust 
on the jaws of the slit, partly to the fact that positions along the slit represent, in 
some measure, heights above the photosphere, and partly to the spicular structure 
of the chromosphere. _ Initially, it is sufficient to denote positions along the slit by 
the arbitrary system of ‘‘levels ’ indicated in Plate 7, as the analysis of each level is 
independent. Only whencorrelating the results from different levels is it necessary 
to know even their relative heights, which are easily obtained. ‘The determination 
of the actual effective chromospheric heights is a difficult problem, but they are 
necessary only for the final interpretation of the results. During the exposure the 
slit, following the Moon, scans the chromosphere at a rate of about 400 chromo- 
spherickm/sec. ‘The ‘‘ seeing ’’ was estimated by Redman to be 2”, corresponding 
to 1400 chromospheric km. Asa result, the recorded emission originates in varying 
proportions from slices right through the chromosphere, covering a range of 
about 2000km of the chromosphere and the (occulted) photosphere. A geo- 
metrical determination of effective heights is therefore impossible. Redman and 
Suemoto (1954), in their study of the metal lines in this spectrogram, compared the 
relative intensities of some prominent lines with their relative intensities on the 
slitless spectrograms obtained by Houtgast (1953) at the same eclipse. In this 
way they estimated that Levels 1, 2 and 3 correspond to heights < 50km, < 1ookm 
and ~6ookm, respectively. According to Houtgast the intensity gradient in the 
chromosphere is much steeper for the CN 3883 A band head than it is for the 
metallines. The effective heights of the CN emission must therefore be somewhat 
less than the above figures. The adopted heights for the three levels, of which an 
analysis has subseqently been made, are: Level 1~okm; Level 2~50km; 
Level 3 ~ 400km. 

In order to obtain the intensities of the CN lines in the flash spectrum, tracings 
of Levels 1, 2 and 3 were made with the Cambridge Observatories’ recording 
densitometer. Levels 2a and 3a proved to be unsuitable because of the weakness 
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of the CN lines and Levels 1a and rb because of their high background intensities. 
The characteristic curve of the plate was obtained from tracings made across the 
step wedge spectra, no change with wave-length being detectable. A mean 
characteristic curve for the microphotometer tracings was then derived, the range 
of measurable intensities being represented by 100 units on an arbitrary relative 
intensity scale. Fluctuations in the dark current and clear glass transmission, etc., 
produce an estimated maximum error of + 2 units in intensities obtained from this 
mean curve. ‘The microphotometer tracings were reduced to direct intensities 
at the Dominion Astrophysical Observatory, B.C., using the Beals (1944) intensito- 
meter. The completed records have an amplitude of 1oin., corresponding to 
100 units on the above scale. At various points they were checked with a point- 
by-point reduction and found to be very satisfactory. The record of Level 1 in 
the region of the CN o-o band head is reproduced in Fig. 1. 
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lic. 1.—Reproduction of the region of the CN 0-0 band head on the direct intensity record of Level 1. 
Theoretical positions of individual lines are indicated below. 














From measurements of the positions of 75 sharp lines, whose identifications 
and wave-lengths were obtained from the Revised Rowland Table of photospheric 
wave-lengths (St John 1928), the dispersion of the plate was found to be linear 
over the range covered by the CN bands. The same lines were identified on the 
intensity records, the dispersion obtained from their measured positions being 
82-86 mm/A, showing the magnification of the records over the plate to be just 
less than 200 times. The maximum error in the position of a line on the records 
was estimated to be + 4mm (+0:°05 A). 

3- The ultra-violet CN bands.—The ultra-violet bands of cyanogen arise 
from the transition between the X?X+ and B®X+ electronic states of the CN 
molecule. The sequence Av=o forms a series of bands degraded towards the 
violet with its head at 3883 A, the head of the o-o band. In 2% electronic states 
the rotational levels with quantum number K are split into two components 
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with quantum numbers J/=K+4. Asa result of the selection rules AK= +1, 
AJ =o, +1, a 22-*Z band consists of a doublet P branch, a doublet R branch, 
and two satellite P-form and R-form QO branches. The satellite branches are 
very much less intense than the P and R branches, so the CN bands are often 
treated as 'X—'X bands (.J=K), each P and R branch doublet being considered 
as a single line at the centre of the doublet. The relations governing the wave- 
numbers of the individual lines in a 'X—'X band are given in the standard texts 
(e.g. Herzberg 1950). The P and R branch relations may be expressed as one 
equation with the “‘ running number ’’m as argument: 


a(m)=A+Bm+Cm? + Dm*+ Em*+ ..., (3-1) 
where for the P branch m= — K”, | 
and for the Rbranch m=K"+1, | 


J 


(3.2) 


and the constants A...F are easily determined functions of the rotational 
constants of the molecule. The most recent rotational constants, due to Douglas 
and Routly (1955), have been adopted. The values of the constants are such 
that for the most important bands of the 3883 A band system, o(m) is a minimum 
at m= — 30, so that a band head is formed in the P branches near P (30). 

The most extensive laboratory measurements of the CN 3883 A band system 
are still those made by Uhler and Patterson (1915). ‘There is no systematic 
difference between their wave-lengths and those recently published by Weinard 
(1955): as the random differences are generally less than 0-005 A, which is only 
one-tenth of the possible error in the intensity records, the accuracy of the 
published wave-lengths is more than sufficient for the purpose of identifying CN 
lines on the records. 

For each of the o—-o...4~—4 bands, the gap corresponding to the forbidden 
line P(o) has been detected, so the quantum numbering of the R branches, and 
the P branches as far as the band heads, is simply obtained by following the 
regular series of lines. The numbering of the returning portions of the o-o 
and 1-1 P branches has been obtained by Smit-Miessen and Spier (1942); 
they compared Uhler and Patterson’s measures with values computed from 
equation (3.1), whose coefficients were determined from the wave-numbers of 
lines already indentified. The identifications have been confirmed by computing 
the coefficients in (3.1) independently from Douglas and Routly’s band constants. 
The P branches of the other bands are too greatly obscured in the flash spectrum 
for identifications to be necessary. 

The positions of all Uhler and Patterson’s lines whose quantum numbering 
had been determined were marked on the intensity records. The general 
concentration of lines in the ultraviolet region of the flash spectrum is so great 
that even the Rowland table of photospheric wave-lengths does not contain them all. 
It cannot be determined with certainty, therefore, that any particular CN line 
is unblended. For this reason, and also because of the irregularities in the 
dispersion of the records, certain identification of lines on the records as CN lines 
is not possible. 

4. Intensity distribution in CN bands in emission.—As the flash spectrum does 
not possess an absolute intensity calibration, the present work is concerned only 
with the relative intensities of the lines emitted as a result of the transition 
B*Y--X* at any particular temperature. Assuming a Boltzmann distribution 
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of molecules amongst the energy states corresponding to an excitation tempera- 
ture 7, the relative intensities of lines emitted as a result of transitions from the 
levels (v', K’) of the upper state to the levels (v”, K”) of the lower state are given by 


Y wv , he + re he 
Ly ge — ore = COy ggg Py exp{ 3'(v ery SKK exp —F,(K ae 


(4.1) 
where C =a constant determined by the transition, the temperature and the total 
number of molecules, 

Oy g’— rx’ = the wave-number of the line, 
P.,_ =the vibrational transition probability, 
Sx =the rotational line strength, and 
G'(v'), F,(K’) =the vibrational and rotational term values of the upper state, 
dependent upon the band constants. ‘The relative intensities of the rotational 
lines in a band of the Av=o sequence are therefore given by 


I g-_ g- = constant of _ geSx-_ ge exp Fy (K’) “at (4.2) 





k1 

The line strengths for the various branches of ?X—*2 and !X—'S bands have been 
derived by Mulliken (1927). When K’>5, the intensities of the satellite Q 
branches of a ?2X—*> band are negligible, and the combined intensity of the two 
components of a P or R branch doublet differs by less than 1 per cent from the 
intensity of the corresponding singlet in a 'X—' band. ‘The intensities of the 
CN lines are therefore given to a sufficient approximation by (4.2), where Sx. ,- 
now represents the line strengths for 'X2—'X bands, i.e. 


for the P branch Sy. ,.=2(K’+1) | _ 
and for the R branch Sy. ¢.=2K’, “3 
and oy, and F,,’(K’) now refer to the centres of the doublets. 

A useful modification of (4.2) is obtained by taking logarithms: 


In {wea} = constant — F,/(K’ =. - (4.4) 
Plotting In {/,._ ,-/o*S,- _,-} against F,'(K’)hc/k, where the /’s are measured 
intensities, a least-squares method yields the slope —1/T of the straight line 
relation, and hence the temperature directly. 

The relative intensities of homologous lines in the v-v and o-o bands can be 
obtained from (4.1). Taking F,,'(K’)=¥,'(K’) and substituting the wave- 
numbers of the respective band origins for those of the individual lines (these 
approximations are justified by the lower accuracy of the measures subsequently 
made on the 1-1... 4~—4 bands and the uncertainties in the values of the vibrational 
transition probabilities), the relation simplifies to 

I(v, K) _ o3_, Py, exp{—G'(v) he/kT} ’ 

T(0,K) ~ of-g' Py» exp{—O(o)he/RT) os 
The transition probabilities adopted are those computed by Pillow (1953). 
The small differences between these and the more recent laboratory values of 
Floyd and King (1955) have no significant effect on intensities calculated from 
(4-5). 
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5. Determination of effective excitation temperatures 

5.1. Preliminary estimates.—As the CN lines are strongest on the intensity 
records of Level 1, this level was selected for preliminary measurements. An 
immediate difficulty is the presence of a continuous background caused partly 
by the continuous radiation from the chromosphere itself, and partly by scattered 
photospheric radiation from the first traces of a Baily’s bead (see Redman 1955). 
The background is distorted in many places by the overlapping of large numbers 
of faint lines and the wings of strong lines, and by errors in the photometric 
reductions. Superimposed upon it are fluctuations due to the photographic 
grain and dust on the plate. Photometric errors due to the use of a mean 
characteristic curve in the reduction to direct intensity can be reduced by drawing 
the apparent rather than the mean continuous background, assuming it to vary 
smoothly between those regions of the spectrum apparently not enhanced by 
overlapping weak lines. ‘The background was determined separately for two 
limiting cases: the lowest possible background, which assumes that the major 
small-scale irregularities are weak lines, and the highest possible background, 
which assumes that these irregularities are all ‘‘ noise’’. An analysis based 
on an estimation of the exact position of the background has no advantage over an 
interpolation between the results obtained from the two limiting cases considered. 

It appears from the intensity records of Level 1 (Fig. 1) that the o-o band 
head is probably affected by self-absorption; the peak intensity of the band head 
is very little greater than the peak intensities of the pairs of lines nearby. The 
simplest method of obtaining a preliminary estimate of the temperature therefore 
appeared to be to find the distribution of intensities amongst the rotational lines 
of the o—o band and make use of (4.2). ‘The intensities of all the lines in positions 
where o-o band lines were predicted were measured with a planimeter, assuming 
the lowest possible background for this preliminary determination. The total 
area of each P branch pair (formed by the near coincidence of lines of the proceed- 
ing and returning portions of the P branch) and the total intensity of the band head, 
in three sections, were also measured. ‘The best method of representing the 
measures graphically is that adopted in Fig. 2, the abscissa of which is a linear 
scale in the running number m reflected at m= — 28, the head of the P branch. 
This makes possible the representation of the intensities of the P branch pairs, 
since the nearly coincident components of the pairs appear as exactly coincident 
inthe figure. From m~ —21 to —35 the separation of the individual components 
of the doublets is of the same order as that of the doublets themslves (see Fig. 1), 
so the band head would not exhibit the “pair ’’ formation even with complete 
resolution of the lines. However, as it is only possible to measure the total 
intensity of the band head, the representation is not rendered invalid, and it is 
convenient for comparing measured and theoretical intensity distributions. 

No account was taken of possible blends, so the intensities in Fig. 2 are 
in effect maximum intensities. As perturbed lines were omitted, the plotted 
points should all lie on or above (for unblended and blended lines, respectively) 
the intensity distribution curve for the appropriate temperature. ‘The intensity 
curves for three temperatures, calculated by means of (4.2), are drawn in Fig. 2 
to scales such that all the plotted points lie on or above the curves. Comparing 
the three curves with the measured intensities, it appears that very few of the 
measured lines are unblended. ‘The distribution for 7000° fits the measures 
best but requires considerable self-absorption at the band head. ‘The conclusion 
that the band head is depressed by self-absorption can be avoided by assuming a 
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temperature of 10000°, but this would necessitate large blending corrections for 
all the strong lines and no correction at all for a number of weak lines. The 
effective rotational temperature therefore appears to be in the region of 7000", 
not as high as 10000 , and probably not as low as 4000”. 

An approximate blending correction to the measured intensities was derived 
from consideration of the average intensities of lines immediately to the red of 
the o-o band head. However, the corrected intensities merely confirm the 
indication of Fig. 2 of a rotational temperature in the region of 7000°. 
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Fic. 2.—CN 0-0 band, Level 1: uncorrected line areas (assuming lowest continuous background). 
In this and subsequent figures the points ‘‘o”’ represent the areas (= intensities) of the corresponding 
lines or “ pairs’’ of lines. The area of each rectangle equals the total area (= total intensity) of the 
corresponding group of unresolved lines. 


Because of the suspected depression of at least the o—-o band head by self- 
absorption, it would be unwise to attempt to obtain an estimate of the effective 
vibrational temperature by comparing the intensities of the band heads. The 
P branches are badly obscured, but the intensities of lines in the R branches of 
the various bands may be compared, utilising the relations (4.2) and (4.5). The 
areas of all the lines in positions where the R branch lines of the 1-1... 4-4 bands 
had been predicted were therefore measured, no attempt being made to allow for 
blends. ‘The measured intensities were plotted on diagrams of running number 
versus intensity, together with the intensity distributions corresponding to those 
already drawn in Fig. 2 for the o-o band. While the 7000” curves agree best 
with the measures, neither the curves for 4000° nor those for 10 000° violate the 
condition that all the measured points should lie on or above the curves. 

5.2. Computed CN line profiles.—It is not possible to deduce reliable CN line 
profiles from the intensity records because all the CN lines are unresolved doublets 
and only a very few appear to be free from blends. However, it is possible to 
compute profiles for an assumed temperature and turbulence of the emitting gas, 
knowing the characteristics of the recording apparatus. It may safely be 
assumed that the effects of collision and radiation damping are negligible for all 
except the strongest chromospheric lines. The emitted CN line with wave- 
length A, will therefore have a Doppler profile with equation 


1(a)=Cexp{ ~ (se + a) a0-ni}, (5.1) 


where jz is the molecular weight, 7 is the kinetic temperature, and the probability 
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of a molecule having a turbulent velocity in the line of sight between € and € + dé 
is proportional to exp {— £?/E,?}dé. A half-width of 0-073 A was obtained for 
all CN lines on computing the profile (5.1) with the following assumptions: 
the kinetic temperature equal to the estimated effective excitation temperature, 
i.e. 7000 , a turbulence = 2-7km/sec, as determined by Redman and Suemoto 
(1954) from metal lines at the identical Level 1, A, constant for all CN 
lines. 

The effect of the apparatus function of the spectrograph and the width of the 
microphotometer slit on the emitted line profile was determined by the method 
of Voigt functions described by van de Hulst and Reesinck (1947). The resulting 
computed profile of a CN line on the intensity records of Level 1 has a half-width 
of 8-gmm (0:11 A). The variation in the computed profile with wave-length is 
negligible for the wave-length range covered by the CN bands. The total 
variation in the computed half-width, due to an assumed probable error of +0°5 
km/sec in the adopted &, and to a temperature range of + 3000° about the adopted 
7000", is less than +10 per cent. The computed CN line profiles are thus 
relatively insensitive to the adopted kinetic temperature, and can be used to 
determine the intensities of the CN lines on the records more accurately. 

5-3. Effective excitation temperature of Level 1.—F¥or the reasons pointed out 
earlier, the predicted line positions may be in error by + 4mm, and the fact that 
the majority of CN lines are affected by blending means that in most cases the 
exact line positions cannot be determined more accurately by inspection. ‘Thus 
in each case lines having the computed profiles were constructed within + 4mm 
of the predicted position such that the resulting line intensity was the maximum 
possible consistent with the feature on the intensity record (the two components 
of the doublets were taken to have equal intensities). For those doublets unre- 
solved by Uhler and Patterson, theoretical doublet separations were obtained 
from a diagram of the separations of the resolved doublets plotted against the 
rotational quantum number K”: the mean curve through the plotted points 
was extrapolated back to K” =o in accordance with the relations governing the 
splitting of the rotational energy levels (Mulliken 1930). 

As the half-widths of all the CN lines are the same, their intensities are 
proportional to their central intensities. Hence the height of the central ordinate 
of cither component of a constructed doublet may be used instead of its area as a 
measure of its intensity. ‘The intensities of the constructed lines of the o—o band 
are plotted in Fig. 3 (heavily blended lines are omitted). Individual lines could 
not be constructed in the unresolved portion of the band head, so its previously 
measured area was converted to the new intensity units by determining the area of 
a standard line profile. 

The intensities in Fig. 3 are once again maximum intensities, apart from six 
lines (two of these constituting a P branch pair) which appear to be unblended 
and whose intensities may therefore be considered accurate. A least-squares 
solution of a plot of In {7/o4K’} versus Fy’(K’ )he/k for the three accurate R branch 
lines R(14), R(49) and R(74), gives T=7300°. The probable error, + 200°, 
does not give such a good idea of the likely maximum or minimum temperatures 
as the scatter + 500°, obtained by assuming each pair of points in turn to be 
correct. ‘The computed temperature distributions for 7300° and 7800°, shown 
in Fig. 3, fit the R branch intensities almost equally well, but the band head again 
appears to be depressed by self-absorption. 








466 D. V. Thomas, Excitation temperatures for Vol. 118 


Comparing the corresponding computed intensity distributions in the 
1~1...4-4 bands with the measured intensities of the constructed R branch 
lines, there appears to be no appreciable difference between the agreement for 
7300 and 7800", which is good in each case. Thus there is no detectable 
difference between the effective rotational and vibrational temperatures. 





= mM T=7800° 


0 T7300? — — — 
Or 

6 > + + = ACCURATE MEASURE 
st (NY 

40 \ 


30+ 


INTENSITY 


20 























-28 “20 -10 ° 10 20 3 40 50 6&0 70 8 90 wo 

-% 46 -56 RUNNING NUMBER 
Fic. 3.—CN o-o band, Level 1: intensities of constructed lines (assuming lowest continuous 
background). 
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Fig. 4.—CN R branches, Level 1: intensities of constructed lines (assuming highest continuous 
background). 
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The effect of changing the assumed height of the continuous background was 
investigated by repeating the above procedure with the highest possible back- 
ground level previously determined. A diagram similar to Fig. 3 was obtained, 
with 7’=7000° and a small probable error. The measures of the R branches 
are plotted in Fig. 4, from which it is clear that a vibrational temperature of 7000° 
requires intensities significantly in excess of the measured intensities for the 
2-2, 3-3, and especially the 4-4 band. There are several possible explanations 
of this. First, there may be a considerable difference between the vibrational 
and rotational temperatures, but this should not be inferred from the above 
evidence alone. Secondly, if the assumed continuous background is too high 
then the intensities of the weaker lines and bands will have been reduced 
systematically, compared with the intensities of the stronger lines. ‘Thirdly, if 
the continuous background has been correctly drawn then the vibrational— 
rotational temperature must be less than 7000’. In this case, the measured 
intensities of the o—o band must either all be too great because of blending, which 
seems unlikely, or the intensities of the strongest lines, in addition to the band 
heads, must be depressed by self-absorption. ‘lhe theoretical intensity distri- 
butions for 4000°, drawn in Fig. 4, indicate that even such a low temperature as 
this might be possible if self-absorption is assumed to be sufficiently great. 

It is clear from these results that the maximum possible excitation tempera- 
ture for Level 1 is 7800°, obtained on the assumption of the lowest possible 
continuous background. Raising the continuous background level necessitates 
a progressive lowering of the excitation temperature and an increase in effects 
which must be attributed to self-absorption. A reliable estimate of the continuous 
background level and quantitative knowledge of the effects of self-absorption are 
therefore required to fix the excitation temperature accurately. 

5.4. Effective rotational temperatures of Levels 2 and 3.—The absence of a 
detectable continuous background on the intensity records of Levels 2 and 3 
is a considerable advantage, but because of the weakness of the CN lines, errors 
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Fic. 5.—CN o-o band, Level 2: intensities of constructed lines. 


in the photometry are proportionally greater. ‘lo minimise such errors, the 
zero intensity in the vicinity of each line was estimated from the original tracings. 
The line profile appropriate to Level 2 was then computed and the procedure 
followed for Level 1 employed to obtain the measured intensities in Fig. 5. ‘The 
close agreement with a rotational temperature of 5000” is possibly fortuitous, 
as the band head exhibits a deficiency, probably due to self-absorption. However, 
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the intensity distribution for 7000° demonstrates that a temperature as high as 
this is most unlikely. The lines with m>75 are too weak to be detected on 
the records, as are all the lines of the 3-3 and 4-4 R branches. A few lines of the 
1-1 and 2-2 bands are measurable, but not enough for a vibrational temperature 
to be obtained. 

The CN lines at Level 3 are so weak that it was not feasible to construct line 
profiles, so the procedure used to obtain a preliminary temperature for Level 1 
from the o—o band had to be adopted. No attempt was made to detect lines 
in other bands. ‘The measured intensities, converted to central ordinate inten- 
sity units, are plotted in Fig. 6 together with theoretical intensity distributions 
for 4000°, 5000° and 7000°. The measured intensities are maximum inten- 
sities, so a temperature of 4500° + 500° seems most likely and 7000° unreasonably 
high. As in the cases of Levels 1 and 2 there are indications that self-absorption 
affects the band head. 
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Fic. 6.—CN 0o-o band, Level 3: measured intensities (uncorrected line areas). 


6. Self-absorption.—Prima facie evidence for the presence of self-absorption 
in the flash spectrum is provided by the conspicuous reversal of intensity in the 
centres of many lines of quite moderate intensity, some of which can be seen in 
Plate 7. However, it must be emphasized that self-reversal can only occur 
where there is a negative excitation temperature gradient, whereas self-absorption 
is inseparable from line emission by a gas in local thermodynamical equilibrium 
(L.T-E.), although in an optically thin layer (by definition) its effects are negligible. 
The theoretical treatment of self-reversal is a difficult problem, but fortunately 
it is not directly observed in the CN lines, so it may be assumed initially that the 
chromospheric CN isin L.'T.E. at a single temperature. Even so, a new approach 
to the problem of self-absorption is necessary because, while previous work has 
been concerned with the profiles or equivalent widths of single lines, the CN 
bands consist of doublets with varying separations and unresolved groups 
of lines. 

In the absence of self-absorption, the ‘‘ideal’’ emitted intensity of each 
component of a particular CN doublet will be given by the relations quoted in 
Section 4. Each component will have a profile given by (5.1), and the total line 
profile will depend on the doublet separation. Self-absorption will modify the 
profile of the emitted line. The intensity records show the emitted self-absorbed 
line broadened by the known apparatus function, but unchanged in total intensity. 
In principle, therefore, the ‘‘ideal’’ intensity of each component of a CN doublet 
can be deduced from the intensity records if the relation between ‘‘ideal ’’ and 
‘“ self-absorbed ’’ intensity is known, The “‘ideal’’ intensity of emission from 
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a cylinder of the chromosphere in L.T.E. at a temperature 7 is 

T/=B,(T)r, (6.1) 
where B,(7') is the Planck function and 7, is the optical thickness of the cylinder 
at the frequency v. ‘Taking self-absorption into account, the intensity of emission 


from the same cylinder is 
1,=B,(T)1—exp (—7,)} (6.2) 


(see, e.g. Unsdld 1955). Hence the “‘ self-absorbed ”’ intensity is related to the 
‘‘ideal’’ intensity by the equation 


1,=B(T)t1 —exp [—1,'/B, (T)]}. (6.3) 


Two difficulties prevent the straightforward application of these principles. 
First, the intensity records do not have an absolute intensity calibration, so 
B,(T) in (6.3) is not known in terms of the relative intensity units adopted. 
Secondly, the removal of the effect of the apparatus is a practical impossibility 
in the case of a band head, and even in the simpler case of an unresolved doublet 
it would require mechanical integration, as the instrumentally broadened self- 
absorbed profile is not a Voigt function. However, if a value for the temperature 
is assumed, a theoretical CN band can be computed completely using the appro- 
priate intensity relations of Section 4 and either measured or calculated wave- 
lengths. Assuming a value for the Planck function, the corresponding self- 
absorbed band can be completely computed using (6.3). The areas of the 
observed lines and line groups on the intensity records (with due allowance for 
blends) should be equal to the areas of the corresponding computed lines, provided 
the assumptions regarding temperature and Planck constant are correct and the 
scale of the computed band has been correctly chosen. Hence, by trial and error, 
it should be possible to determine (approximately) whether the CN emission is 
in accordance with a chromosphere in L.'T.E. and, if it is, to find the value of 
the corresponding absolute intensity scale of the records and the effective 
excitation temperature. This procedure was applied to the CN o-o band at 
Levels 1, 2 and 3. 

In order to compute a theoretical band it is necessary to know the frequencies 
of the lines in the unresolved portions of the band head. These have been calcu- 
lated by Smit-Miessen and Spier (1942), and their values have been adopted for 
the whole of the P branch. Doublet separations were obtained as described 
in Section 5.3. A rotational temperature of 4500° was adopted, this being the 
temperature previously derived for Level 3, where self-absorption must be least, 
and also the temperature obtained by Hunaerts (1947) from photospheric CN. 
The turbulence adopted was that appropriate to Level 1. The computed band 
need not be very accurate because the intensity records are subject to errors of 
many kinds already mentioned. It is therefore justifiable to take the positions 
of the lines correct only to the nearest millimetre (0-012 A) on the records, 
considerably reducing the labour involved in the computations. Portions of 
the approximate band profile were computed in a comparatively short time. 
The portions selected were the band head as far as the 3878 A Felt line, a region 
which appears to be comparatively free from blends, the line P (56) which is also 
well observed, and those R-branch lines having theoretical doublet separations 
of an integral number of millimetres. The computed band, on an arbitrary 
intensity scale, is shown in Fig. 7. 
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To construct a corresponding self-absorbed band, it is necessary to assume 0 
a value for the Planck function in the intensity units of Fig. 7*. From (6.3) it a 
is obvious that if, at the band head, B,,>J,,’, there will be hardly any self- I 


absorption, while considerable self-absorption requires B,,<J,,’, giving By 
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Fic. 7.—Effect of self-absorption on the computed CN 0-0 band profile. (T'=4500°; 74 =2°63). 
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Fic. 8.—Effect of self-absorption on the intensity distribution in the computed CN o-o band. 
(T=4500°; 7 =2°63). 

* Because the intensity gradient in the chromosphere depends on the emitting element, the 
effective height and absolute intensity scale of the Level 1 spectrum are not the same for all elements. 
Hence, when considering the CN emission, it is incorrect to assume that the Planck function is 
approximately equal to the peak intensities of strongly self-absorbed lines such as Hé. 
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on the intensity scale of Fig. 7; this is equivalent to assuming the optical depth 
at the band head to be 7,,=2-63. The corresponding self-absorbed band in 
Fig. 7 was then computed. 

The intensities of corresponding lines in the ‘‘ideal’’ and self-absorbed 
bands were obtained from their areas, and the intensity distributions plotted in 
Fig. 8. A quick glance at Figs. 2, 3 and 5 shows that the computed self-absorption 
produces an effect at least of the right order. While the self-absorption at the 
band head amounts to nearly 50 per cent, the intensity distribution in the self- 
absorbed R branch corresponds very closely to a rotational temperature of 4800. 
This indicates that whatever the true temperature, unless self-absorption at the 
band head is very much greater than 50 per cent (i.e. 7,,>2°6), the rotational 
temperature obtained from the R branch neglecting self-absorption will not be 
very much too high. On this account, agreement of the measured intensity 
distribution at Level 1 with the self-absorbed band computed for an assumed 
temperature of 4500 would not be expected—indeed, it was found to be quite 
impossible to reconcile the two, whatever value of 7,, was adopted. Unless the 
assumption of L.T.E. is seriously incorrect, either the emission from Level 1 
is due to an atmosphere at a uniform higher temperature or there is a negative 
excitation temperature gradient. ‘The lower temperatures already obtained 
for the higher levels show that the latter alternative is probably correct. ‘Thus it 
was not considered worth while undertaking a further computation of the band 
for a higher temperature. 
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Fic. 9.—CN 0-0 band, Levels 2 and 3: comparison of computed self-absorbed intensity distributions 
(T=4500°) with measured intensities. 


Several different values of 7,, were also tried to obtain the best possible 
agreement between the computed self-absorbed intensity distribution and the 
measured distribution at Levels 2 and 3. ‘The best fits, obtained with 7, =3°5 
and 1-8 respectively, are shown in Fig. 9. ‘The good agreement in the case of 
Level 3, when allowance is made for possible blends, shows that the CN o-o 
band on the intensity records is consistent with L.T.E. and a temperature of 4500°. 
The actual value of 7, is not critical; for a range of 7,, from 1°5 to 2°5 there is 
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little variation in the resulting agreement, so small errors in the computed band 
profile will have little effect upon the result. Such errors will arise from the 
inaccuracy of the adopted line positions and the use of the Level 1 profile (the 
true Level 3 profile is wider because of the increase of turbulence with height 
detected by Redman and Suemoto). 

For Level 2 the agreement is satisfactory except in the immediate vicinity 
of the band head. Again the value of 7,, between 3 and 4 is not critical, so the 
discrepancy at the band head cannot be ascribed to experimental errors. The 
CN emission from Level 2 may therefore correspond to a uniform temperature 
somewhat higher than 4500°, or there may be a negative excitation temperature 
gradient. The latter possibility seems the more likely, in view of the temperature 
of 4500° now established for Level 3. 

7. Discussion and conclusions.—Correlating the results for the three levels 
studied it appears that, assuming L.'T.E. and taking self-absorption fully into 
account, a definite excitation temperature can be ascribed only to Level 3. The 
CN spectrum is in agreement with an excitation temperature ~ 4500” at an effective 
height ~400okm. ‘The results obtained from Levels 1 and 2 can be interpreted 
as indicating a negative excitation temperature gradient, the temperature being 
greater than about 4500° at an effective height ~50km, and less than about 
7000° at an effective height ~okm. The value of the temperature gradient 
cannot be determined without detailed consideration of the emission and absorp- 
tion by an optically thick chromosphere with a negative temperature and density 
gradient—a problem which has not yet been attempted. 

The excitation temperatures obtained above may be compared with those 
obtained by the authors mentioned in Section 1. The agreement with Parker’s 
rotational temperature of 4600° + 400° for an estimated effective height < 500 km, 
is very satisfactory. Blackwell’s temperatures of between 5000° and 7500° 
refer to an effective height of 1500km and cannot be directly compared with the 
present results. 

Superficially, Pecker and Athay’s rotational temperature ~4500° at an 
effective height of 1ookm appears to be in good agreement with the present 
results, but detailed examination of their paper reveals otherwise. Pecker and 
Athay’s results are based on an analysis of the region between the o—o band head 
and the 3878 A Fe1 line on slitless spectrograms. The resolution and dispersion 
are greatly inferior to those of Redman’s spectrogram, and the region consequently 
presents an irregular and unresolved profile. The ‘‘ rotational temperature ”’ 
was determined from the peak intensities of the irregularities in the profile, 
the peaks being identified with the lines of the proceeding portion of the P branch. 
This procedure cannot be justified. In the absence of self-absorption, the 
observed band profile depends on the following factors: (1) the relative intensities 
(defining the effective rotational temperature) of the lines of both the proceeding 
and returning portions of the P branch, the latter being just as important as the 
former; (2) the separation of the lines, which decreases towards the band head ; 
(3) the separation of the individual components of the doublets, which increases 
with increasing |m|; (4) the emitted line profile; (5) the instrumental profile. 
With a wide instrumental profile, such as that effective in producing Pecker and 
Athay’s spectra, the combination of these five factors results in a band profile with 
peaks which cannot be identified, either in position or intensity, with the low 
m-number lines. ‘This is clearly demonstrated by Fig. 10, which is identical 
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id with the computed ‘“‘ideal’’ band in Fig. 7, except that the line profile is rectan- 
1e gular and equal in width to that assumed by Pecker and Athay. If their assump- 
1e tions and results are correct, the band profile in Fig. 10 must be almost identical 
ht with that of their eclipse spectrograms. 
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it Fic. 10.—Computed CN 0-0 band profile (T=4500°), smoothed by a wide instrumental profile. 
it 
‘ An attempt was made to determine the rotational temperature from Fig. 10, 

following exactly the procedure adopted by Pecker and Athay: the result was a 
“ negative temperature. The correct conclusion to be drawn from their data is 
s that the region of the o-o band head exhibits an intensity distribution which is 
, very far from being in accord with a rotational temperature of 4500°, when self- 
0 absorption is neglected. It would appear that self-absorption affects the band 
. profile to just such an extent that Pecker and Athay’s wrong analysis, quite 

fortuitously, gives approximately the right answer. 
" The above discussion does not, of itself, invalidate Pecker and Athay’s 


detection of a positive temperature gradient from their spectrograms. However, 
d this was deduced from an increase in the observed general slope of the band 
d profile, an effect which will be produced by the certain decrease of self-absorption 
with height. Their measures do not necessarily indicate a positive temperature 
gradient*. 

> If strict L.T.E. holds throughout the regions considered, as assumed in the 
discussion of self-absorption, then the approximate excitation temperatures for 
the CN emission must be equated to the kinetic temperature at each height. 
There must therefore be a negative kinetic temperature gradient in the chromo- 
sphere, the temperature being defined as ~4500° at a height ~4o0okm, with ‘ 
the value of the gradient uncertain. However, de Jager’s (1957) most recent 
and reliable chromospheric model requires a positive kinetic temperature gradient, 
although the temperature of 4650° at 400km agrees very well with the above 
results. Pagel (1956) has also obtained a positive temperature gradient, his 
temperatures being somewhat higher than de Jager’s, due to the neglect of 
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* Pecker, in a private communication, has subsequently expressed a general agreement with 
these criticisms, but he has suggested that the observed band profile could be caused by the mixing 
of radiation from optically thin layers at different temperatures. However, although temperature 
inhomogeneities may well be present they could not alone explain the band profiles in Redman’s 
spectrogram. 
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turbulence. It is clear that an analysis of the CN emission based upon L.T.E. 
and the general run of de Jager’s kinetic temperatures will lead to a negative 
excitation temperature gradient contradicting the original assumptions. Thus, 
the hypothesis of L.'T.E. is not self-consistent and must be abandoned. 

The opposite extreme to purely thermal emission would be a chromospheric 
spectrum which is simply the scattered photospheric spectrum. At first sight, the 
o-o band head at zero height would appear to support such a hypothesis, as it is 
almost exactly the mirror image of the band head in the Utrecht Atlas (Minnaert, 
Mulders and Houtgast 1940). However, Hunaerts’ (1947) rotational temperature 
for the o-o R branch in the Utrecht Atlas is 4464° + 117°, which does not agree 
with the ~ 7000° deduced originally for zero height. In any case, if all the photo- 
spheric light were scattered by the bottom layers of the chromosphere, there would 
be none left to account for the spectrum at greater heights. 

An intermediate possibility, first suggested privately by Pagel, is that most of 
the CN absorption might occur in the chromosphere. ‘The bottom layers would 
then be excited by dilute radiation from the photosphere with a colour temperature 
of 6700° (Canavaggia and Chalonge 1946), which would determine the distribution 
of molecules amongst the various rotational and vibrational levels of the upper 
electronic state. With the formation of absorption lines in the lowest layers, the 
exciting radiation would be quickly depleted at the vital frequencies, so that 
ultimately the absorption spectrum of the centre of the disk would be controlled 
by the kinetic temperature (assumed to be ~ 4500” at heights above rookm). ‘The 
excitation temperature of the emission from these heights would be the same as the 
kinetic temperature, i.e. ~4500°. The intensity distribution of the emission 
from zero height would correspond to an excitation temperature of 6700’, but 
would be modified by absorption by the cooler overlying layers. Also, if appreci- 
able absorption at the band heads already existed in the high photosphere, the 
rotational levels 20 < K’ < 35 would be under-populated and the emission from 
these levels would be depleted. ‘The final result would be an apparent excitation 
temperature ~ 7000° with depressed band heads, as observed. 

An insurmountable difficulty, which prevents the adoption of the above 
explanation, is that the tangential optical depths in the chromosphere, being of the 
order 1 for strong lines, are altogether too small to account for the radial optical 
depths in the photospheric spectrum, which must also be of the order1. ‘The bulk 
of the CN absorption therefore takes place in the photosphere. A complete 
evaluation of the CN emission and absorption should start in the low photosphere, 
before any CN molecules have been formed, and progress upwards with the 
kinetic temperature, density, pressure, and partial pressure of CN all varying: 
excitation both by radiation and collisions must be considered. Such an analysis 
cannot be attempted at present. 

The interpretations attempted or suggested above all assume implicitly that 
the recorded CN emission can be attributed to a certain effective height or range of 
heights in the chromosphere. It has already been pointed out (Section 6) that the 
effective height of any level depends on the intensity gradient of the line con- 
sidered. The variation of self-absorption with ‘‘ideal’’ line intensity and with 
chromospheric height produces a steeper intensity gradient for weak lines than 
for strong lines of the same element. Hence, the CN band heads must refer toa 
greater effective height than the weak lines, and the absolute intensity scales for 

each must also be different. A comparison of the intensities of CN lines therefore 
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has little meaning, except as a rough approximation. It would be difficult to 
assess the magnitude of the effects of these considerations upon the present work. 
However, it is clear that the result must be a depression of the intensities of the 
band heads and the strong lines relative to those of the weak lines at the same level. 
Thus the chromospheric CN emission may not be inconsistent with de Jager’s 
kinetic temperatures if departures from L.T.E. and the imperfection of the 
observational data are taken into account. 

If the experimental and analytical difficulties could be overcome, the CN bands 
would yield a wealth of information about the transition region between the 
photosphere and the chromosphere. Ideally, one would like to have spectrograms 
of the mottles, and the regions between them, from the centre of the disk to the 
limb—also high-resolution flash spectra with absolute intensity calibrations, taken 
in the best possible observing conditions (from balloons or rockets, perhaps). It 
is unlikely that such data will be available in the near future. In the meantime, 
it is doubtful if much further progress can be made with the interpretation of 
solar CN spectra. 
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Summary 


The radiative capture-cascade equations for hydrogen are solved, taking 
explicit account of the orbital angular momentum degeneracy of the atomic 
levels. All (n,1) levels with n<13 are considered, together with the 
continuum. 

The method of steepest descents is used to obtain asymptotic expansions 
of hydrogen dipole matrix elements involving the quantum number / 
explicitly, and other methods are also employed to obtain useful approximate 
expressions. Hydrogen photo-ionization cross-sections are tabulated for 
all (n,/) levels with n<13, and recombination coefficients are given for 
electron temperatures T,= 104 and 2 x 104 °K, 

The cascade equations are solved for both case A and case B and T,= 104 
and 2x 104°K. ‘The separate / levels for a given n are very far from being 
populated in proportion to their statistical weights. The present calculations 
give, for both case A and case B, Balmer decrements similar to those obtained 
by Baker and Menzel for case B. 

The observed Balmer decrements of six planetary nebulae, when corrected 
for reddening, are in good agreement with the calculated decrements. 





1. Introduction.—Considerable effort has already been devoted to the 
calculation of the statistical equilibrium produced in an assembly of hydrogen 
atoms by processes of radiative capture and cascade (see Baker and Menzel (1)). 
The calculated relative line intensities produced by such an assembly may be 
compared with relative intensities observed in gaseous nebulae. Even when 
the observed intensities are corrected for space reddening, calculations and 
observations may not agree, due to the effect of neglected physical processes, 
such as collisions in hot dense nebulae and line absorption in optically thick 
nebulae. In attempting to assess the possible importance of such processes it 
is important to consider the accuracy of the capture-cascade calculations. 
Previous work has not taken proper account of the orbital angular momentum 
degeneracy of the hydrogen atom. Let N,, be the number density of hydrogen 
atoms in quantum states n/ and put N,, = >) N,,._ Previous workers have implicitly 
assumed JN, to be proportional to the statistical weight (2/+1), giving 


2l+1 
Nu = ne Nu (1 ) 





This relation will be satisfied in thermodynamic equilibrium but may not be 
satisfied when there are appreciable departures from thermodynamic equilibrium. 
The main purpose of the present work is to make calculations in which (1) is 
not assumed. Some of the results obtained will be used in further calculations 
for non-hydrogenic atoms and ions. 
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2. The equilibrium equations.—Consider, first, case A of Baker and Menzel. 
It is assumed that an optically thin nebula surrounds a central star. Ground 
state hydrogen atoms are ionized by stellar radiation. Due to the effectiveness 
of Coulomb collisions, the free electrons will have a Maxwellian velocity 
distribution corresponding to a kinetic temperature 7, (see Bohm and Aller (2)). 
It is assumed that the excited states of hydrogen are populated by radiative 
capture and by cascade from higher levels. One neglects the possibility of 
excited states being populated by absorption of stellar Lyman line radiation 
or by absorption of radiation produced in the nebula. The excited states are 
depopulated by cascade to lower states, the radiation field being considered 
sufficiently dilute to enable one to neglect absorption of radiation by atoms in 
excited states. On equating the number of atoms entering any excited state to 
the number leaving, one obtains for each level with n>2 

aN? + > Ny Ay; = Nj; 2 As (2) 
J>3j 7s 
where j stands for the quantum numbers (n,/), j’<j implying n’<n. The 
number of atoms in state 7 is denoted by N,, the Einstein coefficient for the 
j'~j transition by A;; and the recombination coefficient for level j by «;. The 
electron density N, is assumed equal to the proton density. 

Consider, next, case B of Baker and Menzel, which applies to nebulae which 
are very opaque to Lyman line radiation. One takes account of the nebular 
radiation field by supposing that a transition from level 7 to level 1 with emission 
of a Lyman quantum is counterbalanced by the reverse transition 1->n as the 
Lyman quantum is reabsorbed. Thus for case B, in the summation on the right 
of equation (2), n” takes the values 2, 3, (n—1), and (2) is valid only for n> 2. 
One assumes negligible flux of stellar Lyman line radiation. 

In case A the population of the ground state, N,, is determined by balancing 
the ionization rate against the total rate of capture on to all levels. The population 
of the 2s state is determined on equating the number entering 2s by capture and 
cascade to the number leaving 2s by 2s—1s transitions with two-quantum 
emission (3)* plus the number of collisionally induced 2s-—2p transitions 
followed by Ly, emission (4). j 

If it were assumed for case B that all Lyman radiation produced in the nebula, 
both line and continuous, was subsequently re-absorbed, then N, would be 
determined by equating the rate of ionization from the ground state by absorption 
of stellar radiation, to the rate of capture onto all levels withn>2. These extreme 
assumptions cannot, héwever, be correct since no mechanism for depopulating 2p 
is included. If only radiative processes were considered, the 2p population would 
become so large that absorption of Balmer quanta would be important. In 
practice depopulation of 2p results from the fact that some Ly, quanta do in fact 
escape and from the fact that 2p—2s collisional transitions may be followed by 
2s—Is transitions with two-quantum emission. Whether or not case B conditions 
are ever closely approached in practice can be decided only by solving the Ly, 
transfer problem and by allowing for collisional effects. 

Assuming that in equation (2) all the «; and A;; are known, we see that an 
exact solution of the problem would involve the solution of an infinite set of 
simultaneous equations. Baker and Menzel (1), while ignoring the / degeneracy 
of the levels, solve for an infinite number of levels. In this paper the effect of 


* See also end of Section 6. 
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allowing for the separate / states is investigated by taking into account all (n, /) 
states with m<12 plus the continuum. We solve for both case A and case B 
for T7,,=10000 °K and 20000 “K. 

3. Computation of recombination coefficients and transition probabilities.—The 
Einstein coefficient for the spontaneous transition j->j’ may be put into the form 





3 (87aa,"\ max (J, /’) 
A. =47°v3 : : nl, n'l') (2, 3)) 
iv = 4 ( 302 2I+1 lp( I (3) 
where 
2ne2 
C= = fine structure constant, 


ch 
h? : 
a= —, = Bohr radius 
me 
and where 


p(ni,n'l') = | * P(nl,r)r P(n'l’, r) dr (4) 
0 


is the dipole moment matrix element (in atomic units), P(nl,r) and P(n'l’,r) 
being the normalized wave functions of the (n,/) and (n,/') states respectively. 
Green, Rush and Chandler (6) have compiled extensive tables of |p|?. ‘These 
cover all the bound—bound transitions needed in this work. 
‘The recombination coefficient is given by 
li:w Te or, {* 7 
a; = (2) (mkT,)~9?2(21 + 1 )eli*Te | a;(v)(hv)?e-*"*%e d(hv), (5): 
c\n Jt, 
J 
where I;=hvy= energy of level } hv = energy of emitted photon, and a,(v) is. 
the photo-ionization cross-section which, for the hydrogen atom, may be put 
into the form 


ahiin (=) ay" i(F Jars [Io(ni, 2 — 1) 2 + (1+ 1)lo(nl, R141) ], (6) 





2n" 
where 

o(nl, R21+ 1) = | P(nl,r) 1 F(k?l + 1,7) dr (7) 

0 

is the dipole moment matrix element between the bound state (, /) and the free 
electron state (k?,/+1), k? being the energy of the ejected electron (in units 
of 13-60eV in all subsequent formulae). The normalization of the free wave 
function F(k#/+1,7r) is such that it has asymptotic amplitude k~"?. 

Thus in order to determine the «,, we require a knowledge of the matrix 
elements o(nl,k?/+1) over a sufficient range of k? to complete the integration 
in equation (5). 

Exact analytic expressions for o(nl, k?/+1) have been given by Gordon (6), 
but since these are highly complicated, except in certain special cases, and since 
a large number of such matrix elements is required, simpler approximate 
expressions were sought. ‘To this end, formulae obtained in connection with 
other work on the derivation of an approximate general formula for the photo- 
ionization cross-sections of atoms other than hydrogen (Burgess and Seaton (7)) 
were found very useful. 

At the threshold (k2=0), o and the rate of change of o with k? are given by 
moderately simple expressions. 
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As a first approximation, for small k, it is found that o may be represented by 
y, A(nl,l+1) 
o(nl, k?1+ 1) =0(nl, ol + (2) ‘ (8) 
v 


where f is determined from da/0k?|,,_». Substituting this into (6) gives 


aud») = (=) wei +1) | en ol—1)/ (“) frre 
a | 


Burgess and Seaton (7) obtain 








: 1/2 Pea I’'+2)!"5 —l-1 
L ol’)= ~_ n+l pl’ +2—-2n (Et 
een) ee a ae 
an 
"SDK, 
Pr ee , 
y(nl,l')=2B—1= ime ~ Sager te 5, (11) 


3+) > S Fn, 1) 
where F,(n,/’) is the confluent eeliniiaints function 


FP (U' —n+t—1,2l' +2; 2n) 
and where gy= 1 


U(l+1)—(n—t)(n—t+1) 

2t(nt+l—t+3) atin 

=(l'+n—t+4)( +n—-t+3)q 

— . U(l’+1)—n(l'—n+t—2)] 
ui 
K, t i(m, I’) + ( +.1)(2l +3) Fy. i(n,l’ +1). 

All the functions F;,(n, /’) required were generated by first computing F’, _,(#, /’) 

and F,,_,(n,/') from the series defining the confluent hypergeometric function 


» wees. f a(a+1)2? 
Wie PD="+ Br + BB+) al 


and then using the recurrence relation 
(l’+n—t+3)F,_,(U,n)+2(t—2)F(U,n)—(U—n+t-—1)F,,,(U,n)=0. 


Table I gives |o(nl,ol')| and y(nl,/’) for all levels with n< 12. 
Substituting (g) into (5) gives 


a) 


1 7n ‘ Y , ’ 
fy( Te) = [oy a2( 2p) "YBa pn) + (L410) Bryan) (13) 





“= 














where 


O%4, denotes o(nl,ol+1), yy, denotes y(nl,/+1), and E,,,(x,) is given by the 
incomplete I'-function 


o 2-y 
Bua(s,)= [" & etd, (14) 
in 
where x, =1,,/RT,. 
Putting 10°t=7',°K and substituting for the constants involved gives 


_ 15°788 
n~ tn? . (15) 








uv 





XUM 
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by TaBLe I 
(8) |a(nl, ol’)|/n? and y(nl, l’) 
nl lV |\oal/n? y nol Ul  |olfn? y nol lI — |oal/n? y 
I © I 2°7139 2°6667 8 oO I 1°3551 2°1039 10 g 8 001553 5°6267 
2 0 1 20777 2°0833 I © 0°8638 2°3333 9 10 03107 5°8167 
I 2 1°5847 2:°1082 
I © 075998 2°3333 (t'a—_ oe (0 U Uc ee oe 
‘9) I 2 2°3991 30833 s anak 1 © 08704  2°3333 
oo ee I 2 1°4306 2°1137 
3.0 1 1°8080 2°0317 3. 2 04546 28292 ei eae “ pa 3 
I © 0°7305 2°3333 3 4 1°8244 274446 : : a Bay aman 
0) I 2 2°1916 2°4444 4 3 0'2989 3°1860 ’ aa 47 pan cinl 
2 I 03267 2°8889 4 5 1°7401 2°8246 : ‘ : ‘ae ae 
2 19602 ° orl 6 = 
3 9 3°4444 : : ne ag alee 4. 3 0°4006 =. 29499 
4 0 I 1°6504 2:°0390 4 § 14-7201 2°4469 
6 5 0°'09027 4:2012 - 
1) I 9 9°7920 2°3333 6 7 10504 748 5 4 0°2852 3°2515 
I 2 2°0009 2°2431 7 6 0°03390 cae 5 6 1°6568 2°7342 
2 I 0°4493 2°7470 7 8 o-sgas eae 6 § O'IQI5 3°6155 
2 3 2°0538 2°7917 6 7 1°4884 3°1329 
3. 2 0°19406 3°3542 9 O FT 1°3140 2°1165 7 6 o1189 4°0546 
3 4 1°5525 3°7917 1 © 0°8678 2°3333 7 8 12217 376582 
: , ; I 2 1°5249 2°1077 8 7 006612 4°5835 
$e 5 eo oe 2 1 06654 2°526 8 9 08868  4:3258 
1 0 08255 23333 54 5207 y 4°325 
: 2 38561 sian 2 3 1°6922 2°1826 9 8 0703103 5°2188 
2 3 05359 26657 3. 2 04879 2°7760 9 10 0°5365 5°1515 
. 3 4 1°7800) 2°3581 10 69g 00106 5°9780 
2 3 1°9949 2°5006 c * 3 
4 3 @nees 3°1933 4 3 0°3386 3:0877 IO II 0°23413 6°1515 
3 4 1°8136 3°1333 ; : pie rire 2 © tt 2006 2°1464 
4 3 0°12090 3°7733 s 6 san ; ” I © 08700 2°3333 
l') 4 5 1°2090 4°1333 2.7 ee I 2 1°3925 2°1178 
: 6 5 0O°1270 3°9520 : 
on 6 © 1 1°4653 2°0735 6 7 1:2653 3°6927 , tom? fy 
: , 2°1413 
1 © O°8451 2°3333 7 6 0°06279 4°5432 7 3 
4°543 3. 2 05566 2°6715 
I 2 1°97447 2°1285 7 8 08476 481 z 
se 3 4 16467 2°2264 
2 1 05783 26129 8 7 0°02284 5:2716 eames 3 -Rq88 
2 3 IQIIS 2°3502 8 9 O4112 5°4815 7 2 = 5 pa 
3 2 0°3599 2°9944 . ¢ Se 
3 4 18714 2°98a1 TO © I 1°2792 2°1277 5 4 0°3126 3°1716 
; avnainn aneiiie 1 0 08698 2°3333 5 6 1:°6622 26208 
F : : au i I 2 1°4743 2°1089 6 5§ 0'2193 3°4966 
5 4 0:07764 4°1667 .: oe oo ee Bhs pil sad 
2) - 6 0°9317  4°4722 2 3 1°6351 2°1625 7 6 014484 3°8851 
. : 3 2 O°5152 2°7339 7 8 173308  3°3936 
7 © I 1°4043 2:°0897 3 4 1°7339 = 2°2986 8 7 0°08824 4°3482 
1 0 0°8568 2°3333 4 3 0O°3721 3°0110 8 9g 1°0484 3°9533 
3) 1 2 1°6564 2°1136 4 5 1°7413 2°5339 9 8 004802 4°8998 
2 1 06160 2°5758 5 4 0'2540 3°3497 9 10 0°7293 4°6464 
1e 2 3 18303 2°2662 5 6 176334 2°8851 10 Q O'02199 5°5542 
3 2 O'4129 2°8990 6 5 016077 3°7634 IO II 0°4222 5°4861 
3 4 1°8604 2°5751 6 7 1°4022 3°3702 II 10 0:'007337 6°3264 
1) 4 3 @2gr2 3-3F7O 7 6 o:'09140 4:2678 Il 12 017608 6°4861 
4 5 1°6862 3:0776 7 8 1:0664 4:0076 
5 4 O'13112 3°8544 8 7 o'04401 4°8821 
5 6 1°2850 3°8095 8 9 06771 4°8167 
6 5§ O'05091 4°5442 
5) 6 7 07128 4°8095 


XUM 
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The computation of the incomplete ['-functions for the two electron temperatures 
considered (t=1, 2) was tedious; tables given by K. Pearson (8) were found 
to be useful for some ranges of y, but resort also had to be made to series of 
the form 


oe a oe : ey | 
[i — ie (p—1)(P—2) | (P—1)(p—2)(P—3) | 


We consider now the accuracy of the «,, resulting from the use of the 
approximate expression (g) for the variation of the photo-ionization cross-section 
with energy. Several checks may be made. 

First, the sum over / of the recombination coefficient «,, for a given n should 
correspond to the recombination coefficient given by Baker and Menzel (1). 
On simplification this requirement gives 

16 (Rch\3 g g 
OY fT.) = — (—) & [-E(—x,) em = 36360 & [—E(—x,)]en. (16 
Shu e)= "7 (Fs) S[- Bi(—s,)]erm= 36360 & [-8(—x,)]e». (16) 
This was checked, and it was found that for n less than about 5, the agreement 
was satisfactory, but as increases }/f,,,(7) tends to become too large. 





Tase IL 
ef (t) (t=10-‘T,) 


OL nd) BPE (t) OF ft) 
nm f m 8 n 

t=! t=2 t=1 t=2 t=I1 t=2 

1 © 1757 1183 7 5 916 3°30 10 7 1°79 0°586 
2 0 261 178°3 6 2°23 0°854 8 0560 0°183 
1 600 362 8 o 611 3°92 Q O'III 0°0406 

3 0 87°3 59°! I 15°3 8-70 II Oo 2°49 1°56 

I 227 "37°79 2 22°3 4 I 5°95 3°26 

2 193 101°4 2 ag3 10'9 2 885 4°60 

40 40°5 27°! 4 176 7°16 4 10°2 5°00 

1 107°6 64°8 5 957 3°47 4 942 418 

2 Ira! 63:2 6 3°62 5°23 5 6-91 2°74 

3 61°9 28°6 7 0'787 0:2098 6 4°09 1°44 

5 Oo 23°2 4°7. ? 9 O 4°39 2°79 7 1°97 0°637 
I 58°8 34°9 I 10°84 6:07 $8 0767 0°240 
2 74°5 38°9 2 16:0 8°37 9 0°226 0°0710 
3 55°0 24°9 q 174 8°30 10 0°0435 o'0158 

4 19°9 8°40 4 143 6:07 IZ OO 1°94 1°22 

6 © 13°54 8-86 5 8-88 3°34 I 4°45 2°49 

I 35°2 20°7 6 418 1°43 2 6°78 3°52 

2 47° 24°9 7 1°43 0°473 3 805 3°91 

3 41°8 18°38 8 0289 0°107 4 7°63 3°43 

4 22°8 9°03 10 © 3°32 2°07 5 5°99 2°42 

5 6°52 2°62 I 7°90 4°39 6 3°80 1°38 

7 © 8°85 5°93 2 113 6°13 7 2°04 0666 
I 22°6 13'06 3 13°4 6-40 8 o'919 0°280 

2 32°1 16°8 4 11°65 5°05 Q 0°339 oOrror 
3 31°4 14°2 5 7°94 3°05 10 ©60'0936  0'0289 
4 20°8 8-22 6 4°27 1°47 II 0°0175 0°00632 











1c 
mn 


it 





XUM 
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Secondly, for / = — 1 the exact expressions for o(nl, k®/ + 1) become particularly 
simple (see Appendix), and numerical integration in equation (5) is straight- 
forward. ‘This was carried out for all the (m,n—1) levels considered, and for 
both electron temperatures. 

Thirdly, for / small (+0, 1,2), much more accurate expressions for 
o(nl, k?1+1) are given by Burgess and Seaton (7). ‘hese are of the form 


~ o(nl, kU’) = G(nl,U') [1 + 22k] 7") cos a E + x(n, 1’) 


ta(lst) (2S :) +B r) ) (saa) | (17) 


where G(nl,l'), y'(nl,l'), x(nl,l'), «(1,U') and B(i,l’) are tabulated. In addition, 
the method of steepest descents was applied to the exact expressions for 
a(nl, kl’) in order to obtain asymptotic expansions valid for large n, small k 
and small / (see Appendix, equations (A 18) and (A 19). 

‘These formulae were used to obtain «; accurately by numerical integration 
for 1=o, 1, 2 and for selected values of m. ‘Thus by comparing these results 
with those obtained from equation (12), the errors in the f,,,(7\,) were known, 
for a given n, for /=o0, 1,2 andn—1. For /lying between these values, the errors 
in the f,,(7,,) were estimated by interpolation, at the same time ensuring that 
equation (16) is satisfied. It should be noted that the greatest contribution to 
the summation on the left-hand side of equation (16) comes from the region 
of intermediate /, over which we have to interpolate, hence the condition (16) 
is sensitive to any error in the interpolation. 

‘Table I] gives the values of t®f,,(7',) after correcting in this way. For /<2 
they should be correct to within 1 per cent and for other cases the error should 
never exceed about 5 per cent. 

4. Solution of cascade equations.—_Equations (2) may be put into a form more 
convenient for solution by relating the populations of the levels to those which 
would hold if the system were in thermodynamic equilibrium. We therefore 


put 





N, =N,2(2l-+1) [a] * enb(n, 1), (18) 
where, for thermodynamic equilibrium, b(n, b = 1, and equation (18) reduces 
to the Saha equation for the system. 

Substituting equations (3), (12), and (18) into (2) we obtain, after simplifying, 


T\3 
1 (MLD gr) + & als Neb Kyy= Al Nem) Y Kips (19) 
1 jfi>¥ i” 


where 
- _ max(/l')) 2 1\3 
rae “(al+1) bl (= ~ me) 
Defining ¢; = (2/+ 1)e4b(j) we have finally 
08089 x 10-4 £f,(T,) + a» by Ky, =$; 2 Ky (20) 
yj jy<j 


Equations (20) were solved for ¢,, starting with the highest levels (m= 12), 
and working down through the levels one by one. ‘The values of (7) obtained 


are given in Table III. 
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Tase III 
(a) Case A, T,=104 °K 
A(n, 1) 
mNI ° I 2 3 4 5 6 7 
3 0°00357 
3 0°594 0°0248 0:0618 
4 0802 0'0487 0°123 0°156 
5 0949 ~=— 00665 0168 0246 o°189 
6 1°045 0:0780 = 0198 0°305 0°295 0158 
| I°rl 00855 0218 0348 0367 0262 8 0°107 
8 1°16 00899 0°230 0°375 0°418 0°334 o'189 0:0628 


(b) Case A, T,=2 x 104 


b(n, Ll) 
ml ° I 2 3 4 5 6 7 
2 0'0394 
3 2°68 0'0976 0°194 
4 2°46 0°133 0°279 0°299 
5 2°42 O'Is51 0°327. 0407 80°28: 
6 2°40 o'161 0°356 0°469 0°396 0°203 
q 2°38 0164 0°375 O°517 0°473 0°304 0°'125 
8 2°38 0163 0°386 = 055600539 (0381 0°200 = 0'0703 
(c) Case B, T,= 104 (d) Case B, T,=2 x 104 
b(n, 1) b(n, 1) 

n\l o I 2 l>3 n\l oo I 2 1>3 

3 0°947 0°215 0°0624 4°09 0°846 0°197 

4 1°06 0°307 07125 3°14 0°837 0°283 

5 "13, 0°367 07170 «=Asin 2°82 0°836 0°332 As in 


a 


) "18 o'402 0-201 Case A 
‘21. 0°422 0°221 


‘23. -0°432 «0°233 


2°66 0-829 0°361 Case A 
2°56 0-809 0°380 
2°49 0°785 0389 


eS Aun SW 


— 


7 
8 


5. Results and discussion.—The first point of interest in the results given 
in lable 11] is that they show that the separate / levels for a given n are far from 
being populated according to their statistical weights, as is assumed in Baker 
and Menzel’s treatment; if this were true then the b(n, /) would be independent 
of /. In particular we note that for case A the p levels are grossly under-populated 
relative to the other levels, this depopulation being due to the fact that the 
probability of the transition np->1s is very large. 

The question arises as to how much the 6(n, /) are in error due to the neglect 
of all bound levels above n=12. Let b,(j) be calculated neglecting levels with 
n>N. Results for b,.(m), calculated neglecting /-degeneracy, are compared with 
the b,,(n) of Baker and Menzel in Table IV. Values of 5,,(m) calculated by 
Baker and Menzel for case B, t=2, are also included. 

For case A it is seen that, as would be expected, b,.(n)<b,,(m) for all n and 
that the effect on b(n) due to neglecting the higher levels decreases as n becomes 
smaller. For case B however this is not the case; we see that b,.<,, for all n, 
but that for n<6, b,,(n) is less than both 6,.(m) and b,,(n). This cannot be 
correct, both on physical grounds and from the form of the equations actually 
used by Baker and Menzel. Thus in equation (18) of Menzel and Baker (9), 
all the terms in the summations are positive, and neglecting all levels above say 











(a) 


(b 
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TasLe IV 
Case A Case B 
n T,=10' T,=2 x 104 T,= 104 T,=2x 104 
b(n) b(n) by2(n) b(n) b,() 5,(n) bil) by) 5,(m) 

2 0°00378 0°0456 
3 "0362 + 0°0393 "137. 0°146 0°108 0:°089 O°40I 0°420 0°330 
4 ‘0869 ‘0961 "215 233 *183 +166 ‘448 °466 +404 
5 136 "150 "272 296 ‘245 °233 *487 ‘510 ‘460 
6 176 "202 312 350 ‘290 «= ‘296 “512 «538 «°512 
7 ‘207 "242 339 «= -°391 "321 °34I1 525 °558 °550 
8 "231 ‘277 358 9-422 342 «379 530 «6«°§68 «+577 
9 248 “341 3690-448 353 417 525 °570 ‘600 
10 "258 "340 372 © *470 356 6-448 513 °565 “620 
II 261 *359 366 +488 "448 Cli *489 ; 
12 0°250 0°376 0°344 0°503 O°f4E uss 0°443 

TABLE V 

(a) Case A, T,=104 °K 
b’(n, Ll) 
unl ° I 2 3 4 5 6 | 8 9 
2 0'00386 
3 0°647 0°0270 0'0674 
4 0°882 0°:0536 0°135 0172 
5 1°06 ©0°0745 0188 0°:276 o'212 
6 1'1g 0°0889 0'226 §=0°348 0°336 0180 
7 1°30 O°100 0°225 0°407 0°429 0°307 0°125 
1°39 O'108 0'276 0°450 0°502 O°401 0°'227 0°'0754 
9 1°48 O'115 0'296 §60°487 0°563 0°480 0°309 0°145 0°0415 
10 1°60 0-122 0°317  0°529 0°627 0°560 0°387 0°209 0°0842 0:'0216 
(b) Case A, T,=2 x 104 
b(n, Ll) 

nwt ° I 2 3 4 5 6 + j 5 9 


0°0414 
2°84 0103 0°206 


Ss by 


4 2°66 0144 O'301 0°323 

5 2°66 0166 0°360 0°448 0-309 

6 2°69 0180 0°399 0°525 0°444 0°227 

7 2°74 O'189 0°431 0°595 0°544 07350 0°144 

8 2°81 O'192 0°455 0°656 0:636 0°450 0°236 0:°0830 

9 2°88 o198  0°472 0°706 0°726 0°545 0°320 0°144 0°0445 

10 2°99 0'202 0°486 0°747 o'801 0:°633 0°392 0°199 0'0799 0°0227 
(c) Case B, T,=104 (d) Case B, T,=2 x 104 

b’(n, Ll) b’(n, Ll) 
n\l o I 2 1323 n\Nl o I 2 1>3 


4°34 0°897 0-209 
3°39 0°904 0°306 
3°10 0°920 0°365 
2°98 0°928 0-404 Asin 
2°94 0°930 0°437 Case A 
2°94 0°926 0°459 
2°98 0-933 0-476 
3°04 0°939 07489 


3. 1°03 0°234 0°0680 
4 1°16 0°338 0°137 
5 1°27 O'411 Ov1g! 
6 1°34 0°458 0-229 Asin 
7 1°41 0°494 0°:259 Case A 
8 1°47 0°518 0280 
Q 1°54 0°541 0°298 
10 =: 1'63,:«0'568 -0°319 


0Oo DN AN SPW 


_ 





XUM 
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n= 12 simply involves cutting off all summations at that point. Hence it would 
appear that for case B, the 6,(m) calculated by Baker and Menzel must be in 
error. ‘The method by which they extrapolate from a finite to an infinite number 
of levels seems to be the most probable source of error. 

A rather rough correction for the omitted higher levels can be made to the 
b(n, 1) by multiplying by the ratio C,,=6,,(m)/b,.(n). Table V gives the values 
of b’(n,l)=C,,b(n,1) obtained. ‘These figures are probably not very accurate 
for n and / large, especially since for case B the b,,(m) are in error, so that it was 
necessary to assume the ratio C,, to be the same as in case A. However, for 
moderately small m, the errors should not be more than about 5 per cent. Such 
accuracy should be sufficient for most astrophysical applications. 

An exact solution of the problem taking into account all excited levels would 
be very difficult, not merely because an infinite number of levels must be 
considered, but also because of the fact that, as m increases, at some stage 
transitions due to collisions must become important. ‘That this is so may be 
seen from the fact that we must have continuity between the highly excited bound 
states and the continuum of positive energy states. Thermodynamic equilibrium 
prevails for the continuum states because energy transfer due to collisions is much 
more important than transfer due to radiative processes (Bohm and Aller (2)). 
Thus as m increases we must reach some stage where b(n,/)=1, due to the 
collisional redistribution of excited atoms between states of different m and of 
different /. It should be noted that the observations of Aller, Bowen and 
Minkowski (10) on NGC 7027 confirm that such effects probably occur among 
the highly excited states of He'. Menzel and Baker do in fact show that their 
b,->1 as n-> oo, but this is merely a consequence of the fact that they have taken 
b=1 for the continuum and that the population of a given line is influenced by the 
population of the levels above it. Their b, tends to unity too slowly as n tends 
to infinity. 

Consider a nebula which is optically thin (case A) and for which 7, is too 
small for collisional excitation from the ground state to be important. ‘The 
probability of Ly, emission by an atom in the 2p state is 6-25 x 10° sec~! and, 
for an electron density N,=10'cm~%, the probability of a collisional 2p—>2s 
transition is only 1-7 sec"! (4). Probabilities of collisional n/->nl’ transition 
have not been calculated for n> 2, but it may be expected that for small » such 
processes will also have probabilities much smaller than radiative transition 
probabilities. ‘This is the situation assumed in the present calculations. As n 
increases, the first collisional effect may be that n/--nl’ transitions become 
important and that the method of calculation used by Baker and Menzel will 
be more nearly correct. It may be only at still larger values of n that collisional 
transitions between energy levels become important. 

6. Astrophysical applications—-The main astrophysical interest in this 
problem lies in the spectrum emitted by the system. 

The energy emitted per cm* per second in transitions j->j’ is given by 

Ey = N,hvAjy. 


From equations (3) and (18) we have 


2 
87rd, 


a I I r r 
Eyy=N2[2nmkT,,]-%%4n2 (=) (Rhc)! ( 5) (al+ 1)e™"b(n, IK, 


where R is the Rydberg wave number (= 109 737 cm~). 
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uld : N2/1 : 
ht This gives E,,=2-417x 107*8 7a (= - =) $;Kj; ergcm= sec—. 
ber 
Hence the intensity of a line nn’ is given by 
the 
N2 I I 
Jues Tune = 2°417 x 107%8 p2 > (7 i =) bj Kyy. (21) 
rate 
was Of particular interest are the intensities J,. of the Balmer lines. As is 
for customary, the relative intensities of the lines referred to H, as unity were 
ich calculated, and the values of %,, =/,,5/I4. are given in Table VI, where the results 
obtained by Baker and Menzel are also given, for comparison. In order to be 
uld able to obtain the absolute intensities, the absolute intensity J,, of the H, line 
be is given in Table VII. An alternative way of expressing the absolute intensity 
age of the H, line is by giving the effective value of b, for the 4-2 transition, i.e. 
- > (21+ 1)b(41) Agar 
ind oa = 
um . 2, (2/+ 1) Aor 
ich 
»)). The values of this quantity are also given in Table VII together with the 
the values of 6, obtained by Baker and Menzel. It should be noticed that Baker 
of and Menzel’s values are consistently lower. 
nd Also of interest are the relative intensities of the Paschen and Balmer lines 
ng I.,3/In2, Which we denote by (P/B)™. These are given in Table VIII, together 
oie with the Paschen—Balmer ratios obtained if the / degeneracy of the levels is ignored 
on (denoted by (P/B)™). It should be noted that (P/B)™=(A) 3¥,3)/(AnoYna) 
he is independent of electron temperature and of whether case A or case B applies. 
ids This is not the case for (P/B)™, as is seen from the table. 


In comparing calculated and observed Balmer decrements the first step is 
me to correct the observed decrements for interstellar reddening (Berman (11)). 
he It has been shown by Aller (12) that the corrected observed decrements for 
Balmer lines other than H, are in satisfactory agreement with Baker and Menzel’s 
case B calculations. The present calculations give decrements which, for both 
on cases and both temperatures, are similar to the Baker and Menzel case B results 
ch and are therefore in agreement with the observations of Aller (12). In more 
recent spectrophotometric work, (10) and (13), H, intensities have also been 


on 
obtained. Let J,(A) denote the observed intensities and [,,(A) the intensities 


in 
ne corrected for reddening. Then 

ill log I,,(A) =log Ip (A) + f(A), 

al where f(A) is given in Table IX, the values being taken from Fig. 2 of the paper 


by Whitford (14). The units of f(A) are chosen such that f(A)=o for the 
is ( wave-length of the H, line and f(00)=-—1. For any given nebula we define B 
to be the factor by which the observed absolute H, intensity must be multiplied 
in order to correct for interstellar absorption. With the units adopted for f(A) 
we then have c=logB. The observed intensities for six nebulae*, (10, 13), 
have been corrected for reddening by calculating c assuming the ratios 
I,(H,)/I,(H,)=0°50 and I,(H,)/Z,(H,)=0'285 (see Table VI) and adopting 





* We consider all objects for which Hg, Hg, H, and Hg intensities were measured either by 
Aller, Bowen and Minkowski or by Minkowski and Aller, with the exception of NGC 1535, for 
which the correction for atmospheric extinction appears to be uncertain. 
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Tase VI 


Relative Balmer intensities 


Case A 


T,= 10° 
Baker 
FF, and 
Menzel 
2°48 I'QI5 
I I 
o'5oI 0°576 
0°288 0°374 
o'181 0'255 
O°122 0°182 


Case A, T,= 10‘ 
Case A, T,=2 x 104 
Case B, T,=104 
Case B, T,=2 x 104 


(P/B)o 
Case A 

T,= 10° T,=2 x 104 - 
0°343 0°297 °° 

*399 *350 

"414 “369 

"424 “382 

"430 "394 

"432 "403 

"435 "409 


T,=2 x 10° 


Baker 
F. and 

Menzel 
2°36 1°984 
I I 
O'510 0°560 
0296 0°353 
0188 = 0236 
0128) = 0165 


TABLE VII 


Absolute intensity of Hg 


Tas_e VIII 


Paschen—Balmer ratios 


Vol. 118 


Case B 
T,= 104 T,=2 x 104 
Baker Baker 
fF. and . Pe and 
Menzel Menzel 
2°62 2°50 2°53 2°59 
I I I I 
0489 O'S 0°493 0°50 
0°276 O'31 0280 = 0"30 
0°172 ~=0°206 O°174 O°192 
O'lII4, 0°143 O°II7 ~=0°130 
La b, 
by (Baker and 
Menzel) 
0°132 0-0961 
0°312 0'233 
0*200 0166 
0°494 0°404 
Case B (P/BY™ 
= 104 T,=2 x 104 
260 0°224 0'277 
‘311 ‘273 ‘294 
"329 ‘294 "300 
*340 *3°7 *302 
°347 "319 *304 
“351 °327 *304 


354 334 *395 








fo fed feed fee 
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the mean of the two values of c obtained. It should be noted that the observed 
H,/H, ratio is not used in estimating the reddening correction. ‘Table X gives 


TABLE IX 


Corrections for interstellar absorption 

A fA) A fA) A fA) A fA) A fA) 
3200 +0°447 4400 +0°122 5400 —0'128 6600 —0°358 7800 —0'526 
3400 +0°395 4600 +0:°068 5600 —o'171 6800 —0°390 8000 —0'549 
3600 +0°342 4800 +0°'016 5800 —o0'212 7000 —0'420 rr cas 
3800 +0°288 4861 0*000 6000 —0:'251 7200 —0°449 oO — 1000 
4000 +0'233 5000 —0°'035 6200 —0:'288 7400 —0'476 
4200 +0°177 5200 —0:'083 6400 —0°324 7600 —0'502 


TABLE X 
Observed Balmer decrements corrected for reddening 

NGC NGC NGC NGC IC J Mean 

7027 7662 2392 2022 351 320 observed Theory 
Hi. 2°38 = .2"78 Ss 319 1°75 2°94 2°27 2°55 2°36-2°62 
Hg 1*0O 1'00 1°00 1°00 1°00 1*0O 1°00 1‘00 
Hy, 0°49 O50 0°47 0°52 0°52 0°48 0°50 0°50 
H, o'290 «00'28——s«OO'3 0°27. O27 0°90 0°29 0'285 
B 17°4 2°6 1°6 4°9 16 1°8 


(16-6) (2°9) (2:1) = (3°5) (19) ~— (1°65) 


the corrected intensities and the H, extinction corrections B. ‘The mean 
corrected observed H.,/H, ratio is seen to be in satisfactory agreement with theory. 
Also given, in parentheses, are the values of B obtained if we now also assume 
the ratio /,(H,)//,(H,) =2°50 to calculate a third value of c, and the average 
over the three values of c is taken. 

When an extinction correction similar to that of Table X is adopted, it is 
found by Aller and Minkowski (15) that the observed Paschen intensities for 
NGC 7027 are too weak. With the present results for (P/B)” the discrepancy 
would be somewhat increased. ‘This discrepancy may be due to physical 
conditions being different from those assumed in the theory, but, as pointed out 
by Aller and Minkowski, it may be a consequence of the underlying continuum 
having been drawn too high. 

The contribution to the continuous spectrum of a nebula due to 2-quantum 
emission from the 2s state (see Spitzer and Greenstein (3) and Seaton (4)) is 
proportional to the quantity 


No. of atoms arriving on 2s state 
No. of recombinations on all levels with n> 2° 





X 


TABLE XI 


T, X 

Case A 104 0°13 
2x10! O'15 

Case B 104 0°34 
2x 104 0°37 


- 
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Spitzer and Greenstein, using case B, estimate that X must lie between 0-30 
and 0-35 and adopt a mean value of 0-32. The present calculations enable a 
more precise estimate of X to be made. The results are shown in Table XI, 

Acknowledgments.—I would like to express my sincere thanks to Dr M. J. 
Seaton for his constant guidance in this work, and to Dr L. C. Green for his 
kindness in sending the results of his computations prior to publication. The 
present work was carried out during the tenure of a Ministry of Education 


scholarship. 
Appendix 


Asymptotic expansions for dipole moment matrix elements.——Menzel and 
Pekeris (16) have used the method of steepest descents to obtain asymptotic 
expansions of the dipole matrix elements for bound—bound, bound-free and 
free-free transitions in hydrogen, in all cases, no account being taken of the 
I degeneracy of the levels, i.e. all expressions were implicitly averaged over the 
azimuthal quantum numbers of the initial and final states. In this appendix 
we apply the method to the matrix elements of transitions involving the quantum 
number / explicitly. 

We consider first, transitions between bound states, later generalizing to 
bound-free and free—free transitions. 

For the transition (n,/)+>(n',/+1) we have, in atomic units (Gordon (6)) 


—y)n-i-1 7 1 1 P 142 fp! siti 
p(al,n’l+1)= § 1) (n'+1+1)!(n+l)! ( 4nn "ee ") + - 














4(2l+1)! N (n’—l—2)!(n—1—1)!\(n'—n)P? n'+n 
(A1) 
where 
m ili widen. 
yak (I+2 n',l+1 njal+2; *) 
n'—n\? , : _ —4nn' 
(=) (1 n’,1+1 nsal+2; 4). (Az) 


From the relation 


oF (a,b; c; x)= (1 —x)°-4-,F,(c—a,c—b: c; x) 
we have 





n' + n\%n'+n) 
n’ *) * 
where, writing x= + 4nn'/(n' —n)?, 





_* 
T's (=) oF (l+n',l+1+n; 2l+2;x)—.F,(l+2+n',l+1+0; 2l+2; x). 

















(A3) 
Now, for a>1 and 1<c<a we have, for a, 5, c, integers, 
su cy Se ice, cs 
oF; (a, b ’ C3 x) — T'(6)P(a) dx*-i [x (1 x) ] 

= P(o— c+ 1)I(c) I pa-e - 

= 1(b) 2m1 ims (t —x)*(1 aad t)bt+i-e 

_ TP(b-c+1)P(c) 1 (1—2)"-° : 

7 I'(d) 2ni $ (1 —x—z)tg?ti-¢ dz, (A 4) 


the path of integration in the z-plane being a negative loop enclosing z= (1 — x) 
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but not z=0, since b+1>c in the hypergeometric functions required. Since 
a>c the contour may enclose z=1 if necessary. 

Also, since we always have b> 0, the integral in (A 4) vanishes over a circle 
of infinite radius. Thus the above contour may be deformed into a positive 
loop enclosing z=o but not z=1—x. 

Substituting (A 4) into (A3) and putting B=(n'+n)/(n'—n), a=n/n', we 


obtain 
_ (n—1—1)!(2l+1)! (B2-1\ 1 ff (—2)"**(B?—-2?) 
Y= (n+1)! B? ) ani} 2” 2(p—z)yr te dz. (A5) 
Putting this into a form suitable for applying the method of steepest descents 
(see Copson (17) or Erdelyi (18)), 


_ (n—1—1)!(2l+1)! (BP -1\ 1 2!(B? — 2?) 
y= (Fr) ah lee a | 
x ean z—(1/a) In [(1—z)/(8* —z)}} dz. (A 6) 


Writing f(z)=Inz—(1/«) In [(1 —2)/(8?—2)], we have df/dz=o whenz=—f. At 
this point f(— 8) =i7+(1+1/«)InB. 














The paths of steepest ascent and descent starting at z= —f are given by 
putting the imaginary part of f(z) equal to ix. Writing z=pe"’, this gives 
1 — pe’ 
a(60—7)=arg am } (A7) 
i.e. _! " 9 (p?—1)sin@ 
p 5| @ +1}ent+ lea 


+ £{[ (p+ 1)e080+ EOS 4p} (A8) 


tan (am —«@) 

Equation (A 8) defines two loops, symmetrical about the real axis, one starting 
at z= —f and ending at z=1, the other starting at z= —f and ending at ?. 
The negative real axis is also a solution of (A7). These contours are shown 
in Fig. 1. 
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Defining 
1(z)=f(z)—in— (: + *) In B (Ag) 


we see that r(—f)=o, and that 7(z) is real along the paths defined above. 

The required paths of steepest descent are those along which 7(z) is a 
monotonic increasing function. ‘To determine these, we expand 7 as a power 
series in u=(z+), obtaining 


raen-2G) 1G) + Ce) G) +E) G) 
(itzets) (3) tice. (Ar) 


GS) 


hence if 7 is positive argu~ —}47, 47, 7 and if 7 is negative argu~ —37, 0, 37. 
Thus the directions along which 7(z) increases are as shown by the arrows in 
Fig. 1. 

We see that a suitable steepest descent contour is given by B—A; also, this 
path encloses z=o but not z=1—x (= f?), as required. 

Equation (A 6) may now be written 


Was (n— I—1)!(2/+1)!(B?—1)e-™ 
(n+l)! pete ot 


* [Fase (FR) Jerre (Arn) 


Now, inversion of the series (A 10) gives 


2 hal 2 4 
5 no 1B +2 (1-2) p— 3 (aarp (Bao) pe... (A 12) 





When |u| is small, 














2800 


oreo) |-Gaseses) CS) 


-1)+1 sil _ 2(l+1) , (l+1) , 
= GPR se ~ Bri” BE” 


_ (B,+1— fl) ast 
(41) Gay ts |; 


m/3 
prs (225) (iat), 
a2 
we have 


Loaraaae les) | Bereerac ee eel 


+ = * [4+ 51—02(51+6)] (; =) +0(nh. (A 13) 











Hence, noting that 
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Substituting (A 13) into (A 11) gives 
~inn W(3) ER) (n= 1— 1)! (20+ 1)! (B= 1p te (* =)" 


7 (n+l)! (B+1)%1 12n” 


S {r+ 1259 1(4) [4+ 51—a2(5/ + 6)] +0(Z)I. (A 14) 


30 TG) (a8 en 





Y=(-1)He 








Thus finally we have 


nl, n'l+1)|=1-30 BAe Se n\ts , 
lp( I 3 44 (n’—l—2)!(n+l)! n' n23(1 — n?/n’?)>/8 


. {r+o728 [4 + 5/—(n/n’)°(5/+ 6)] +0()}. (A15) 














(1 = n?|n'*)23n2 3 


The expansion is obviously valid when n> 1, n'’>n and /<n, or when n'>1, 
n>n' and l<n'’. 

Since the expansion is valid for the second set of conditions, we may obtain 
from (A15) 


Ln'l—1)|=1 (n+l)! (n'—1)! n'\!-2 . 
|p(nl, n'l— 1)|=1-3044 lot 1)!(n’+l—1)! (=) n®8(1 —n®/n") 
g [s/+ 1 — (n/n’)*(51—1)] ta}, (A 16) 


x <1 —O'172 eericanatn 

{ 7 (1 — n?/n'?)2/3n78 

with the same conditions on n, n’ and / as previously. In Table AI, values 
of |p|? obtained from (A 14) and (A15) are compared with those obtained from 
the exact formulae. 














TasB_e A I 

|p|? Error 
n l n’ t Exact Asymptotic (per cent) 
5 ° 10 I 2'273 2°233 1°8 
10 ° 20 I 3°195 3°181 0°44 
5 I 20 2 O'2101 0'2269 8:0 
15 I 20 2 89°61 88°35 1°4 
10 I 20 ° 1°562 1*530 2°0 
I5 2 20 3 107°5 105°3 2°0 


We now consider the case of bound-free transitions. By a consideration of 
the hydrogen bound state and continuum radial functions it is easy to show that 


1/2 
|o(nl, k?I’)| = | ees | le{ml, (¢/k)I’]|. (A17) 
Hence from (A 14) and (A15) we have that 


I+1 1/2 
I] (1 +7r?k?) nis 
|o(nl, R21 + 1)|= 16348] = 
| u (Pn) (1 + 2h)58 


r=0 


x {1 40-1728 4+ Sh+ wRNSI+ 6)] | (-)} (A 18) 


(1 4. n*k2 )? 32:3 
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and 
l 1/2 
Il (1—1?/n?) 

re 


Il (1 +7°k?) 
r=0 


— [5/+1+mn?k?(5/—1)] 1\\ 
x {1 01728 (PRA +Of- yr (A 19) 
the expansions being valid for n> 1, 1/k> 1 and /<n. The factor (1 —e-2*/*)-12 
in (A 17) has been omitted since it is very nearly unity over the permitted range 
of k. 
In Table A II values of |o|/n? obtained from (A 18) and (A 19) are compared 
with those obtained from exact formulae. 


nil/6 


=_— = 
|o(nl, R21 — 1)| = 1-63.48 G4 








Tasie A II 
| o| /n? Error 

n l k? I’ Exact Asymptotic (per cent) 
4 ° ° I 1°6504 1°653 o'2 

8 ° ° I 1°3551 1°3557 0°04 
12 ° ° I 1°2226 1'2229 0°02 

8 I ° 2 1°5847 1°618 2°1 
12 I ° 2 1°3925 1°406 I°o 

8 I ° ° 0864 0850 1°6 
12 I ° ° 0870 0864 o'7 


Lastly we consider the case of free-free transitions. If the required matrix 
element is 


x(Ay7l), Ro*lg) = I F(R, 1,7) F(hq2lp, 1) dr, 
0 
where the continuum functions F(k,?/,,r) and F(k,2/,,r) are normalized to 


asymptotic amplitude k,-”? and k,~? respectively, then it can be shown 
(cf. bound-free case) that 





i. ¥ 
—l,,—l, 
anf 1 \32 o(; DE )| 
Ix(A171, Rola) = Po (cz (ree aC ~e-t ia . (A 10) 





We thus obtain 


1+1 
iT (1 +72h,2) 
k,?l, k.21+1)|=2-0489} t= 
Ix( 14, )) 409 it (x 47°) kp 8(1 — k2/k,2)>* 
r=0 


‘a {1 ~01728 [4+ 51—(Re/k,)?(51+ 6)] iy? + O(k,)} , (A21) 











(1 —R,?/k,?)?8 
which is valid for k, <1, k,/k,> 1 and /<1/k, or for k, <1, ko/ky> 1 and 1<1/kg. 
In the course of this work, the asymptotic expansions obtained by Menzel 
and Pekeris (16) were checked. There appears to be an error in the third term 
of all their expansions, e.g. in equation (A 16) of their paper the factor 


(: + Ht at +a) 
15 
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should be (1—4«?+.«4). Hence equation (1.38) for example should read 
~ 7, 08728 t+n?/n?) (1 —$n’2/n? + n’4/n*) 
&1= [ (1 mee n'®/n?)2/3(n’ 2/8 . 0496 (1 = n'?/n2)4/3(n’)*3 Ta 
similar modification being necessary in the formulae for g,, and g,,;, the factor 
+ being replaced by —$. As a check, in Table AIII we compare exact 








— 


Taste A III 
n’ n Aexact A,symptotic 
M. and P. Corrected 
5 10 3°4261 x 107? 3°15 x10% 3-4217x 101% 
10 20 3°6672 x 1076 3°556 x1078 3-6660 x 1078 
20 100 1°3087 x 10% 1°3069 x 107° §=1-3087 x 104° 


values of the quantity A=[,F\(—n’, —n+1; 1; x)]}?—[,F\(—n'+1, —n;1; x)}’ 
(see equation (A 16) of Menzel and Pekeris (16)) with values given by using 
Menzel and Pekeris’s expansion and with those using the correct asymptotic 
expansion (cf. table VI of (16)). We see that the corrected expansion gives much 
better results especially when m and n’ are small. 


Department of Physics, 
University College London, 
W.C.1: 
1958 May 26. 
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Summary 


The paper describes the measurement of 17 solar and vacuum-arc wave- 
lengths in the region 3940-4450 A, and brings to 51 the total number of solar 
red-shifts obtained by the method of circular channels. This group of 17 
lines is of interest in that it includes 6 iron lines from the same multiplet: 
a*F-y°F®; a resonance line of iron: 4375°9 A; and a line of relatively 
great strength in the chromosphere, due to Tit: 4443°8A. The red-shifts of 
these lines do not accord with any of the predictions which have at present 
been made for them. In fact, for the 51 observed lines, the general increase 
of solar red-shift with line strength is still the only observable regularity. 





Introduction 


Although the theoretical aspects are as yet uncertain, the major difficulty in 
interpreting the difference between solar and laboratory wave-lengths probably 
still arises from the scarcity of observations with small and known limits of 
uncertainty. In illustration we may point to the Allegheny Observatory—Bureau 
of Standards red-shifts quoted in Paper V of this series (5) for lines of the same 
strength (Section 2.2); these were selected as the most homogeneous available 
but they show a completely anomalous variation with wave-length. While this 
may be a genuine observational feature of red-shifts, it strongly suggests scale 
errors, different for the arc and solar wave-lengths in the A.O.-B.S. measures. 
Burns and Meggers themselves mention the possibility of scale errors in their 
work. We believe that such errors are substantially reduced by the technique 
we are now using for measuring wave-lengths (Paper IV (4)), and that it is 
therefore worthwhile to redetermine wave-lengths and red-shifts for representative 
solar lines over a wide range of line strengths and spectral region. ‘The Oxford 
measures, so far, have included only three lines in wave-length regions shorter 
than 5080A; in this paper we give results for 17 lines in the blue and violet 
between 3944-4450A. ‘The lines include some of special interest from the 
point of view of theoretical interpretation; these are noted in the comments 
following the observational results. 

In Table I are recorded the lines whose solar and arc wave-lengths have 
been measured. The lines fall observationally into two groups, 3944-4143A 
and 4376-4449 A, which are discussed separately in Sections 1 and 2. Various 
relevant data are included in this table. For convenience, the values of 5/A for 
the lines are shown plotted against log (W/A x 10°) in Fig. 1, where 4 is the residual 
red-shift equal to observed red-shift minus relativity shift, W is the equivalent width 
of the line and A the wave-length; both W and A are measured in the same units. 
The numbering is continued from Paper V of this series. The open circles on 
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the figure refer to lines 1-34 (already published), and the filled circles to the 
new observations. Members of the same multiplet are shown joined. 


TABLE I 


Lines measured 


Wave- Ex. pot. 
No. El. Intensity length Multiplet Transition (low) log (W/A x 10°) 
(Rowland revision values) (Reference : 10) (volts) (W taken from 6) 

35 ~=OA‘d 15 3944°0A I 3°P°-4°S 0°00 2°217 

36 Al 20 3961°5 I 3°P°-4°S o-'o!l 2°356 

37 Fe 10 3969°3 43 a®F-y*F° 1°48 2°057 

38 Fe 30 4045°8 43 a®F_-y®F° 1°48 2°473 

39~—ssé Fe 20 4063 °6 43 a®*F_-y°F° 1°55 2°324 

40 Fe 15 4071°7 43 a®F-y°F° 1°60 2°248 

41 Fe 10 4132°'1 43 a®F-y°F° 1°60 1°964 

42 Fe 15 4143°9 43 a®*F_-y°F° 1°55 2°025 

43 Fe 6 4375°95 2 aD-z'F° 0°00 1°520 

44 Fe 15 4383°6 4! a®F-25G?® 1°48 2°388 

45 Ca 5 4434°97 4 4°P°-4°D 1°88 1°614 

46 Ca 4 4435°7 4 45P°_48T) 1°88 1°456 

47 Fe 6 4442°3 68 a’P-x5D® 2°19 1°554 

48 Fe 3 4443°2 350 b®P-x*®D° 2°85 1°386 

49 ~+‘Tit 5 4443°8 19 a’D-z*F° 1'075 1°464 

50 Fe 6 4447°7 68 a’P-x5D° 2°21 1°534 

51 "ea 2 4449°1 160 a®F-y°G° 1°88 I‘114 


1. Strong lines of iron and aluminium in the region 3944-4144A 

1.1. Procedure.—In the region 3944-4144A, eight strong solar lines were 
measured: six of these are members of the iron multiplet, Fer 43, and the other 
two are the interlocking lines 3962 A and 3944A of Alt. These lines are listed 
in the first part of Table I. The observations were made during 1957 August 
and September by S.N. 

The method of circular channels was used in the determination of the solar 
and laboratory wave-lengths and the final experimental procedure followed 
Paper III (3). The solar lines were observed at the centre of the disk over a slit 
height 1/10 of the diameter of the solar image. The spectrograph was employed 
in its three-prism, diaphragmed form giving a dispersion of about.o-63 A/mm 
in this region. 

For the laboratory lines, a vacuum-arc with vertical pole-pieces was used 
and a rotary oil pump maintained the arc pressure at 4m of mercury; the 
pressure was indicated by an absolute aneroid type of gauge. With a pole-piece . 
separation of 10 mm, currents of about 8 amps were used. Only light from the 
central part of the arc was used so as to avoid any possible complications due to 
pole-effect. Mild-steel pole-pieces were inserted for the iron lines, and carbon 
pole-pieces plugged with aluminium alloy for the aluminium lines. The latter 
were thus obtained without any indication of self-reversal. 

Ilford Rapid Process Panchromatic R.40 plates were used for the solar work, 
and I.R.P.P. and H.P.3 plates for the arc. The plates were measured on the 
modified Hilger machine fitted with a rotating stage. Both the machine and the 
method of reduction have been fully described in the earlier papers. 
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1.2. Corrections to the observed wave-lengths.—The solar wave-lengths are 
corrected for the Earth’s diurnal rotation and orbital motion. The arc wave- 
lengths have been reduced to their vacuum-arc values using Babcock’s empirical 
relation between the depression d (cm-") of a spectral term per atmosphere, and 
the excitation potential V (volts) of the term. This relation is d= 1-15V?—1-93V 
for the triplet terms of iron (7), and the pressure shift is given by the difference 
between the depressions for the two terms involved in the transition. The same 
relation was also used for the aluminium lines because no other information 
is available. In any case, the corrections were very small and of the order of 
o-ooo1 A. 

Both solar and arc wave-lengths had to be corrected for the dispersion of air. 
If conditions differ greatly from N.T.P., this correction is important because 
of the large difference in wave-length region ot the lines measured and the primary 
standard (Cd 6438-4696A). Pressures and temperatures were recorded during 
the observations, and the measured wave-lengths were reduced to standard 
atmospheric conditions by using H. D. Babcock’s tables (8). The corrections 
were of the order of 0-0005 A. 

The phase change corrections were obtained from plates showing the line 
4046A of Hg198, the source being a water-cooled mercury—198 isotope lamp. 
Because of the small wave-length range, the phase change corrections for other 
wave-lengths were determined with sufficient accuracy by extrapolation, the 
slope of the phase change—wave-length curve being known from previous work. 

Although fused silica etalon plates of small wedge angle (2}’) were used, 
a correction for the differential deviation of the two cones of rays, standard and 
line, was made to the wave-lengths ((1), Section 2.2). This correction was of 
the order of only 0-0006 A. 

1.3. The observed wave-lengths—The measured wave-lengths in the region 
3944-4144 A are given in Table II. In column 3 the width of the etalon spacers 
are given. Two spacers, 1} mm and 2} mm, were used for the Fe 1 line No. 38 
and the results are in good agreement. For the solar lines the difference for the 
averages of the spacers is 0-oo11 A and for the arc lines o-o001 A. In the case 
of the line 3969 A, which lies in the wing of the H line of Cat+, 1 mm and 13 mm 
spacers were used, the wave-length differences being 0-0009 A, solar, and 0-0035 A, 
arc. The latter difference, especially, in the arc lines is somewhat disappointing 
and not easy to explain, but a 1 mm spacer, which was only brought into use 
for the H and K lines, must strain any interferometric procedure to the limit. 


TABLE II 
Wave-lengths and red-shifts in region 3944-4144 A 
Sun Arc 
Spacer No. of No. of | Red- 
No. El. mm Wave-length Std. dev. plates Wave-length Std. dev. plates shift 
35 Al 14 3944°0136A +0000! 6 3944°0070A +0:'0005 6 6-6mA 
36 Al 14 3961°5320 ©0005 6 3961°5210 00004 6 Ir‘o 
37 Fe 1,1} 3969°2650 0°O015 16 3969°2581 0°0005 12 6°9 
38 Fe 14,24  4045°8185 00005 12 40458100 00002 12 8°5 
39 Fe 14 4063°5980 0°0009 8 4063 °5932 0'0003 8 48 
40 Fe 1 4071°7468 o-ooll 8 4071°7383 0°0007 8 8°5 
41 Fe 14 4132°0688 0°0004 8 4132°0597 0°0003 8 g'I 
42 Fe 1} 4143°8758 00010 8 4143°8667 0°0009 8 Q'I 
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The relativity shift in this region is 8-6 mA, and it can be seen trom Table II 
and Fig. 1 that all six members of the Fe multiplet, irrespective of their strengths, 
show approximately this value, with only the line 38 seriously different. The 
mean red-shift for the two aluminium lines is also closely equal to the relativity 


shift. 
2. Fe, Ca, Ti lines in the region 4375-4450A 


2.1. Observations.—Observations for wave-lengths at the centre of the solar 
disk and in the vacuum-arc for the nine lines 43-51 were made by the method 
of circular channels, before the investigation referred to in Paper III was 
completed. The light patches for the standard and unknown sources were 
coincident on the interferometer plates but not of the same size. The spectra 
were obtained by Miss Adam in 1953 August, using the spectroscope in its 
three-prism form giving a dispersion of o-95 A/mm. For most of the plates 
the spectroscope collimator was diaphragmed to its usual 3 cm strip (11), but 
in order to reduce the arc exposure time for the faint iron lines and the titanium 
lines the diaphragm was removed. ‘To assess any change in the measured 
wave-lengths arising from this change in the optical train, the iron arc was 
observed through both the diaphragmed and the un-diaphragmed spectroscope. 
(The original vacuum-are with horizontal pole-pieces was used for these nine 
lines, with the usual operating conditions; pressures were of the order of 
4 atmosphere. ) 

The cadmium pattern was photographed on the plate in the form of a simple 
cross, as in Paper I, without the ‘side arms’’, which were subsequently added 
following the investigations of Paper III. As is shown there, this leads to possible 
errors in wave-length due to irregularities in the interferometer plates, but the 
red-shifts obtained with the same plate-orientations for the solar and arc 
wave-lengths, are not systematically affected. In our results in Tables III and IV, 
the groups A, B, C and D, correspond to different plate-orientations and the 
red-shifts have always been obtained from photographic plates in the same group. 

2.2. Measurement and reduction.—T he interferometer plates were the originals, 
of crystalline quartz with a wedge angle of 15’. The large angle necessitated a 
correction (as in Paper I) for the shift of the pattern centre between the regions 
6438 and 4400A. The modification of the measuring machine, whereby the 
photographic plates can be measured from very near the true pattern centre, 
has made the application of this correction a much simpler matter than as 
described in Paper I. ' The displacement of the pattern centre due to the wedge 
angle can be calculated as in I, and it suffices to incorporate the component of 
the displacement parallel to the dispersion with the displacement 2x’ (see (2)) 
ot the Cd pattern centre from the centre of rotation of the table. The maximum 
error arising from this simplified procedure was estimated at 0-0002 A in our 
results. 

The reduction of the plates otherwise followed the normal scheme. The 
phase-change corrections for reflection at the silvered plates of the etalon (silver 
coatings of similar density to those in Paper IV) were taken from Table I of 
Paper IV. The solar wave-lengths were corrected for Doppler effects due to 
the Earth’s rotation and radial motion relative to the Sun, while the arc 
wave-lengths were corrected for the small pressure shifts according to Babcock’s 
formulae (7) for iron lines—also applied to calcium and titanium lines. 
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Gale and Adams (9) have found that although pressure shifts for ionized 
titanium lines are generally quite large that of the line 4443-8 A is comparable 
to the shift for iron lines. Assuming that the variation of pressure shifts in the 
Tit line with excitation potential followed that of the Fe lines at low pressure, 
a small value was deduced for the necessary correction. ‘The pressure shift 
corrections made are shown in Table III. 

2.3. Wave-lengths and red-shifts—The observed solar and arc wave-lengths 
for each group are recorded in Table III. As mentioned above, the absence of 
‘‘side arms’”’ on the plates does not allow us to estimate errors in the wave- 
lengths due to irregularities in the etalon plates, and this has, to a certain extent, 
rendered our absolute values of wave-lengths uncertain. The (A solar-A arc) 
values taken from each group should however remain unaffected by these 
variations; these red-shifts are recorded in Table IV. In ‘Table III, the 
asterisked wave-lengths were obtained from plates taken with the spectroscope 
undiaphragmed. It was found from the iron lines that (A diaphragmed—A un- 
diaphragmed)=0-0035 A. This correction, which is presumably due to the 
change in area of the etalon plates traversed by the light, was applied to the 
‘‘undiaphragmed wave-lengths’’, prior to the calculation of red-shifts and 
absolute wave-lengths. The mean absolute values for the solar and arc wave- 
lengths are also recorded in Table III; the standard errors are rather large due 
to the differences in wave-length from group to group. The red-shifts between 
these absolute values are slightly different from those in Table IV but it is felt 
that the latter have more significance since errors due to alterations of the etalon 
plate orientation are eliminated. 


Tasce III 
Observed wave-lengths in the region 4400 A 


Element Fe _ Fe Ca Ca Fe Fe ‘Tit Fe Ti 
wave-length (in A) 4375°— 4383°- 4434°- 4435°- 4442°- 4443°- 4443°- 4447°- 4449°- 
Pressure shift (arc) 0°0002 0°0004 0°0005 0°0005 0:0006 0:0007 0:0004 0°0006 0:0005 


Group 
A Solar (4) ‘9641 6862 -3488 ‘+1999 6°8083 °7273 
Arc (6) 9556 6786 +3395 "1950 7168 
B Solar (3) ‘9659 «6867 «+3461 ‘1992 ‘8070 7242 
Arc (2) "3396 +1939 "7152 
Cc Solar (5) ‘9656-6871 +3505 ‘1997 ‘8087 7261 ‘1472 
Arc (2-6) 9568 6799 °3412 ‘1954 “7181 
*-3376 *-1926 *-7993 *°7143 **1402 
D Solar (10) +9454 : 
Are (5) *9296 
Solar (6) "5581 
Arc (8) "5453 
Solar "9454 °5581 ‘9652 ‘6867 +3488 ‘1996 -8083 +7260 +1472 
Mean +oo10 +0008 +0015 +0013 +0025 +0013 +0015 +0024 +0012 
wave- 
lengths Arc 9296 «*5453. 9559 67889-3305 «1943. “8028 «+7163 +1437 


+oo10 +0011 +0010 +0007 +0012 +0014 +0005 +0017 +0009 
The figures in brackets indicate the number of plates measured. 
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TaBLe IV 
Red-shifts in mA (region 4400 A) 


Wave-length 


(in A)  4375°95 4383°6 4434°96 4435°69 4442°35 4443°20 4443°81 4447°73 4449° 
Group 
A 8:5 7°6 9°3 4°9 10°5 
B 6°5 5°3 90 
Cc 8-8 7°2 9°4 4°0 5°9 8-2 3°5 
D 15°8 12°8 
Mean* 15°8 12°8 8-6 7°5 8-9 4°6 5°9 9°3 3°5 


* The mean red-shift was calculated by weighting according to the number of plates taken. 


2.4. Comments on results.—In this group of nine lines the red-shifts of the 
lines 'Ti+ 4443-8 A and Fe 4375-9A are of special interest; the former both on 
account of the fact that it is an ionized line and also because of its relatively strong 
emission in the chromosphere, and the latter because it is a resonance line. 

The seven lines 45-51 were taken for the most part in the same solar exposure. 
This is of some importance should secular effects be responsible for the large 
scatter amongst observed red-shifts; for these measures of seven lines on the 
same plate eliminate any uncertainty which may arise from the comparison of 
lines measured at different times. The most obvious characteristic of this group 
is the steady increase of red-shift with line strength. Schréter (13) predicts 
a difference between the red-shifts of neutral and ionized lines, since the lines 
are formed at different levels in the solar atmosphere. Also, Plaskett (12) has 
pointed out that the chromospheric emission line must be superposed on the 
photospheric absorption line and so any asymmetry in the former will lead to 
an asymmetry and apparent wave-length shift in the latter. Our present data 
are hardly conclusive, but this study of seven lines 45-51 does not reveal any 
exceptional red-shift of the Ti* line (see Fig. 1, showing 6/A against log (W/A x 10°)). 

The resonance line of iron, 4375-9 A, has an exceptionally large red-shift. 
(The Einstein shift for this region is 0-0093 A.) This isin contrast to the measured 
resonance lines of calcium (4226-7, 6572-8 A), aluminium (3961-5, 3944 A) and 
titanium (5173°8 A), whose red-shifts are more in accord with some line-strength 
relationship. ‘The Utrecht Atlas shows no blending lines likely to disturb the 
iron resonance line, and it is therefore probable that the abnormally large 
red-shift is real. 

Conclusion.—The results given in this paper bring to 51 the total number of red- 
shifts measured by the method of circular channels. ‘The measurements now extend 
over the whole visible spectrum and the line strengths range from lines of 
Rowland No. 1 to Rowland No. 30. This statistical approach establishes only 
one recognizable regularity in the shifts, namely the increase in shift with line 
strength (Paper II). For the 51 results as shown in Fig. 1, the correlation 
coefficient between 8/A and log (W/A x 10°) is 0-458, with a significance level 
better than o-oo1. The y on x regression line for the distribution is shown on 
the figure. 

The present paper, VI, concludes the series of papers on interferometrically- 
determined red-shifts at the centre of the disk. It is proposed now to investigate 
special lines of theoretical interest, without restriction to method and not limiting 
the measurements only to the centre of the disk. 
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THE QUANTUM DEFECT METHOD 


M. #. Seaton 


(Received 1958 May 27) 


Summary 


The Quantum Defect Method uses interpolated or extrapolated quantum 
defects to determine the asymptotic forms of atomic wave functions. The 
method may be used in the calculation of atomic transition probabilities and 
photo-ionization cross sections, in electron-ion collision calculations and also in 
connection with solid state problems. 

The paper gives a summary of previous work on the fundamental ideas 
of the method and presents some new results for positive energy states and for 
the normalization of bound-state wave functions. Some applications of the 
method are discussed. 





Introduction.—When calculating atomic properties one may consider semi- 
empirical methods which make use of known energy-level data. One approach 
is to construct potentials which will reproduce a number of observed levels 
(1, 2, 3) but this has the disadvantage that each individual case has to be treated 
separately and that the calculations may be very laborious. Hylleraas (4) 
and Bates and Damgaard (5) have developed more general methods for 
transition probability calculations. ‘lhe transition integrals are tabulated directly 
as functions of the effective quantum numbers* v,, defined by 


x2 


Ty = —74 (1) 


Vul 
where 7’, is the energy, measured in cm~', required to remove the n/ electron, 
R= 109737 cm~! is the Rydberg wave-number and z is the residual charge 
(z=1 for neutral atoms). 

The quantum defect p,,, defined by pw, =n—v,,, usually varies slowly in a 
spectral series and tends to a definite limit y,, at the spectral head. For 
collisions between electrons and positive ions one requires the phase differences 
between the wave functions for the electron in the ion field and the corresponding 
wave function in a pure Coulomb field. It has been shown previously (6) that 
the zero energy phase is equal to zp,,. This result has been used in various 
electron-ion collision problems (7, 8, 9, 10). 

Similar methods have been developed for the study of solid-state problems 
(Ham (11)); for these problems one requires the asymptotic forms of the wave 
functions for negative energies other than those corresponding to eigenvalues 
of the free atom. ‘The use of interpolated or extrapolated quantum defects will be 
referred toasthe Quantum Defect Method. Inthe present paper the general theory 
of the method is developed by considering a simplified mathematical problem. 
Some of our analysis is a summary of the work of Ham (11) and of some-much 


* For effective quantum numbers we use v in place of the more familiar r*. 
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earlier work by Hartree (12). These authors were concerned with the negative 
energy states. In the present paper new results are given for the positive energy 
states and also for the normalization of bound state wave functions. We discuss 
the use of the method in collision problems and also its practical limitations. 
In a later paper (Burgess and Seaton (13)) the Quantum Defect Method will 
be used for the calculation of photo-ionization cross-sections. 

1. Valence-electron states——We consider an N-electron atom with nuclear 
charge Z and put z=(Z—N+1). An approximate radial function for a 
valence-electron may be obtained on solving an equation of the type 


ad U(l+1) 

dro 
where ry is in atomic units and E and V are in Rydberg units (13-60eV or 
109737 cm~!). The potential V(r) is such that rV(r)=—2Z for r small and 
rV(r)=-— 22 for r large. 

For E<o solutions of (2) which are everywhere bounded and continuous 
exist only for the discrete eigenvalues E,,, of E. The quantum number m may be 
defined by stating that: m may take on the values (/+1), (/+2),...; m=(/+1) 
corresponds to the lowest eigenvalue; £,,,, is greater than E,, for n’>n. An 
equivalent definition is that the number of nodes in P(E,,,/; r), excluding the 
origin and infinity, is (n—/—1). We put E,,= —2?/v,/? and p,,=n—v,,. Then 
/,,;=0 for the pure Coulomb field, V(r) = —223/r. 

For E>o we put E=k?. The solution P(k?,/; r) which is bounded at the 
origin will have asymptotic form 


» } . 
P(k?,1; r) ~ const. x sin| hr — Slr + zim (2kr) + arg P(t — ;) +3(5) | 


Vir +E] P(E,I;r)=0 (2) 








(3) 
The phase 6, is zero for the pure Coulomb field. 
2. Coulomb functions.--We consider the equation 
2 
<-Sye8 = +E] y= fo) (4) 
where A is not necessarily an integer. inte 
E I 
p=sr, a ala (5) 
we obtain 
d= X(A+1), 2 =| 
— -——>—" + - -5 ] y=0. 6 
E Pp a] ° "7 


For «<o we put «=v with v real and positive and for «>o we put x=ty 
with y real and positive. Then y=2/k. 

Entire analytic functions.—Let y(«,A; p) be a solution of (6). Following 
Ham (11) we shall say that y(«,A; p) is an entire analytic function of «= —1/k? 
if, for p infinite, y(«,A; p) is a single-valued analytic function of € in the entire 
e-plane save for the point at infinity. Such a fynction may be represented by 
an expansion 


y(«A; p)= ¥ f(A p) (7) 


p=0 


34* 
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absolutely and uniformly convergent for «<|é|< 2, p<|#|< oo, Conversely, 
y will be an entire analytic function if such an expansion exists. 

For the theory of the Quantum Defect Method we require, for the case of 
A=l=0, 1, 2,..., two linearly independent solutions of (6) both of which are 
entire analytic functions of «. We also require the asymptotic forms of these 
solutions for p large. 

Solutions of the Coulomb equation.—We consider the following solutions of (6): 





alos P)= Tey ezy Maasaal2Ple (8) 
Yo(k,A; p)=i(k, —A—1; p) 

= Fay Ma -a-ral2ole, (9) 
ge as A(x, eae ae Valk, A; P) (10) 
yalesA5 p)=IalK,A5 p) — Gre Alya(1e 25 p), (11) 
YslA; p)=W,, .s19(2p/*)- (12) 


The functions M,,,,. and W,,,,12 are defined by Whittaker and Watson 
(14) and 


A(x, A)= - eat (13) 
G(k,A)= 2 AA ES Ae.) 
Absa) fub(« +A+1) +(x —A)—21n («)} (14) 


where (x) =d [In r(e)yide. 
The definition of y, may be written 
yalee,A5 p) = (2p)? e-**D(A-+ 1 —K, 20-425 2p/k) (15) 
where 
OP +1 ~n hea 2 es ee - 
The series for ® and for exp(—p/«) converge absolutely and uniformly for 
lp/k|< |B] < o. Since P(A+1—«+0)/'(A+1—«) is a polynomial in « of order o 
one may express y, as a series of powers of (1/«). By Kummer’s transformation 
(14), ¥i1(—«,A; p)=94(«,A; p) and therefore only even powers of (1/«) will 
occur. Therefore y, is an entire analytic function of «. In the same way one 
may show that y, is an entire analytic function of «. 
For A41 (l=0, 1, 2, ...), y, and y, provide two linearly independent solutions 
of (6) but on letting A-/ one obtains 





2 PA+1-—«+o)P(2A+2) (72)’. 


Yo(x,!; p)= — A(x, ly (x, 1; p). (17) 
Since 
-” IK, A5 p) = (2p)"*F a, 44(V/(8p)), (18) 
“Tim ya(sA; p)=(2p)"°J -2,-aCv (80), (19) 
lim A(«,A)=1, (20) 


in the limit of large « (17) corresponds to the Bessel function relation 
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J o44,= —J_y_;. The Weber function is defined by 


cos m(2A+1)J 9,4, —-J _sa_1 

sin 7(2A +1) o 

and by the limit of (21) for Al. Then./,,,, and Y,,,, are linearly independent 

for all A. Wannier (15) introduced the function y, by analogy with (21). From 
equations (18) to (21) one obtains 

Tim ya(0,s p)=(2p)!*Va, av (8p): (22) 


Defining y3(«,/; p) as the limit of equation (10) for A>/ one obtains 


yas ls p)= {Ald n(KAs p)+ Z yalesAs p) 





Yoni ia 


tyntets es Als. (23) 


Since y,(x,A;p) and y(x,A;p) are entire analytic functions for all A, dy,/dA 
and dy,/dA are also entire analytic functions. Also, since 
+1)(k+l—1)...(«-l 
Fe ei (24a) 


_ [e®@ =P] fe? (I-12)... [e® = 1] (248) 


Ke 








and A(x,0)=1, A(x, /) is a polynomial in 1/«? of order /. But 0A(x,A)/dA cannot 
be expressed as a convergent power series in 1/«? and therefore y,(x,/; p) is 
not an entire analytic function of «. The entire analytic function y,(«, /; p) is 
obtained (11) on subtracting from y3(«,/; p) the term in 04/0. 

The required entire analytic functions are therefore y,(«,/; p) and y,(«,/; p). 
The expansions of these functions are discussed in the Appendix. 

Properties of the function G(«,1l).—It is convenient to put 


GK, 1) = G(x, l) +130 («, 1) (25) 
where Y and # are real. From (14) G(v,/) is real for v real and positive: 


therefore 
H (v, ly=o0 (v real and positive). (26) 


For y real and positive G(iy,/) is complex: using 


(x) — (1 — x) = — cot (7x) 


one obtains 


H (ty, l)= Pies _ (y real and positive). (27) 
Relations involving the Whittaker function.—It is convenient to express ys 
and y, in terms of y, and the Whittaker function y;. Using the relations 


['(2\ +1) 





W,.asile= T(-a-«) atie + Taticn) 7 (28) 
(14) valid for A#/ and 
C(x)P(1—x)= > = (29) 
one obtains 
- id a oF 
¥5~ Sin (2m) resees Patt =}: (3°) 
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We use (30) to eliminate y, in (10) and then let A>/. This gives 


yale 5 p)= — AK I) cot (mee)yy(4 15 p)— ~FF*=") yle,15p) (31) 
and from (11) 
yale 5 p) = — [A(x 1) cot (mx) + G4, DIva(0s 15s p)— “FAH y(t p). (32) 


Asymptotic expansions.——In the limit of p>, y,(«,/; p) ceases to be a single- 
valued function of x. We therefore consider k=v+?y with v>0, y>o. The 
asymptotic expansion for y, is then (16) 


kitle—itml—«) Gen exp [p/« +137(/—«)] 





yilx,l; p) ~ 





ota T(i+1-«) 
_ (2p/«)* exp [—p/k —137(/—«)] 
P(+1+K) \. (33) 


It is seen that y,(v,/; p)>0o for poo, except when v=n=(/+1), (/+2),...; 
when v=n, I'(/+1-—v) is infinite and y,(n,/; p)-o for p>. 
Substitution of « =7zy in (33) gives 


l: 2y +1 ea7yvl2 
yilty, 4; p) poo P(+1—iy)| sin (x) (34) 
with 
“ — $lr+yln(2p/y) + arg C(/+ 1-7). (35) 
Using (29) one may show that 
, +1 A (iy, 1) 
P(l+1—-#y) P= “4 
IP (+ 1-H p= (36) 
and therefore 
. 2y 1/2 I —e-27y 12. 
’ l; ned — oy ee x}. 7 
valints) ~ (2) "| See | sin) (37) 
The asymptotic form of the Whittaker function is : 
Ys(k,l;p) ~ (2p/ejre’*. (38) 
p> 
Therefore, for all v, y;(v,/; p)>o for poo. From (38) one obtains 
D(il+1-—ty) ,. 5. 2y\"2[ A(iy, 1) ]? _.. 
ty oti e) ~ (2) Laem | - 


and, using (25), (27), (32) and (36), 


valintie) ~ ~ (2) { [ae] Me dsin y+ [ AZ] cos 0} 





(40) 
and therefore 


1/2 ‘ , ! 
Laliyts p)+ Her Dyalints ol ~ ~ (22) | ASD] cose). (4x) 





p> := a 
3. The modified Coulomb field—We consider 
qe ad i ei 
ter p)+2 - =| F (x,1; p)=0. (42) 


Let po, p, be such that 0 < po, py <p, < © and let u(p) be such that: pu(p) analytic 
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for p<po; u(p) piecewise continuous for py<p<p,; u(p)=o for p>p,. Let 
F («,1; p) be a continuous solution of (42) satisfying 


F (x, 1; 0)=0, (43) 
F (1; p)=yiles 1; p)—Ble)yalest;p) (p> pi). (44) 
It is shown by Ham (11) that A(«,/; p) will be an entire analytic function of 
e=—1/k*. Since y, and y, are entire analytic functions it follows that B(e) will 


be an entire analytic function and from the fact that y, and y, are real for «¢ real 
it may be shown to follow that A(e) will be real for « real. Therefore B(e) may be 
represented by an expansion 
Ble)= & Bye? (45) 
p=0 
convergent for |e]<|@|< 0 and the coefficients B,, will be real. 
The negative energy states—From (32) and (44) 


F (1,1; p)={1 + B(e)LA (1) cot (we) + G(s, Piya J; p) 


+ ple) Ett =) 


aK! 


¥s(x,/;p) (p>pr) (46) 
and since G(v,1)= G(v, ), 
F (vil; p)={1 + B(e)[A(v, /) cot (xv) + G(v, 1) ]}}y4(v, 1; p) 
Cd+1- 
+ Ble) wh ak be 
TV 

Let ¢«, = —1/v,? be an eigenvalue of «; then -F(v,,/; p)->o for p>oo. If u(p) 
is not identically zero we may assume v, #n and therefore y,(v,,/; p)> for 


p>o. The condition F(v,,/; p)->o for p> then requires that the coefficient 
of y, should vanish in (47); this requirement is (11) 


ys(v,l;p) (e>pi) (47) 


BlEn)= — [A(rn, !) cot (av,) + F(v,, 1] (48) 
or, in terms of the quantum defect defined by 
- He, )=2— ry, (49) 
it is 

BCE, )=[A(v,, 1) cot (mH(En)) — Grn A). (50) 


Suppose all the eigenvalues ¢, to be known. Then A(e,,) is known for an infinite 
sequence of values of ¢,, and, in consequence of the possibility of the expansion (45), 
B(e) is uniquely determined for all finite «. A quantum defect (e) may therefore 
be defined for all finite « by 

Ble) =[A(x, 1) cot (mu(e))— G(x, 1) (51) 
together with the condition that (e) should be a continuous function and should 
satisfy (49) at the eigenvalues. Since Y is real and B(e) and A(«,/) are real 
for real ¢€ it follows that* (e) will be real for real e. 

The »-defect.—Letting v->m where m=o,1,2,.../ we find that Y(v,/) 
remains finite but that A(v,/) tends to zero. However f(— 1/v*), as determined 
from the eigenvalue spectrum, does not normally tend} to — 1/Y(m, 1). Therefore 
as, v->m, (—1/v?) must tend to an integer in such a way that A(v,/) cot (zp) 

* The definition of u(¢) used by Ham (11) is obtained on replacing G by G in (51). For € real 
and negative this is identical with our definition but for ¢ real and positive u(e) as defined by Ham will 


be complex. 
+ This may be shown (11) by solving the radial equation for various functions u(p). 
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remains finite (rr). This is confirmed by examination of observed energy 
levels. Fig. 1 shows the values of for the nd series in Ki and in Cair; in the 
former case ->o and in the latter p—>1 as v2. 

Ham introduces the y-defect defined by 


A(x, 1) cot 7 = cot 7. (52) 
We consider that (e) is a continuous function and that y(0)=u(0). Since 


does not necessarily tend to an integer as y->m we may expect that 7 will usually 
vary more slowly than y». Values of y are shown in Fig. 1. 
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The nd series of K 1. The nd series of Cat 
Fic. 1.—Quantum defects and n-defects. 


The positive energy states.—Equation (44) may be written 


F =(1+B9)y,—Blyyt+ Fy] (p > p1)- (53) 
Using the asymptotic forms (37) and (41) we obtain 


Flintie) ~ (2)" lean | Ble) {| ar + Hl 1 | sin) 
+| 28 cos(s)b (54) 





1—e~*y 
F (iy,; p)~ (2) C(c) sin [x +8(€)] (55) 
where 
_ [A/(1-e*”)] 
tan$= (CDEKAL (56) 
and 





cacayre] HS Pad | ta] +eaie—etmyph™. 657) 
Substitution of (1/8)+ Y= A cot (zp) from (51) gives 


tan d(e) = eee (58) 
and 
C= (2)!?[A(1 —e-”) 2B { [cot mp]? + [1 — er] -2} 12, (59) 
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To the phase d(e) defined by (55) we may add any integral multiple of 7. 
Consistent with (58) we may define 5(e€) uniquely by requiring that 5(e) be 
continuous and that 

5(0) =zp(0). (60) 


In most applications we may assume 2z7y to be large. Neglecting* exp (—27y) 

in (58) and (59) we obtain 
5(e)= me) (27y>1), (61) 

and C=(2A)!?8/sind (27y> 1). (62) 

The zero-energy phase and the number of bound states.-For an attractive 
potential V(r) which, for large 7, tends to zero more rapidly than the Coulomb 
potential a general theorem (17) states that 6(0)=7N where N is the number 
of bound states obtainable with the potential V. For a potential tending 
asymptotically to the attractive Coulomb form the number of bound states is 
infinite. In such a case one may define N(e’) as the number of bound states 
with ¢,, less than some small negative value <’. Then N(e’) will increase by unity 
when ,(e’) increases by unity. One may thus regard p(o) as the increase in the 
number of bound states which results from the non-Coulomb part of the 
potential u(r); this despite the fact that (0) is generally not an integer. The 
relation 5(0)=7p(o) may therefore be regarded as a generalization of the 
relation 5(0)=7N. 

4. Normalization of the radial functions.—-We put 





F =F («,1;p) 
F' =F («’, I; 
(x p) . (63) 
aa 3% 
dp 
Then from 43) 
gl p W(1+1) F' 
D?- - +u(p += -— F' dp =o. (64) 
Integrating by parts and using the equation aaa by F one obtains 


(FDF' —~F'DF) + ( 7 =) | $F" ene (65) 


Assuming D¥ and D¥’ to be finite for p=o it follows from (43) that 
(F¥ DF' —F'DF ) is zero at the lower limit and therefore 


° ea, 1  (F'DF —-F DF’) 
| FF’ dp= -a (66) 


We consider this expression in the limit of poo. 
Equation (66) is the standard expression used for the normalization of positive 
energy radial functions. One obtains (18, 19) 


{ °F (iyo ls PYF (iy Ii p)dp= (Cle) P8(e-€) (67) 


where 8(¢—e’) is the Dirac 5-function. 








* The condition (27y)>1 is equivalent to k<(2z7z) where k? is the electron kinetic energy in 
units of 13-60 eV. 
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To obtain the normalization for bound states we put 


F =F \(v,l; p) (68) 
and 
F' =F ,=F (vy, 1; p) (69) 
where e, = —1/v,? is an eigenvalue. The integral 
A= | °F 2dp (70) 
0 


is evaluated as * 











«hi lim {7 . ads (71) 
PO vy (1 1/v* i 1/v,”) 
Using (33), (38) and (47) we have 
F ~ ap’ e~*!” + bp-” ef! (72) 
where 
P(/+1—v) /2\’ 
= ak (=) : (73) 
TV v 
v+1(2/v)” 
wat aa) : 
b={1+ B[A cotav+G}]} Titi) (74) 
Since F,, is an eigenfunction, 
Fy, ~ay,prueP''n (75) 
where 
P(l+1 2\"n 
a, = Ble,) =) (2). (76) 
TV»! Vy, 
Substitution of (72) and (75) in (71) gives 
A, iii 24,2 (77) 
ov one 
Using (74) and (51) one obtains 
3b] __ tA (Yay IYB(n)0n!*2(2/04)—"L (rn) (8) 
av |e», sin® (nu(e,))PU+1 —»,) 7 
where 
> 0 
C(v)ar+ 5. (79) 
v 


From (76), (77) and (78) the normalization integral is 


_ yn? B(E, )A(v,, I)C(v,) 
An ain mle) ss 
The normalized radial function P,,,(r) is taken to satisfy 


| ” P,2(r)dr=1 (81) 


where the integration is over r and not p=zr. Making a convenient phase choice 
we put 





g! | es )" a 
P.,,= ‘e ; (82) 
n 





* Ham (11) quotes an expression similar to (71) and states that it can be evaluated using the 
Quantum Defect Method. He does not, however, give the explicit expression for the normalization 
factor as obtained in the present paper. 
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From (47) and (50), 
FF = Ble, )P(L+ 3—v,) 
— 7 





SA yor tsp) (p> rr) (83) 
and therefore 


P.n=K(ry,,)ys(Yns ls p) (p> px) (84) 


K(v,, l) = 32 [C(v,, )v,71'(v, > l +1 P(r, = N))-**. (85) 


If all the eigenvalues are known dy/dv may be calculated and the exact 
normalized eigenfunctions determined for p>p,. It may be noted that 


where 


20 
i(v)=1+ & =1+ 55. (86) 
de 
For large v, ou/de remains finite, and therefore {(v) tends to unity. The 
normalization factor then becomes 


K(v,,/)=2"?[y,70(v, +1+ 10 (v,-)" —— (v, large). (87) 


With v, =n this is the exact normalization factor for the pure Coulomb field 
(for which »=oand{€=1). Hartree (12) has shown that the factor (87) satisfies 
a necessary recurrence relation for values of v, differing by integers. This is 
consistent with the exact result since the v,, may be considered to differ by integers 
only when op/dv is negligible. 

Numerical results for the bound-state normalization factor.—Hartree (12) 
and Bates and Damgaard (5) have shown that the normalization factor with (= 1 
is a good approximation when v,, is not too small. In illustration of cases for 





TABLE I 
The normalization factor K for 2p radial functions of atomic oxygen ions 
K 
Jon State Vo r A _ 
Eqn. (85) with Eqn. (85) with Exact value 
J +2 
{a1 {= (v2 (ve ) 
v(V2+ 1) 
O° 2p’ *P o'890 imaginary o'61 0-78 
O* 2p° *P 1:296 0°38 0°67 0°67 
OF 2a 7 1°483 0°48 o'71 o’71 
oF gp F 1°672 O'sI 0°69 0°65 


Note that the values of v, are obtained from the energy parameters in the radial equations 
and not from measured energies. 


which this approximation is not sufficient we consider the 2p radial functions 
for atomic oxygen ions calculated by Hartree, Hartree and Swirles (20). Fig. 2 
shows the quantum defects for the 2p and 3p states of O**. It is seen that pu 
goes to unity for v=1 and that p(e) is a nearly linear function of e for v between 1 
and v,. Assuming an exactly linear function we obtain 

pip = Sa a3). (88) 

V9(v2 +1) 

Table I gives the normalization factors as calculated with ¢=1 and as calculated 
using (88) and also the exact values determined from numerical evaluation of 
the normalization integral (81). 
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Fic. 2.—Quantum defects for O** np. 


5. Discussion of the Quantum Defect Method.—The mathematical theory 
considered so far is concerned with central fields which tend to the attractive 
Coulomb form at large distances. A different approach would be needed in 
order to develop an exact theory applicable to many-electron atoms. Some 
of the results of the present theory may, however, be used directly. 

Consider an atomic spectral series for which the ion is left in the ground 
state at the series limit. From the observed energy levels quantum defects maz 
be obtained and extrapolated to positive energies. Just beyond the series limit 
the free electron will have a kinetic energy which will be too small for collisions 
to produce excitation of the ion. A single real phase parameter will then suffice 
to describe the asymptotic form of wave functions corresponding to definite 
angular momenta. In such cases there can be little doubt that accurate phases 
can be determined from the extrapolated quantum defects. 

For higher energies the situation is mere complicated. The asymptotic 
form of the wave function may be described in terms of a scattering matrix S, 
the off-diagonal elements of this matrix determining the probability of the free 
electron producing transitions between states of the ion. If all atomic energy 
levels were considered, including those corresponding to excitation of two or 
more electrons, it may be expected that information about the complete scattering 
matrix could in principle be deduced. The more complete theory required has 
not yet been developed. We consider some applications of the present form 
of the Quantum Defect Method. 

Extrapolation procedures.-For the modified Coulomb field we have shown 
that the mixture parameter A(e) can be expressed as a convergent power series 
in « and that the quantum defect may be defined for all € as a function of f. 
We have not shown that the quantum defect itself may be expressed as a 
convergent expansion in powers of « but it may be shown (12) that » may be 
expressed as an asymptotic expansion in powers of «. 
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In practical applications information concerning the eigenvalues is always 
incomplete. Quantum defects obtained trom observed energy levels may 
usually be represented by polynomials in « containing only a few terms. ‘These 
polynomials may be used for extrapolation to positive energies. The range of 
extrapolation is limited by incomplete knowledge of the exact eigenvalues and 
also by the fact that the simple theory is invalid for positive energies such that 
inelastic collisions may occur. 

It may be convenient to extrapolate the y-defect when 7 varies more slowly 
than p. 

Comparison of calculated and empirical phases.—Such comparisons have been 
made in previous papers (6, 7, 8, 9, 28). In favourable cases, phases obtained 
from solutions of the continuous state Hartree-Fock equations (21) are in very 
good agreement with extrapolated quantum defects. The reason is that the 
Hartree-Fock equations for continuous states may be derived from a variational 
principle for the phase just as the equations for bound states may be derived 
from a variational principle for the energy. In cases for which the Hartree-Fock 
phases are less satisfactory it may be expected that configuration interaction is 
important; an example of such a case is the continuum corresponding to 
extrapolation of the 4snp series in Ca (7). 

Calculation of elastic scattering cross-sections.-We consider ions with non- 
degenerate ground states. For elastic scattering of slow electrons the appropriate 
theory is then that for a modified Coulomb field. For scattering through an 
angle @ the cross-section per unit solid angle is (22) 


1(8) = [aof()P (89) 


where a, is the Bohr radius and where 


(0) = 4. S (21+ 1)P\(cos 6)[e%+9 — 1] (go) 
21k 1=0 
and* 
$y) =arg P(/+ 1—ty). (91) 


It is convenient to put 

1(8) = 1,(8)- (8) (92) 
where /,(@) is the cross-section for the unmodified Coulomb field. The latter 
is given by 


T.(9) = laofe(®)P (93) 
where f,(@) is obtained on putting 6,=0 in (go): 
f.(0)= > (21 +1)P(cos )fe%*— 1]. (94) 
1 
This series may be summed (22) to give 
= ba . . 1 
f.(9) al 2k sin2 (0/2) exp [27($o +yIn sin 58)] (95) 
and 
1.0)=| a | (96) 
AY) | Resin? (Oa) | 7 


* We recall that 5,(y) is the phase for the /-wave due to the non-Coulomb part of the potential 
and that y=2/k. 
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which is the Rutherford formula. From (go) and (94) we obtain 





f(0) =f.(0) + = & (al+ x)Prerife%— 1] (97) 
and, using (95), 
iG =1I+ semen p3 (21+ 1)P,(cos 6)e1sin 8, (98) 
with 
p(6,7)=2($1— $0) — 2y InsinO)2 +8. (99) 


It may be noted that 


e2il6y—t0) — T(/+1-ty)[(1 + ty) 
C(/+1+7y)C(1 -ty) 


_ (1-ty)(2-1y)..- (/-ty) 








(1 +ty)(2 +iy)... (1 +4y) tid 
and therefore 
l 
(%:—$0) = — 2. arctan (y/m). (tor) 
From (99) we obtain finally 
I(0)=X2+4 Y? : (102) 
with 
all 
X=1+ — > (2/+ 1)P,(cos 8) cos p; sin 8, 
ai (103) 


. 2 
Y= 2 sin* (0/2) > (21+ 1)P,(cos 6) sin p, sin 8, | 
Y l 


It follows from these expressions that .4(0)=1 for all k* and that, for k?=0, 
JI (#)=1 for all 6. 

Fig. 3 shows results obtained for scattering of electrons by Nat ions. The 
phases* have been obtained using formulae given by Ham (11) for extrapolation 
of the y-defects. The cross-sections are seen to be significantly different from 
the corresponding Coulomb cross-sections. Complicated interference effects 
are particularly marked for small values of k. 

Inelastic collisions——For ions with degenerate ground states the Quantum 
Defect Method may be used to obtain accurate cross-sections for electron-induced 
transitions between ground-state fine structure components. ‘The method has 
been applied (10) to the transitions C+ 2p4/.-2p3o and Sit 3p4/.-3Ps3j2;_ these 
are of importance for the cooling of the interstellar gas. 

For other transitions in ions produced by electron impact it has not yet proved 
to be possible to deduce cross-sections from quantum defect data alone. One 
may, however, make calculations in which the various parameters are so adjusted 
as to be consistent with extrapolated quantum defects (10). 

The Quantum Defect Method may be used to obtain continuum wave 
functions having asymptotic forms which are virtually exact. Such functions 
have been used for the calculations of photo-ionization cross-sections (13). 


* Only s- and p-phases have been taken into account. The higher order phases are small and 
will have little effect on the cross-section. 
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(97) 

(98) 

(99) 

100) 

IOI) 

102) 

103) ° 4 ad 

Fic. 3.—Elastic electron scattering by Nat ions. £ (0) is the factor by which the Coulomb 
cross-section must be multiplied, @ being the scattering angle. k? is the kinetic energy in Rydbergs 

=, (13°60eV). Note that .4(0)=1 for k=o. 

The | APPENDIX 

tion Expansions of the Coulomb functions in powers of the energy.—The expansion 
rom of the regular Coulomb function y, has been discussed in a number of papers 
ects f (see (23, 24, 25, 26, 27)). One may first obtain 

um ' Vi(k, F; p) a = a,(k, 1)(2p P+? 9.1 + »(4/(8p)) (A 1) 

p=0 
“i where (26) 
” @y= 1, a, =0, ay = (1+ 1)/4K?*, 
lese 
ts ; . (A2) 

wed ay= 4p (21+ p)a,_»— a, —3) (p>3) 

yne Putting 

ted a,(«, l) = bs aq, p(l)n*4 (A 3) 

q 
one obtains 

ave 

one Q@gy=0 unless 2q9<p<3q, 

3). | Q,o=1 a, o=(1+1)/ — 

o,o= 1, 12> 4, ee <5 
“a 7 12 +.  (A4) 
I 
aq,p= gp 2! +P) tet, p-2- 40-1, ps) (q21, p23) 
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‘The expansion for y, is 


ee) 3q 
Niles; p)= -_ Gq, p(/)(2p PPT ar41+0(V/(8p))- (A5) 
a= p=24 


Various expansions for irregular Coulomb functions have been proposed. 
It has been shown (see Section 2) that the function y, cannot be represented by 
a convergent expansion in powers of the energy. It may, however, be represented 
by the asymptotic expansion 


Q 3q 
yal, 1; p)= A(x, 1) & «4 2 40, vp NO Varsr+o(V/ (8p) +O), (A6) 

q=0 p=: 
(25, 26,27). On expressing A(«,/) as a polynomial in «~*, (A6) gives an 
asymptotic expansion in powers of x~*. The function G(«x,/) (defined by 
equation (14)) may be represented by an asymptotic expansion in powers of «~? 
(12,27). Substitution of the asymptotic expansions for y, and for Gy, in 
¥4=Y3— Gy, gives a convergent expansion for y, (Ham (27)). For numerical 
calculations the asymptotic expansion (A 6) is often convenient. 


Department of Physics, 
University College, London, 
W.C.1: 
1958 May 26. 
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GIANT STARS OF TYPE II 


C. B. Haselgrove and F. Hoyle 
(Received 1958 March 15) 


Summary 


Integrations along the Type II giant sequence are given that reach an 
evolutionary stage at which there is an appreciable energy yield from nuclear 
reactions in helium. The new integrations seem to be in close agreement 
with observation. 





1. Results.—Seven evolutionary steps beyond those given in a former paper 
(1956 a) have been made for a star of mass 2°512 x 10%° gm., the initial composition 
being the same as before. 

Results are set out in ‘lable I below, the solutions being numbered so as to 
follow consecutively with those obtained previously. ‘The symbols appearing in 
the table have the following interpretation : 


R =total radius, 

M = total mass, 

L = luminosity, 

Tott = effective temperature, 

| =the contribution to the luminosity from nuclear reactions 
that convert helium to carbon, oxygen, and neon. 

P, = central pressure, 

Ty = central temperature, 

Ma. = mass of helium zone, 


M.. © x.= mass of carbon-oxygen-neon core. 
T ¢ and Ty are in deg. K, all other quantities being measured in c.g.s. units. 
The method of solution was that given in an earlier paper (19565). 


TABLE I 

Solu- 

tion log,R logy, L ate He logioT upp logig Po logio To Muye/M Meg o, nel logio LiLo 
17 11°756 35°208 3°708 21'215 7°590 O°195 1°630 
18 11°936 35°524 Per 3°698 21°425 7°615 0218 a 1°946 
19 12°126 35°851 gate 3°084 2I°ss7 7°917 o'241 ee 2273 
20 12°329 36°200 ee 3°675 21°770 07°767 0°264 ae 2°622 
21 12°652 36°757 ss 3°648 22°141 7°889 0°322 ets 3°179 
22 11930 35°521 34°537 3°700 =. 20°574. 8087-0341 ee 1°943 


23 12126) §«635°851 9 35°264 83°684 20°714 8124 o°312 O°105 2°273 


The least accurate feature of the integrations lies in the treatment of the term 
de/dt, appearing in the equations for the structure of the star (cf. 19566). The 
estimation of this term for any integration requires a difference to be taken between 
the quantity « for the current integration and the value for the previous solution. 
Thus, for example, the value of de/dt for solution 20 depends on the difference 
between the value of ¢ for Solution 20 itself and the value for Solution 19. Nowin 
operating a machine of limited capacity (the present integrations were performed 
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on EDSAC 1) it is not possible to store the previous solution in any great detail, a 
circumstance that of necessity introduces some error. The extent of the error 
can be seen from Fig. 1, which shows that the points scatter around a mean line, 
the extent of the scatter being about +0-1 mag. Evidently an error of this order is 
adequately small. 

The concentration of metals was taken to be so low that in all cases there was no 
appreciablecontribution fromthe metals totheelectron pressureat the photosphere. 


-3-5 
3.4 
+ 3-3 
3-2 
3.1 
+3.0 
2.9 
.2-8 
- 2.7 
- 2.6 


2-5 


—logio LILO 


2.4 
L2.3 
+ 2.2 
2-8 
- 2-0 
L 1.9 
+ 1-8 
ri.7 
L 1.6 


- 1.5 








i i i I. L L will 4 


3-71 3.70 3.69 3.68 3.67 3-66 345 3.64 





"logis Ts 


Fic. 1. 


2. Discussion of the integrations.—The outstanding feature of the table is the 
marked change from Solution 21 to Solution 22. Material at the centre of the 
star is degenerate in Solution 21, the density being 6-00 x 10°gm/cm*. But no 
solution with a degenerate helium zone could be found that extended appreciably 
further than Solution 21. The reason for our failure to find such a solution can 
readily be understood from the following argument. 

It will be observed that the central temperature 7, rises steadily from 
3°89 x 10 deg. K at Solution 17 to 7-75 x 107 deg. K at Solution 21. This increase 
is due to a temperature gradient that maintains an outward flow of energy, the 
energy supplied by a general gravitational compression of the helium zone. Now 
any appreciable extension of the evolution beyond Solution 21 demands a central 
temperature above 8 x 10’ deg. K, when the energy yielded by nuclear reactions in 
helium exceeds the effect of compression. ‘his leads to an accentuated rise of the 
central temperature, which in turn increases the rate of nuclear energy generation 
still further...and soon. Evidently a rapid evolution of the star is to be expected 
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in which the helium zone becomes expanded to lower density. This expansion is 
necessary to reduce the opacity sufficiently to allow a more ready outflow of the 
energy generated in the helium. 

This is just what is shown by Solution 22. The central density in this solution 
is reduced to 4:05 x 104gm/cm', and the helium is no longer degenerate. 

Although it would have been very desirable to trace in detail all stages of the 
rapid evolution from Solution 21 to Solution 22, this was not possible in EDSAC I, 
since a prohibitively lengthy calculation would have been necessary. Our 
expedient was to “‘ pick up ”’ simply by trial and error the evolution at a stage after 
the expansion of the core had taken place. One important assumption was 
involved in this procedure: that there was no mixing of the inner helium with the 
outer hydrogen during the period of rapid evolution. Since there are no incipient 
signs of any such mixing at Solution 21, this seems an entirely reasonable 
hypothesis. 

The fraction of the total luminosity arising from energy production in helium 
increases appreciably between Solutions 22 and 23. An attempt to follow the 
evolution beyond Solution 23 showed that L,,, increased rapidly to a value close 
to L, and that a further instability in the structure of the star then ensued. Since 
in this case there were some signs that a mixing of hydrogen and helium might be 
important, no attempt was made to pick up the evolution beyond the second 
instability. More detailed computations with a substantially faster machine will 
be needed to follow the evolution still further. 

3. Comparison with observation.—All the solutions given in ‘lable I lie on 
what is usually called the giant branch of the representation in the Hertzsprung- 
Russell diagram. ‘The relation between L and 7 \,, for Solutions 17 to 21 agrees 
very closely with earlier calculations by Hoyle and Schwarzschild (1955). The 
values of M,,,/M for these solutions are about 20 per cent less than those obtained 
by Hoyle and Schwarszchild, however. 

Solutions 19 and 23 fall at almost exactly the same point of the Hertzsprung- 
Russell diagram, showing that apparently there is no unique interior structure for 
stars on the giant branch with L/L > 100. 

Provided that a unique structure is assumed at lower luminosities, in particular 
at Solution 17, we can compare the value of M,,,/M obtained for this solution with 
that deduced by Sandage (1957) from the observed luminosity function of the 
stars of the globular cluster M3. Making a slight adjustment forthe initial hydro- 
gen concentration 0-9309 used in the present calculations, and interpolating 
for logy, L/L@=1°630 between Sandage’s values, we obtain the following 
comparison :— 

My,./M at logy L/Lo = 1-630 


Present paper O°195 
Sandage’s semi-empirical value 0°209 





A second very close agreement with Sandage’s results follows from estimating the 
amount AM of the envelope that is added to the helium zone between 
logy) L/Lo = 1-630 and log, L/L ag =2°273. According to Sandage, AM/M=o0:121. 
From the way that this estimate was made, it must be interpreted as the total 
burning of the envelope material that takes place within the specified range of 
luminosity. Now from Table | it will be seen that burning occurs at two distinct 
stages. The first stage is given by the range from Solution 17 to Solution 19, and 
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the second stage by the range from Solution 22 to Solution 23. Writing AM,, 
AM, for the amounts of the envelope material that is burned in these two stages, 
we have from the table :— 


AM, AM, 


ar =0'046, aT =0°076. 
Our calculations therefore give 
AM _ AM,+AM, 
— = + = 0122, 


M M 


in almost exact agreement with Sandage’s estimate. It would clearly be very 
desirable to extend the present comparison to higher luminosities, but this un- 
fortunately requires an extension of the calculations beyond Solution 23, and such 
an extension is not yet available. 

Solution 21 gives the highest value of L. If we interpret this solution as giving 
the top of the giant branch, the magnitude of the brightest Type II stars becomes 
—2°8. This value follows from the equation 


M= + 4°62 —2°5 log,),Z/Lq + bolometric correction, 


where the bolometric correction at log 7',,,= 3°648 was taken as 0°53 mag. ‘This 
upper limit of — 2-8 is in excellent agreement with extensive data obtained by Arp 
(private communication) for a number of globular clusters. 

It has already been stated that the present calculations were carried out on the 
basis that electrons derived from metal atoms do not make an appreciable contri- 
bution to the electron pressure at the photosphere. ‘This case has a special 
importance in relation to the search that is being made for stars of abnormally low 
metal content. The effect of a small metal content to is push the top of the giant 
branch towards the right of the Hertzsprung-Russell diagram (Hoyle and 
Schwarzschild, loc. cit.). Hence if a giant sequence can be observed that falls 
near the sequence of solutions given in the above table (a sequence not pushed to 
the right in the H—R diagram), it is possible to infer that the stars in question have 
no more than a negligible trace of metals. This criterion of the presence or absence 
of metals may well prove to be more sensitive than a direct spectroscopic criterion. 

It is a pleasure to acknowledge our gratitude to the Director of the Mathe- 
matical Laboratory, Cambridge for provision of the facilities that allowed these 
calculations to be made. 


Department of Mathematics, St. John’s College, 
The University, Manchester 13: Cambridge. 
1958 March. 
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STABILITY OF POLYTROPIC GAS SPHERES 


W. B. Bonnor 
(Received 1958 May 23) 


Summary 


It has previously been shown that an isothermal gas sphere is gravitation- 
ally unstable to small perturbations in the pressure at sufficient distance from 
the centre. This work is here extended to polytropic gas spheres: it is 
shown that this type of instability does not affect spheres with indices n given 
by 0o<n<3, but that for »>3, a sphere larger than a certain critical size is 
unstable, and tends to collapse towards the centre. 





1. Introduction.—It has recently been shown (1, 3, 5) that the isothermal gas 
sphere exhibits an unexpected type of instability. To understand how this arises 
one may think of a fixed, very large mass of gas, at uniform temperature 7, subject 
to an experiment investigating Boyle’slaw. When the gas is in equilibrium with a 
very large volume V it satisfies the law in the form 

pV =NAT, (1.1) 
N being the number of molecules present, k being Boltzmann’s constant, and p 
the external pressure applied to maintain equilibrium. If p is gradually increased 
so that the gas passes through a series of equilibrium states, it is necessary to add to 
(1.1) a term representing the gravitational attraction of the molecules of gas on 
each other. Finally, at a certain critical pressure and volume (dependent on N 
and 7), it is found that dp/dV becomes positive, and the mass of gas unstable, witha 
tendency to collapse towards the centre. 

A result of this phenomenon is that a sufficiently large isothermal gas sphere is 
unstable, and it is possible that the instability may be of importance in explaining 
the formation of stars of Population I (5). It therefore seems of interest to find 
out whether a similar type of instability applies to polytropic gas spheres, which are 
used as stellar models, and this question is investigated here. 

It is known that certain complete polytropic gas spheres must exhibit an 
instability predicted by the virial theorem, but it has been shown (1, 5) that the 
conditions for instability according to the virial theorem are not necessarily the 
same as those for the instability described above, although the two phenomena are 
of the same general type. It does indeed turn out (Section 3) that for polytropic 
spheres the criteria for the two instabilities are different; in fact, the newly- 
discovered phenomenon gives criteria for the stability of incomplete polytropes, to 
which the virial theorem in its ordinary form does not apply. 

2. Perturbations in polytropic gas spheres.—A polytropic gas sphere is called 
complete if it is surrounded by vacuum, and if the polytropic relationship applies, 
with the same index , throughout its volume. _ It is called incomplete if the relation 
applies (with uniform m) within a sphere which is subject at its boundary to an 
equilibrium pressure provided by some other material. The following work will 
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apply to spheres which are complete or incomplete in the foregoing sense. One 
can also envisage annular regions of gas in polytropic equilbrium, but these are not 
dealt with here. 
Consider a spherical mass of gas in equilibrium, and let p(r) and p(r) be the 

density and pressure at radiusr. ‘The hydrostatic equation is 

dp __ pGM(r) 

> = —— ’ (2. I ) 
where M(r) is the mass inside radius 7, and G is the constant of gravitation. The 
polytropic relation is 


p=Kpatun) , (2.2) 
where x and mare constants. Using (2.2), and also 
M(r)=47 | " gr dr, (2.3) 
0 


we find that (2.1) gives 


I d pit In i 47G 
a5 \e s ar P )} = =" (2-4) 


Egn. (2.4) determines p(r) up to two arbitrary constants. 
Let us introduce new variables 6 and € by the formulae 


p=Aée", E=r/a (2.5) 
where 
k(1+n) 742 
“ | seas | é (2.6) 
A being a constant; then (2.4) gives 
1d dé 


This is Emden’s equation of index n. The constant A does not appear in it, and 
so may be taken as arbitrary. The solutions of (2.7) with 


6=1, — at €=0 


df 


correspond to polytropic spheres, of central density A. 

Now consider a sphere S, of radius ry inside and concentric with the polytropic 
sphere, and containing mass M. Let us contemplate a small variation in ry such 
that 

(i) M remains constant ; 
(ii) eqn. (2.1) is satisfied throughout S, after the variation, and at the new 
boundary of S, the pressure is that appropriate to equilibrium ; 

(iii) eqn. (2.2) is satisfied during the variation. 

We may imagine that the variation takes place in the course of fluctuations of 
the molecules of the polytrope or, if the latter is incomplete and of radius 74, from 
fluctuations in the material outside. 

From (2.3), (2.5) and (2.6), we find 


&o 
M = 4nd | &20” dé, 
which gives, with the use of (2.7), : 
M = — 4nad€,? (2), , (2.8) 
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Here, and elsewhere, a suffix o means that the value of the quantity concerned is 
to be taken on the boundary of Sy. Now put 





“= BA-G—Uny2 (2.9) 
1/2 
so that B= ow , (2.10) 
then (2.8) becomes 
M = — 4rf3\-A-3im2 € (Z) : (2.11) 
dé} 


As M is to remain constant as r, (or &)) varies we have 
d dé 
=$M = —arG2A-G-sim2) © | ge @ . 
+ ampen-an—wma(s — 3) ,7 ( z), dA. (2.12) 


The change in A allows for the slight alteration in central density which occurs 
during the variation. From (2.5) and (2.9) we have 


af)=8 (2). * = 519+ “(1 - = ) a+ (2.13) 


Substituting (2.13) into (2.12) and using (2.7) and (2.11), we find after a calculation 


_ fi—3/n)M | (i/n—1)% 
n={ + - hoa. 


871A y2p9 (2-14) 





From (2.2) we have 


i. “iF (1+ iy + 3: (2.15) 


te (B) 


so that after a calculation, in which (2.8) and (2.13) are used 

M Bi (1—1/n)M 
4rrAUn Bay 2 9 Bar B2x, +n) 
Substituting this into (2.15) we find with the help of (2.5) and (2.14) 


I Spo 3 (: - 3 lin Q2 > 2 
Po 3%o [(:- 3) m+ 1) 477po%o° | 47 Po!" Br, 


Further 


80) = — 





=3ant (14 :) B% pyi ttn y4 — (1-+n)(1- :) M2. (2.16) 
Simplifying this further, and using (2.2) and (2.10) we eventually find 
(n—3)G M? 


(2.17) 





(;F) en (n+ 1) Po io (n + 1) 82poro' 

dV)» 3(m—1)Vo_  n-3 M 
N—1 3p) Vo 

If we put G=o and take the density as uniform, we find 


(35), a t2se. 
aV n V,’ 


which is the formula obtained at once by differentiating (2.2) and using pV = const. 
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Formulae (2.16) and (2.17) are those required for the discussion of stability 
in the next Section. If in (2.17) we let n> oo, we get a formula for the isothermal 
gas sphere which agrees with that in (1). 

3. Stability.—It is easily shown that the polytropic gas sphere will be unstable 
to small variations in py and V, if (0p/0V) is positive. It follows at once from (2.17) 
that (dp/dV), is negative for all r, if 


I<n<}3. 


This means that a polytropic sphere with an index in this range is stable to the 
type of perturbation described in Section 2. This result is true whether or not the 
sphere is a complete polytrope. 
The case in which 
O<n<I (3-1) 


is best dealt with by considering (2.16). With this range of m the right-hand side 
of (2.16) is positive, and (dp/dV), will be positive if the term in square brackets on 
the left is positive, that is, if 


_3 - 3 
(: 3) m+ (2 t) amour >o, 


that is, if 47poro? (1—n) > (3-2) M. (3-2) 
If we write (3.2) in the form 
I-n  p 
an? x’ (3-3) 


p being the mean density within radius rp, we see that (dp/dV)) cannot be positive, 
with n satisfying (3.1), if the density diminishes as the radius increases. Because 
if the latter condition is fulfilled the right-hand side of (3.3) is greater than unity. 
whereas the left-hand side is less than unity. We conclude that no physically 
interesting case of instability arises if m satisfies (3.1). 

If we ignore, as being of no physical interest, spheres for which n<o, there 
remain those with 

n>3- 


If these are complete it follows from an application of the virial thoerem (2) that 
they are unstable. It was emphasized in the case of the isothermal gas sphere 
that this theorem applies only if the gas is surrounded by vacuum, and if this 
is not so the situation may be quite different. The work in Section 2 holds for a 
polytrope which is incomplete provided that it is surrounded by gas at the appro- 
priate equilibrium pressure. We can use this work to find how the stability of 
incomplete polytropes depends on the radius for > 3. 

To do this we need to study how the sign of (dp/dV), in (2.17) depends on ry. 
We can do this most easily by transforming to the variables and 0, and using 
Emden’s calculations (4). For > 3 it is clear from (2.17) that (dp/dV), has the 
sign of 











n—3 GM? 

~ n+ 1 Bxporot 
a M ’ 
N—T1 3pV, 
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and by using (2.2), (2.5), (2.6) and (2.8), we can transform this expression into 


— 9 dé n+) 
‘3 (3), Mo 
NR—-3, 19 -,/4 

r+ 234,10, (2), 
For given n, this expression can be worked out as a function of & by the use of 
Emden’s tables. One finds that for small &, (i.e. small 7)) (3.4) is negative, but 
when &, is increased to a certain value, é, (crit.), the numerator changes sign. For 
this value of &, (¢p/2V), becomes positive and the sphere becomes unstable. By 
arguments such as those used for the isothermal gas sphere, one can easily show 
that for all larger radii the polytropic sphere must be unstable, and that the 
instability must result in a tendency to collapse towards the centre. 

Table I shows (in the second column) the values of &, (crit.) for polytropes 
with various values of n (> 3). These have been calculated from Emden’s tables 
by interpolation, and are given to the nearest 1/10 of a unit of £. The appropriate 
values of r) can be obtained from (2.5) and (2.6). In the third column of Table 
I, €, indicates the greatest possible value of &, that is, the value of € for which 
the pressure vanishes. 


; (3-4) 








TABLE I 
n &, (crit.) & 
4 3°5 14°999 
4°5 3°2 32°143 
4°9 30 169°4678 
5 30 6) 
6 27 fo) 


It follows that incomplete polytropes with radii less than those referred to in 
Table I are stable to the type of perturbation considered in this paper. (One 
cannot of course assert from this that they would be stable to other perturbations. ) 
It is clear from the work of McCrea (5, p. 563 et seq.) that the instability envisaged 
is of the same general type as that predicted by the virial theorem, but the method 
used here enables one to take account of the presence of a surrounding medium. 
One may therefore take this work as supplementing the well-known criterion 
obtained from the virial theorem that complete polytropes are unstable if n > 3, and 
as showing that the criterion cannot be extended to incomplete polytropes of sizes 
less than that given by &, (crit.). 


Queen Elizabeth College 
(University of London), 
Campden Hill Road, 
London, W.8: 
1958 May 22. 


References 


(x) Bonnor, W. B., M.N., 116, 351, 1956. 

(2) Chandrasekhar, S., An Introduction to the Study of Stellar Structure, Chicago, 1939, Pp. 52. 
(3) Ebert, R., Z. Astrophys., 37, 217, 1955. 

(4) Emden, R., Gaskugeln, Leipzig and Berlin, 1907, 

(5) McCrea, W. H., M.N., 117, 562, 1957. 





528 Errata . Vol. 118 


ERRATA 


M.N., 113, No. 4, pp. 455-467, 1953: 

J. G. Porter and D. H. Sadler, The Accurate Calculation of Apparent Places 
of Stars. 

The statement “ sec 6 may be replaced approximately by tan 5” in para g (c) 
on p. 460 applies only to magnitude and not to sign. ‘The substance and con- 
clusions of the paper are unaffected by this error; but the formulae for the tabula- 
tion method of correction (c) apply only for northern (positive) declinations. 

For southern (negative) declinations: (1—sin 6) in the coefficients of the 
residual second-order terms No. 12 should be replaced by (1+sin 6); and the 
signs of the terms 2gh in the definitions of J, J’ in equations (14), (15) on p. 460 
should be negative. Separate tabulations of the values of J, J’ are therefore 
required for northern and for southern declinations. 

We are grateful to Dr ‘IT. Lederle of the Astronomisches Rechen-Institut, 
Heidelberg, for pointing out this error to us. 

J. G. P. 
1958 June 25. D. H. S. 


M.N., 117, No. 6, p. 581, 1957: 
List of Fellows elected: 
for Matthew Leslie Billis read Matthew Leslie Bellis. 
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PE3IOME JIOKJIAJIOB, B MEPEBOJLE HA PYCCKHM A3bIK 


MPEHMYINECTBO HEZABHCHMOM ACTPOJIABUH 
JIA OCHOBHOM ACTPOHOMHUM 


JIexyua umenn J>Kopmoxa Jlappuna, mpountaHHan mpodeccopom Axnyps J,aH>KOHOM 9 Tro Max 
1958 r. 


MSYUEHHE OVIPEDEJIEHHA ATMOC@PEPHOTO COTEPXKAHUA COJIHLIIA 
JI. Mazeacmoyn 


Vsyuaetca onpeyesieHve COJIHCUHbIX aTMOC(PepHbIX cCOeprxKaHHit C MOMOLIbIO HaGsIOJaeMbIx 
MHTCHCHBHOCTeli cmaObix JiMHHii CpayHrodepa. Mero BecoBbix (byHKUMii, Nomb3yroumiica 
pesyJIbTaTamMH HaOsOaeMOro MOTeCMHeCHHA Kpas JMCKa, pasOupaetcA B Aetanax. IIpeynaraetca 
ONOMIHUTeIbHbIM MeTO—MeTOs rpagventa IInaHkka, KOTOpbI He HyKaeTcH B HaOsNOJeHHAX 
NOTeMHeHHM K KpaxiM JMCKa, HO BCelleIO 3aBHCHT OT (Pu3H4UeCKHX XapaKTePHCTHK paccmaTpH- 
BaeMOi MOJeIM COHeUHOM aTMocdepprl. 

B cratbe jelaeTcA CCbIJIKa Ha CJlyuaii a30Ta (MOTeCHUMaJI MOHM3allMu = 14,54 BOJIbTa), P1IaBHbIM 
o6pa30M JJIH CpaBHeHHA pe3yJIbTaTOB, MOJIYYCHHbIX OOOHMHM BbIMIeyKasaHHbIMH MeTOamMH. B 
CTaTbe paCCMaTPHBaIOTCA JICBATb JIMHHM MOTIOUeHHA M WaloTCA ONpeyesieHHA COsep»xKaHHA a3s0Ta 
No Kako u3 HHx. OOHapyxKeHO 3HaUHTeJIBHOeC pacxoOrKeHHe M@X*KYy pesyJIbTaTaMM NMo1yueH- 
HbIMH OOOMMH MeTOAamMH. M3 TeopeTHyuecKoro pa360pa CJlejlyeT BbIBOM, UTO B CIyYae a3z0Ta, 
a HW BOOOUIe JIA BCeX 9JIEMCHTOB BbICOKOrO NOTCHUMasla HOHM3AalMH, MeTOJ, BeCOBbIX dbyHKuHi 
COBEPUICHHO HeHaéKeH HM HaO MOUb3OBaTBCH MeTOJIOM rpayMeHta IInaHKka. 

PaccmaTpuBaloTcaA TPH aTMOC*epHbIX MO/(CJIM HU, WIA a30Ta, BLIBOAATCA CpeqHHe pesysILTAaThI, 
NONyUeHHbIe NOCIeEHHM MeTOZOM (IpHHHmMaaA log An=12,00): mMogenb Knaaca—log An=7,71, 
MOJe1b Butensu—log Ax =7,81, Mozenb Cyuxapta—log Ax=7,93 (41d CpaBHeHHaA: mo Xionaeptcy 
9,02, a no Yucénpyy 8,61). 


TEMIIEPATYPbI BOSBYX]EHHA JIA CN-H3JIYYEHHA HYOKHEN 
XPOMOC®PEPBI 


Jl. B. Tomac 


MuxpodotometpupopaHa oOsIacTh CHcTemMbI TlosIoc CN mp 3883 A, ABsIAIollelica 3amMeTHO 
OCOGeHHOCTHIO CIICKTPOrpaMMbI BCIIbILIKH, MosyueHHoH PeymaHoM MpH 3aTMeHHH B 1952 Pr. 
Tlonocpr CN cCHJIbHO HCKaKeHbI COBMAaJ[CHHAMH, HO OKa3aJIOCb JJOCTaTOUHO HCHCKa>KeHHbIX JIMHHH, 
YTOObI MOXKHO ObINIO H3MeEPHTL pacnpeseseHve poTal[HOHHOM MHTeHCHBHOCTH B MosIoce O—O Ha 
Tpex BbICOTax HWKE ~ 400 KM. CaMonorsolujeHie 3HauHTesIbHO, B OCOOeHHOCTH y rpaHHiy MosI0C 
HM MpousBeyeHa OWeHKa ero BIMAHHA. W3mMepeHuA yKa3bIBaloT Ha OTPHLaTeJIbHbI rpayHeHT 
TeMMepaTypbI BO3OyxKLeHHA; TEMIepaTypa Ha BbICOTe B ~ 400 KM OKa3bIBaeTCA ~ 4500", a BEJIMUHHA 
ee rpagMeHta HeompeyereHHa. Merny tem CN-n3s1yueHie MOXKeT OKa3aTbCA B COOTBETCTBHK C 
NOMOXKUTEIBHBIM [pas[MeHTOM KHHeTHUeCKOM TeMMepaTypbI B MocuesHel MOee xpomocdepbI no 
ne Arepy, mpu ycuIoBuM yueTa OTKIIOHCHHA OT JIOKAaIbHOrO TePMOJIHHAMHYCCKOTO PaBHOBeCHA HM 
HECOBEPLICHHOCTH pe3yIbTaTOB HaGsIOJ(eHHH. 


CIIEKTP PEKOMBHHALINH BOJJOPOIA 
A, Bépoowecc 


JlaHo pellieHie KacKafHbIX yPOBHCHHH paMaTHBHOrO 3axBaTa JI BOJOpOsla, NPHHHMaA pu 
9TOM BO BHHMaHHe BbIPOXKIGHHOCTh OPONTaIbHOrTO yrIOBOTO MOMCHTa aTOMHbIX yPpOBHeil. 
PaccmoTpeHbI Bce ypoBHH (7, /) c n< 13, TakKxKe KaK M KOHTHHYYM. 

IIpumensvetcA MeTO, MepeBasla [VIA MOMYYeHHA aACHMIITOTMUCCKHX pa3JIO*KeHH J{MMOJIbHbIX 
MATPHUYHDIX IIEMCHTOB BOLOPOMa, ABHO BKINOYUAIOMIMX KBaHTOBOe YHCIIO /, a IPM MOMOLLM {{pyrHx 
METOJIOB BbIBCCHbI MOJIeSHbIc MpHONM3HMTeIbHbIe BhIpaxkenus. IIpuBeyenbr TaOmuUbI Ceuennii 
doToHOHH3allMH BOOpOa JIA BCeX ypoBHeil (7, /) cn<13, a TalkoKe KOod(pHUNeHTOB peKOMOnNHalMH 
JH Temlepatyp dIeKTpoHOB T,=10* uw 2.10! °K. 

JJaHbi pelieHuA KackajHbIx ypaBHeHHii {IA OOoMx curyuaes A uw B npn T,=10* n2. 10° K. 
Hacenenne OTACJIBHBIX YpoBHeil / JIA JaHHOLO n OTHIOJIb He MPOMOpUMOHAaIbHO HX CTaTHCTHUeCKOMy 
Becy. Ilo aHHbIM pacueTaM NosyualoTcA {1a OOoux cryuaes A u B GasbMeposckue {eKpemeHTHI, 
NONOGHbIe J[eKPeCMeHTaM, BbIBe/CHHbIM bokepom Hu Meutesnem ja cayuan B. 

Ha6soyaempie GasIbMepOBCKHe JICKPeCMCHTHI WICCTH M1aHeTaPHbIX TyMaHHOCTel, C NompaBKow 
Ha TIOKpaCHeHHe, XOPOWO COBMAaalOT C BbIUMCJICHHbIMH J|CKPeMeCHTaMH. 
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MHTEP®PEPOMETPHYUECKHE H3MEPEHMA JUIMH BOJIH 
VI. H3MEPEHHA B OBJIACTH 3940-4450 A 


C. Huxoasc u C. B. M. Karoo 


OnucaHo M3MepeHHe 17 COJIHCUHbIX M BaKYYM-AYPOBbIX {VIMH BOJIH B OOJIaCTH 3940-4450 A. 
OGujee 4NCIO COJIHEUHbIX KPaCHbIX CMeII[CHHii, H3Me€PCHHbIX MeTOJOM KPYIFOBbIX KaHaJIOB, 
cBeyeHO K 51. OTa rpynma WHHHii MHTepecHa NoTOMy, UTO OHA COMePeKHT 6 JIMHHH >KeTe3a H3 
OWHOrO MyIbTHIIIeTa: a*F—y*F°, pesoHaHCHyt0 JIMHH1O 2KesIe3a: 4375,9 A HW OTHOCHTeJIbHO CHJIBHYIO 
JIMHMIO B xpoMocidbepe, Mo mpuunHe Tit: 4443,8A. Kpacnbie cMelljeHHA STHX JIMHHH He 
corslacyloTCd HH C OJHHM H3 CYIUeCTBYIOIUMX MpegnonoxKennit. Jia 51 HaOmoyaemMol JIMHHH 
oOujee BOSpacTaHMe COJIHeYHOrO KpaCcHOrO CMeIIIeHHA C MHTeCHCHBHOCTbIO JIMHHH ABJIAeTCH 
€HHCTBCHHOM HaOsOaeMOi 3aKOHOMePHOCTBIO. 


METOJ,L KBAHTCBOTO JIEDEKTA 
M. JIoic. Cumon 


B metoge KBaHTOBOrO ecbeKTa MPHMCHAIOTCA MHTePNOJMPOBAaHHble WIM IKCTPanoOsIMpOBaHHbIe 
KBaHTOBBIe J[e(eKTbI WIA ONpeesIeHHA ACHMIMTOTHUeCKHX BHJOB AaTOMHBIX BOJIHOBbIX (pyHKIUMit. 
QTOT MeTOX MO?KHO IIPHMCHATE MPH BbIUMCICHHAX BEPOATHOCTeli ATOMHBIX MepexOjIOB MH CeueHHii 
(oTOHOHH3al|HH, TIpH pacueTaX CTOJIKHOBeHHi Me7KAY OJICKTPOHAaMH HM MOHAMH, a TaKoKe pu 
pelieHHu mpobsiem TBepyoro cocTOAHHA. 

JlaHHast CTaTbA COfepoxKHT OO30p MmperxKHUX paOoT, KacalOUMXCA OCHOBHBIX wWelleii MeTOJIa H 
WaeT HECKOJIBKO HOBbIX PeC3YJIBTATOB, OTHOCAIIMXCH K COCTOAHHAM MONO7KHTeIbHOM 9HEPruH U 
K HOPMHPOBKe BOJIHOBbIX (PYHKUHii CBASAHHBIX COCTOAHHH. O6cy2kqeHbI HEKOTOPbIe MpHJI0OxKeHHA 
MeToyIa. 


SBESQQbI-THPAHTbI THIIA II 
K. B. Xaszeaepoye u PD. Xoitiae 
JlaHvbI HHTerpaljHH MO MOcweOBAaTeIbHOCTH 3Be3-ruraHToB Tuna II, npuBojsaume K cTaqHH 


IBOJHOWUMH, B KOTOPOM BLIXO] SHEP 3 AMePHbIX peakYHii MoO resHio 3HauNTeeH. PesybTaTHl 
9THX MHTerpallMii BAAHMO XOporloO cormacyloTCcA Cc HaOTIOeHHAMH. 


YCTONMUMBOCTE MOJINTPOIMNMUECKHX LWAPOB FASA 
B. B. Bonnop 


‘ 


Panee 6bis10 JJOKa3aHO, YTO H3OTEPMHUeCKHH Wap rasa ABJIAeCTCH PpaBUTal|MOHHO HeyCTOMYMBbIM 
IIPH MaJIbIX BOSMYU[CHHAX aBJIeHHA Ha JOCTaTOUHO GOObLIOM paCCTOAHHH OT WeHTpa. OTH 
BbIBOJIbI MpHJlaraloTcA 3/{eCb K NOJMTPOMMYeCKHM LWiapaM rasa: OKa3bIBaeTCA, UTO 9Ta HeycTo;- 
YUMBOCTS HE BJIMACT Ha Wapbl, MOKa3saTeJIb M KOTOPbIX yJOBJIeETBOPAeT YCOBHIO O<n < 3, HO WIA 
n>3 M pagmyca Wapa, MpeBocxosALlWero ONpeyesIeHHYIO KPHTHUECKYIO BeIHUHHY, Wap HeycToOiuHB 
M CKJIOHEH K CIJIKOWIMBaHHIO 110 HalipaBJIeHHtO K LeHTpy. 
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